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1  Summary 


Large  spacecraft  such  as  the  space  shuttle  perturb  the  natural  environment  by  releasing  con¬ 
taminants  into  the  atmosphere  through  outgassing,  chemical  dumps  and  thruster  firings.  A  fluid 
model  has  been  developed  to  describe  the  large  scale  structure  of  the  contaminant  plasma  cloud 
The  model  is  two  dimensional  with  the  structure  along  the  magnetic  field  being  averaged  out.  The 
model  predicts  that  an  important  parameter  is  the  parallel  length  scale  of  the  cloud  since  this 
determines  the  effective  length  over  which  electrons  can  be  drawn  to  short  out  the  ions  flowing 
perpendicular  to  the  magnetic  field.  In  addition,  a  significant  amount  of  work  has  been  done  on 
the  small  scale  kinetic  processes  that  occur  in  the  contaminant  cloud  and  that  cannot  be  explained 
by  the  fluid  model.  These  small  scale  processes  include  instabilities  in  the  cloud,  ion  acceleration 
by  the  cloud  and  anomalous  ionization  in  the  cloud. 

2  Personnel  Involved  and  Outline 

This  work  has  involved  thesis  research  by  four  graduate  students.  Mr.  N.  Gatsonis  produced  the 
two  dimensional  fluid  model  of  the  plasma  cloud.  His  work  resulted  in  two  papers  in  the  Journal  of 
Geophysical  Research  which  are  included  as  appendices  to  this  report.  In  addition  he  received  the 
S.  M.  degree  on  this  topic.  Mr.  T.  Mogstad  worked  on  a  theory  of  low  frequency  instabilities  for 
the  plasma  cloud.  He  obtained  the  S.  M.  degree  for  his  thesis  which  was  entitled  “  Stability  Thec.y 
of  Plasma  Clouds  about  Large  Space  Vehicles  ”.  Currently,  Mr.  R.  Biasca  is  working  on  a  doctoral 
dissertation  on  the  possibilities  of  anomalous  ionization  in  the  contaminant  plasma  cloud  and  Mr. 
D  Rivas  is  working  on  a  doctoral  dissertation  exploring  other  kinetic  effects  in  the  contaminant 
plasma  cloud.  These  include  a  theory  of  ion  acceleration  and  high  frequency  instabilities  in  the 
cloud.  The  aim  is  to  explain  some  of  the  experimental  observations  which  'annot  be  explained 
by  a  fluid  theory.  In  addition,  Mr.  N.  Gatsonis  is  working  on  a  doctoral  dissertation  on  a  three 
dimensional  theory  of  the  plasma  cloud.  Finally,  Dr.  Michael  Gerver  ■•pent  some  time  working  on 


1 


this  prgram  and  contributed  to  the  work  of  Mr.  David  Rivas. 

The  outline  of  the  report  is  as  follows: 

1)  Kinetic  plasma  mechanisms  in  contaminant  plasma  clouds  (D.  Rivas  and  M.  Gerver) 

2)  Anomalous  ionization  in  contaminant  plasma  clouds  (R.  Biasca) 

3)  Low  frequency  instabilities  in  contaminant  plasma  clouds  (T.  Mogstad) 

4)  Three  dimensional  fluid  theory  of  contaminant  plasma  clouds  (N.  Gatsonis) 

5)  Two  dimensional  fluid  theory  of  contaminant  plasma  clouds  (N.  Gatsonis) 

3  Kinetic  Plasma  Mechanisms  in  Contaminant  Plasma  Clouds 

3.1  Introduction 

In  low  earth  orbit  the  space  shuttle  releases  neutral  water  or  ions  through  outgassing,  water 
dumps  or  thruster  firings  and  thus  perturbs  the  ambient  ionosphere.  Neutrals  and  ions  undergo 
chemical  reactions  and  form  a  contaminant  plasma  cloud  within  the  perturbed  environment  of  the 
shuttle.  Various  authors  (  Hastings  et  al.,  1989;  Eccles  et  al.,  1988)  have  studied  the  composition 
and  evolution  of  this  cloud  by  using  two-dimentional  fluid  models.  These  studies,  however,  did  not 
adequately  model  the  small  scale-length  structure  of  the  cloud.  In  this  report  we  therefore  extend 
the  theory  of  the  shuttle’s  plasma  cloud  by  studying  some  of  the  kinetic  processes  that  have  been 
observed  to  occur  in  orbit. 

We  are  mainly  interested  in  explaining  observations  of  high  energy  ambient  0+  ions  (Paterson 
et  al.,  1989;  Pickett  et  al.,  1985;  Shawhan  et  al.,1984)  that  were  detected  in  the  cloud  by  the 
Plasma  Diagnostic  Pa-kage  (PDP)  flown  aboard  the  third  Space  Shuttle  flight  (STS-3).  In  these 
experiments  O +  ions  with  energies  up  to  10’s  of  eV  were  observed  in  the  cloud  and  up  to  100  eV 
during  primary  and  vernier  thruster  firings.  Since  these  populations  of  ions  have  energies  at  least 
one  order  of  magnitude  larger  than  the  expected  5  eV  ram  kinetic  energy,  they  could  significately 
affect  the  structure  of  the  cloud.  The  mechanism  that  accelerates  these  ions  is  described  in  the 
section  that  follows. 
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3.2  Fermi-like  Ion  Acceleration  Mechanism 


The  acceleration  of  charged  particles  to  very  high  energies  is  also  a  problem  of  interest  in  the 
study  both  of  cosmic  rays  and  hot  plasmas.  Fermi  has  proposed  a  method  for  accelerating  particles 
to  cosmic  ray  energies.  In  the  proposed  mechanism,  a  charged  particle  is  moving  in  a  magnetic  field 
between  two  interstellar  clouds.  If  the  magnetic  field  in  the  clouds  is  assumed  to  be  greater  than  in 
the  intervening  region,  the  particle  is  trapped  between  two  magnetic  mirrors.  Such  trapping  will 
occur  of  course,  only  for  particles  whose  velocity  is  inclined  to  the  magnetic  field  at  a  sufficient 
angle.  The  clouds  comprising  the  two  mirrors  are  assumed  to  be  moving  toward  each  other  at  some 
relative  velocity.  A  charged  particle  now  gains  energy  on  each  reflection  from  the  mirror.  This 
acceleration  may  be  computed  from  the  constancy  of  the  longitudinal  adiabatic  invariant  (second 
adiabatic  invariant). 

In  our  work  we  explain  the  mechanism  that  accelerates  the  0+  ions  to  high  energies  by  extending 
this  Fermi  acceleration  concept  to  the  case  of  the  shuttle’s  plasma  cloud.  We  show  that  the  ambient 

ions  can  be  accelerated  by  interacting  with  the  irregular  electrostatic  potential  structure  of  the 
cloud.  Some  of  the  electrostatic  potential  irregularities  of  interest  to  us  are  those  due  to  the  partial 
shielding  of  the  motional  electric  field  V XB  by  the  dense  plasma  cloud  and  those  resulting  from 
the  electrostatic  noise,  instabilities,  electron  density  enhancements,  and  pressure  increases  observed 
in  such  events  as  flash  evaporator  system  releases  and  primary  and  vernier  thruster  firings  (Pickett 
et  al.,  1985;  Shawhan  et  al.,  1984). 

The  mechanism  that  fermi-accelerates  the  ambient  ions  can  be  described  as  follows.  Due  to  the 
orbital  motion,  from  a  frame  of  reference  moving  with  the  shuttle,  the  ambient  0+  ions  are  seen 
to  flow  with  a  5  eV  ram  kinetic  energy.  Those  ions  in  front  of  the  cloud  will  go  through  the  cloud 
and  interact  with  its  irregular  potential  structure.  In  cloud  regions  with  large  potential  gradients, 
specifically  with  potential  increases  greater  than  5  eV,  the  ions  will  not  have  enough  kinetic  energy 
to  overcome  the  potential  ’barrier’  and  will  therefore  bounce  off  it.  Now,  if  we  examine  this  event 
from  a  stationary  frame  (fixed  to  the  ground),  we  will  see  that  in  this  frame  the  ion  energy  increases 
Figure  1  illustrates  this  situation,  it  can  be  seen  that  this  situation  is  analogous  to  the  case  of  a 
ping-pong  ball,  which  is  initially  at  rest,  being  overun  by  a  brick  wall  (or  barrier)  moving  at  the 
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orbital  speed.  By  linear  momentum  conservation  it  can  be  shown  that  after  the  first  bounce  the 
ion  will  gain  a  velocity  equal  to  two  times  the  orbital  velocity  ( V ,).  Moreover,  further  increases 
in  kinetic  energy  will  be  possible  since  after  the  first  bounce  the  magnetic  field  will  deflect  the  ion 
trajectory  and  make  it  bounce  off  other  potential  barriers.  As  illustrated  in  Fig.  1,  in  each  bounce 
the  ion  velocity  will  be  incremented  by  as  much  as  2V,  until  the  ion  kinetic  energy  will  be  larger 
than  the  potential  barrier. 

3.3  Kinetic  Theory 

A  rigorous  theoretical  analysis  of  the  processes  described  in  the  previous  section  would  require 
the  derivation  of  a  kinetic  model  of  the  ambient  0+  ions  interacting  with  the  plasma  cloud.  One  of 
our  major  tasks  would  then  be  to  determine  how  the  electrostatic  potential  structure  of  the  cloud 
(which  to  zeroth  order  is  already  known  from  fluid  models  (Hastings  et  al.,  1989))  is  affected  by 
the  accelerated  ion  population.  Specifically,  we  would  be  interested  in  determining  whether  the 
electrostatic  potential  gradients  due  to  the  shielding  of  the  motional  electric  field  are  indeed  the 
driving  forces  that  accelerate  the  ambient  0+  ions.  If  so,  we  would  then  calculate  the  resulting 
self-consistent  structure  of  the  electrostatic  potential  with  high  energy  ion  populations  taken  into 
account. 

The  derivation  of  a  ‘closed  form’  kinetic  model  of  the  cloud,  however,  is  a  rather  difficult  if 
not  impossible  task.  A  reasonable  kinetic  model  of  the  self-consistent  shielded  potential  structure 
of  the  cloud  would  require  consideration  of  collisions  and  charge  exchange  reactions  between  the 
various  species  (four  neutral  and  three  ion  species)  and  this  would  complicate  the  mathematics 
enormously.  It  is  therefore  useful  to  borrow  some  of  the  tricks  and  short  cuts  that  are  often  used 
in  solving  other  problems  similar  to  ours.  In  the  next  section  we  shall  review  some  of  the  standard 
techniques  that  are  used  in  the  derivation  of  kinetic  models  for  strong  field  gradient  problems 
where  the  gradients  are  in  the  magnetic  field  and  collisions  can  be  neglected.  Then  we  shall  apply 
some  of  these  techniques  to  our  problem  and  we  will  derive  an  approximate  constant  of  motion  for 
particle  motion  in  strong  electrostatic  potential  gradient  regions  which  will  be  useful  to  analyze  the 
processes  that  occur  there. 
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3.4  Kinetic  Solution  to  Strong  Field  Gradient  Problems. 


There  is  a  great  deal  of  interest  at  present  in  acceleration  processes  in  large  gradients  in  the 
magnetic  field  which  usually  occur  in  critical  regions  such  as  in  shock  structures,  in  the  magne¬ 
topause  or  in  other  sheet  regions  such  as  in  the  magnetotail.  These  regions  are  important  because 
they  are  often  bounderies  separating  different  plasmas,  and  also  because  they  are  frequently  the 
location  of  particle  acceleration. 

The  standard  technique  used  in  solving  these  strong  B-field  gradient  problems  self-consistently 
is  the  integration  of  the  Vlasov  equation  in  terms  of  the  single  particle  constants  of  the  motion.  In 
these  problems  we  therefore  seek  solutions  of  the  Vlasov  equations  for  ions  and  electrons,  which  are 
expressible  in  terms  of  the  available  constants  of  motion,  and  are  suitable  to  represent  the  particle 
distributions  in  the  strong  field  gradient  configuration.  It  should  be  pointed  out  that  this  is  not 
the  same  as  providing  general  integrals  to  the  Vlasov  equation.  For  instance,  in  the  case  of  the 
two-dimenitonal  Vlasov  equation,  the  general  integral  of  this  equation  in  five  independent  variables 
(two  spatial  and  three  velocity  components)  is  expected  to  contain  an  arbitrary  function  of  four 
independent  arguments,  not  of  two  only  (which  is  the  number  of  available  constants  of  motion). 
Therefore,  these  standard  techniques  are  restricted  to  the  study  of  a  subclass  of  possible  solutions, 
in  the  hope,  of  course,  that  it  is  relevant  to  the  actual  problem. 

In  the  case  of  the  magnetotail,  many  authors  (Grad,  1961;  Harris,  1962;  Bird  and  Beard,  1972) 
showed  how  solutions  could  be  obtained  by  using  the  Poisson  and  Ampere  equations  after  expressing 
the  plasma  distribution  functions  in  terms  of  the  constants  of  the  motion.  In  the  adiabatic  limit 
the  constants  of  motion  in  the  magnetotail  are  the  particles  total  energy  and  canonical  momentum 
component  perpendicular  to  the  field  gradient.  In  two-dimensional  field  gradients,  however,  the 
canonical  momentum  is  no  longer  constant.  Nevertheless,  Propp  and  Beard  (1984)  and  Beard  and 
Cowley  (1985)  used  the  same  techniques  to  calculate  ion  orbits  in  a  two-dimentional  model  of  the 
tail.  In  the  latter  two  papers  the  authors  found  that  adiabatic  theory  predictions  were  generally 
valid  in  spite  of  ‘violation  of  the  adiabatic  limit’.  On  the  other  hand,  Roggers  and  Whipple  (1988) 
were  able  to  derive  a  more  accurate  model  of  the  two-dimentional  structure  of  the  magnetotail 
current  sheet  by  deriving  a  new  constant  of  the  motion,  the  generalized  first  adiabatic  invariant. 
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which  is  an  extension  of  the  magnetic  moment  for  regions  of  large  magnetic  field  gradients. 

Finally,  going  back  to  our  problem,  we  should  recall  that  in  the  analysis  of  the  acceleration  of  O  ' 
ions  in  the  shuttle’s  cloud  we  are  interested  in  strong  gradients  in  the  electrostatic  potential  (w:th 
the  magnetic  field  constant  everywhere).  Since  any  constant  of  the  motion  that  can  be  found  for  the 
particle  motion  in  such  regions  would  be  very  useful  in  analyzing  the  physical  processes  that  occur 
there,  we  shall  derive  a  generalized  first  adiabatic  invariant  for  the  case  of  the  cloud  by  using  the 
same  method  that  Roggers  and  Whipple  (1988)  have  used  in  their  analysis  of  the  geomagnetic  tail. 
The  procedure  for  obtaining  generalized  adiabatic  invariants  was  invented  by  M.  Kruskal  (1962) 
and  we  shall  therefore  denote  it  as  the  Kruskal  method.  Latter  on,  we  shall  integrate  numerically 
for  the  0+  ion  orbits  and  for  the  electrons  we  shall  use  the  generalized  adiabatic  invariant  along 
with  energy  conservation  to  predict  and  model  the  electron  behaviour. 


3.5  The  Generalized  Adiabatic  Invariant  for  Electrons  in  Regions  of  Strong 
Field  Gradients. 


3.5.1  The  Kruskal  Method 


Before  describing  the  Kruskal  method  let  us  first  of  all  recall  from  classical  mechanics  that 
whenever  a  system  has  a  periodic  motion  the  action  integral  f  pdq  (taken  over  a  period)  is  a 
constant,  where  p  is  the  generalized  momentum  and  q  is  the  generalized  coordinate.  If  a  slow 
change  is  made  in  the  system  so  that  the  motion  is  not  quite  periodic,  the  constant  in  the  motion 
does  not  change  and  is  then  called  the  adiabatic  invariant. 

The  Kruskal  method  (M.  Kruskal,  1962)  is  a  systematic  method  for  finding  asymptotic  solutions 
of  (and)  any  adiabatic  invariant  associated  with  coupled  differential  equations  having  the  form 


(1) 


where  X  is  a  N-dimentional  vector  and  e  is  a  small  parameter.  Furthermore,  ^  =  /(.Y,0)  has 
closed  (periodic)  solutions  that  are  closed  in  X-space. 

The  difficulty  in  treating  eq.(l)  is  that  /  is  in  general  both  large  and  (through  its  dependence 
on  .Y)  rapidly  varying  in  time. 
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The  Kruskal  method  is  a  transformation  procedure  for  splitting  off  the  rapid  time  variation  of 
X  into  a  single  angle  variable  4>{X)  leaving  an  (N-l)-dimentional  vector  Z(X)  which  varies  slowly 
By  construction,  Z  and  4>  have  equations  of  the  form 

~  =  eh(Z,e) 
d<t> 

m 

The  important  feature  of  the  transformation  is  that  4>  does  not  appear  in  h  or  u,  so  that  Z  varies 
smoothly  with  time  (with  a  time  derivative  of  the  order  of  the  small  parameter  e)  rather  than 
showing  oscillations  at  the  frequency  of  the  <p  variable. 

Equations  (2)  are  the  new  equations  of  motion,  equivalent  to  Eq.(l),  and  thus  the  solutions  of 
Eqs.(2)  when  converted  back  to  X-space  are  solutions  of  Eq.(l).  There  exist  closed  curves  (‘rings') 
in  X-space  on  which  Z  is  constant.  The  4>  variable  is  periodic  about  the  ring  and  specifies  position 
on  it  (phase  angle).  The  real  trajectory  differs  slightly  from  a  ring  and  it  is  not  closed  in  X-space 
If  equations  (1)  are  canonical  equations,  derived  from  a  hamiltonean,  then  there  exists  an  adiabatic 
invariant  in  addition  to  the  asymptotic  solutions  (in  e)  of  the  set  (l).  The  invariant  is  an  integral 
over  a  ring,  not  over  any  segment  of  an  actual  particle  trajectory: 

J  ■=  p{Z,4>)  —  Zj^).d<t>  (3) 

J  ring  0<p 

p  being  the  canonical  momentum  and  q  the  generalized  coordinate.  For  an  adiabatic  invariant  to 
exist,  it  is  not  necessary  that  eq.(l)  itself  be  of  canonical  form  only  that  it  be  possible  to  convert 
it  to  that  form. 

The  transformation  X  to  ( Z,4> )  is  actually  made  by  way  of  intermediate  variables  (V',  i/),  where 
Y  is  constant  on  the  closed  unperturbed  (i.e.  for  e  =  0)  curve  (’loop’)  and  i/  is  some  phase  angle 
about  the  loop.  The  distinction  among  actual  trajectories,  loops,  and  rings  can  be  further  classified 
as  follows:  an  actual  trajectory  (i.e.  for  e  £  0)  is  generally  not  even  a  closed  path  in  the  N- 
dimentional  space  of  X.  In  the  special  case  (unperturbed  case)  where  e  —  0,  the  trajectory  in  A'  is 
required  to  be  closed  in  order  to  apply  the  Kruskal  theory  Only  in  this  case  a  loop  roinsidrs  with 
a  ring 
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3.5.2  Derivation  of  the  Generalized  Adiabatic  Invariant  for  the  Electrons. 

Let  us  now  use  the  Kruskal  method  to  derive  a  generalized  adiabatic  invariant  for  regions  with 
strong  electrostatic  potential  gradients.  We  define  our  coordinate  system  to  be  moving  with  the 
shuttle,  with  the  z-axis  pointing  towards  the  wake,  the  y-axis  pointing  towards  the  motional  electric 
field  direction  VXB  and  the  z-axis  pointing  rwards  the  B-field  direction. 

In  order  to  appl,  generalized  adiabatic  theory,  the  motion  .1  electric  field  V  XP  nust  be  a&_urned 
to  be  small  and  will  be  denoted  by  e.  Since  the  dimentionless  quantity  which  oust  be  small  is  fh 
ratio  (s  /  B)  drift  speed  (which  in  our  case  is  equal  to  the  orbital  velocity,  ..e.  r  Km/sec)  to  the 
particle  thermal  velocity  {vthe  =  132.  Km/sec,  »  >„  =  0.8  Km/sec  ),  this  derivation  will  ^e  valid 
only  for  the  electrons. 

We  define  <I>(z)  as  the  electrostatic  potential  end  assume  that  it  only  changes  strongly  wi.h  z. 
This  potential  does  not  include  the  potent. a  e  3tiib'’t'on  due  to  the  motional  el°ctric  field  VXB 
).  The  equations  of  motion  then  are 


dv.  e  _  t  __  ,  , 

if  =  m V*(l)  (4> 

di>V  e  '  0\ 

iT  =  --B| 

where  the  charge  e  equals  to  -e  for  the  electrons.  These  equations  are  of  the  form  of  eq.(l)  with 
X  —  (z,vx,vv).  When  e  =  0,  the  particle  does  not  drift  gradually  accross  the  shock  in  and 

all  the  three  components  of  X  are  periodic  as  required. 

The  next  step  in  the  Kruskal  procedure  is  to  use  the  constants  of  the  zeroth-order  equations 
(i.e.  eqs(l)  with  e  =  0)  as  starting  points  for  constructing  the  ’Z-variables’  that  are  useful  when 
£  is  finite.  In  our  case,  these  constants  are  the  generalized  transverse  energy  W  and  the  canonical 


y-rnomentum  p 


W  —  K  4-  2 e$/m 
p  =  pv  =  rnvv  eAv 


where  K  -  u*  and  Av  =  xB.  The  quantities  W  and  p  are  called  ’T-variables’  in  Kruskal’s 
terminology ;  it  can  be  shown  that  these  variable  '>ave  time  derivatives  of  order  e.  The  phase  angle 
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v  may  be  defined  by 


where  <  v„  >t  is  the  gyro  averaged  y-component  of  velocity,  and  u  is  an  appropiatly  defined 
phase  angle  for  the  gyromotion.  The  advantage  of  the  Z-variables  is  that  they  are  strictly  drift 
parameters,  varying  slowly  and  smoothly  as  the  particles  drift  in  the  i-direction.  It  can  be  shown 
that  the  time-derivatives  of  the  Z-variables  are  non-oscillatory  and  of  order  e.  Moreover,  to  zeroth 
order  in  e  the  phase  angle  <j>  is  equal  to  u. 

The  particle’s  motion  accross  the  layer  can  be  described  in  terms  of  the  drift  parameters  and  the 
generalized  invariant.  The  generalized  invariant  is  constant  throughout  a  particle’s  orbit  through 
first  order  in  e  even  in  the  region  with  the  strong  gradient  in  the  electrostatic  potential.  In  our 
case  the  action  integral  eq.(3)  in  terms  of  the  Z-variables  reduces  to 

J  —  i  pzdx 

J  nng 

In  order  to  estimate  this  constant  of  the  motion  for  a  very  strong  potential  gradient,  which  we 
shall  model  as  a  step  function  in  the  potential,  we  first  define  Q  =  2eA where  AO  is  the 
potential  jump.  We  must  distinguish  particles  that  reflect  from  the  discontinuity  from  those  that 
cross.  Reflection  occurs  when  sin2 1/0  <  Q  for  eAO  >  0  (a  repulsive  potential  jump).  Here,  u  is  the 
phase  angle  tan_1(-ur/vv)  as  mentioned  above,  and  i/o  is  the  value  of  u  as  the  particle  is  incident  on 
the  discontinuity  from  the  upstream  side.  The  angle  v  jumps  discontinously  for  both  reflected  and 
transmitted  particles  because  vz  changes  discontinously  as  a  consequence  of  the  infinite  acceleration 
at  the  potential  jump.  The  invariant  is  given  by 

m2W  , 

J  =  — —  i/q  -  sin  i/o  Cv>s  , 

eB 
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for  reflected  particles,  and 


mJtV  ,  (^o  ~  sin  t'o  cos  vo) 

—7-I  Q 

(l-gjj.r-cos-Mp/mC^ci-g))1/2!} 

B 

cos  t/Q(sin2  vp  -  Q)1/2, 

B  1 

for  transmitted  particles.  The  reflected  particles  will  continue  to  ’bounce’  against  the  discontinuity, 
drifting  tangentially,  and  gradually  gain  energy  until  they  are  able  to  overcome  the  potential  barrier 
and  propagate  across  it. 

3.6  Acceleration  of  Ambient  Ions  by  the  Cloud  Shielding  of  the  Motional 
Electric  Field 

In  the  reference  frame  moving  with  the  shuttle  and  in  the  absence  of  any  electric  field  pertur¬ 
bation  an  observer  will  see  the  unshielded  motional  electric  field  Em  =  V  X  B.  However,  within 
a  plasma  cloud  density  perturbation  the  electric  field  and  consequently  the  drift  velocity  will  be 
determined  by  the  balance  between  Pederson  and  Hall  currents,  diamagnetic  currents,  "pickup” 
currents  as  well  as  parallel  currents  carried  by  the  AJfven  wings.  The  self  consistent  interaction  of 
these  currents  will  create  a  polarization  field  which  will  shield  the  imposed  motional  electric  field 
Consequently,  the  electric  field  within  the  plasma  cloud  will  be  |  E  |<|  Em  \  and  the  ions  will  drift 
backward  (in  the  shuttle  frame)  with  a  speed  which  is  less  than  the  orbital.  In  the  case  that  the 
motional  electric  field  is  entirely  shielded,  the  plasma  cloud  will  be  stationary  or  it  will  be  moving 
with  the  shuttle  velocity  in  a  fixed  frame. 

Caledonia  et  al.  (1987)  showed  that  experimental  measurements  of  the  contaminant  plasma 
cloud  were  consistent  with  an  ion  residence  time  in  the  vicinity  of  the  shuttle  of  40  ms.  Theoretical 
work  on  the  plasma  cloud  about  the  space  shuttle  started  with  the  suggestion  that  the  plasma  cloud 
could  be  highly  polarized  leading  to  a  long  residence  time  for  the  ions  in  the  shuttle  frame  (Katz 
et  al.,  1984).  A  model  for  the  unsteady  motion  of  the  plasma  cloud  was  developed  by  Hastings 
et  al.  (1988),  where  it  was  shown  that  for  the  densities  around  the  shuttle  the  polarization  of  the 
cloud  was  highly  dependent  upon  the  ion  density.  Recently,  Hastings  et  al.(1989)  have  extended 
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this  model  (1988)  to  include  the  coupling  of  the  cloud  to  Alfven  waves  which  can  carry  current 
parallel  to  the  magnetic  field.  Another  study  of  the  plasma  cloud  which  takes  into  account  the 
effect  of  Hall  and  Pederson  currents  was  undertaken  by  Eccles  et  al.  (1988).  The  Hall  currents  are 
shown  to  rotate  the  polarization  <  hat  it  is  not  antiparallel  with  the  motional  electric  field. 

In  the  following  sections  we  wii.  „  ine  whether  a  cloud  that  shields  the  motional  electric 
field  could  have  a  self  consistent  potential  structure  that  could  fermi  accelerate  the  ambient  0~  ions 
to  the  energies  observed  in  the  experiments.  We  will  also  make  an  effort  to  estimate  the  electrostatic 
effects  of  the  accelerated  particles;  therefore,  we  will  try  to  determine  a  cloud  potential  configuration 
which  will  be  self  consistent  with  all  the  cloud  and  ambient  species  including  the  accelerated  ions. 
Our  approach  to  the  problem  will  be  to  use  Hastings  et.  al.  (1989)  2-D  fluid  model  to  model 
the  cloud  and  neutral  ambient  species  and  to  treat  the  ambient  O +  ions  by  using  kinetic  theory 
techniques.  This  will  therefore  allows  us  to  take  Fermi  acceleration  effects  into  account.  In  the 
following  section  we  start  by  analysing  the  parts  of  the  2-D  fluid  model  relevant  to  our  problem. 

3.6.1  Two-dimensional  Fluid  Model  of  the  Cloud 

Hastings  and  Gatsonis  (1989)  developed  a  two-dimensional  fluid  model  for  the  motion  of  the 
contaminant  cloud  perpendicular  to  the  magnetic  field  lines.  They  assumed  that  the  contaminant 
cloud  consisted  of  ions  such  as  0  +  ,  HjO+ ,  and  neutrals  such  as  O,  H,  OH  and  H^O. 

Their  numerical  solution  examined  the  effects  of  Alfven  wave  coupling,  neutral  water  density,  ion 
temperature  and  initial  conditions  in  the  cloud  motion.  Their  results  show  that  in  low  density 
neutral  water  clouds  (<  109  cm-3)  the  shielding  of  the  electric  field  was  small.  On  the  other  hand, 
in  neutral  water  clouds  with  densities  of  interest  for  shuttle  conditions  (~  1010  cm-3)  the  shielding 
was  predicted  to  be  of  the  order  of  the  motional  electric  field  and  the  ratios  of  the  line  averaged 
densities  of  the  ions  was  consistent  with  experimental  data. 

Hastings  et  al.  (1989)  have  solved  the  fluid  equations  for  a  5  Km  by  5  Km  square  region  using  a 
uniform  101  by  101  cartesian  mesh.  It  is  useful  to  review  the  algorithm  that  they  have  used,  latter 
on  we  shall  modify  it  to  include  kinetic  effects.  This  algorithm  is  as  follows: 

1)  With  the  density  profiles  nWj0+(z,y),  nH:i0+(x,y)  and  n0+{x,y)  given  at  some  time,  inte- 
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grate  the  current  balance  equation 

V  J  =  0  (5) 

to  determine  the  potential  <f>{x,y). 

2)  Use  the  obtained  potential  <p(x ,  y)  in  the  momentum  equations  to  determine  the  ion  velocities. 

3)  Use  the  obtained  ion  velocities  in  the  continuity  equations  to  advance  the  densities  by  one 
time  step.  Then,  go  back  to  (1). 

It  is  also  useful  to  examine  the  current  balance  equation  of  Hastings  et  til.  (1989).  This  will 
enable  us  to  rigorously  determine  the  physical  parameters  that  contribute  to  the  shielding  of  the 
motional  electric  field.  Due  to  the  small  collisionality  assumption  in  the  model  (i//f2  <S  1,  which 
restricts  the  model  to  water  densities  below  10u  cm’3),  the  electron  current  canceled  out  with 
the  zeroth  order  ion  current  (Hall  current)  and  thus  the  current  balance  was  made  between  the 
first  order  perpendicular  ion  current  and  the  parallel  current  carried  by  Alfven  wings.  The  current 
balance  equation,  which  in  Hastings  model  is  solved  for  <f>,  can  be  put  to  the  familiar  form: 

V  •  (eV<t>)  =  pezt  (6) 

where  e  and  pezt  are  independent  of  <p.  This  two  dimensional  equation  has  already  been  integrated 
along  the  Alfven  characteristics;  pezt  contains  terms  due  to  the  diamagnetic  ion  currents  and  to 
currents  that  arise  from  the  friction  forces  between  ion  and  neutral  species.  The  parameter  c  which 
can  be  thought  of  as  a  dielectric  constant  has  the  form 


where  v0  is  the  Alfven  velocity  and  P,  is  the  effective  ion-neutral  plus  reactive  collision  frequency. 
The  first  term  in  this  equation  represents  the  current  carried  by  Alfven  waves  and  the  second  the 
pederson  current. 

The  parallel  interaction  length  to  perform  the  parallel  integration  over  can  be  estimated  as 

cloud  (8) 

vL 

where  rL./ouj  is  the  radius  of  the  cloud  in  the  perpendicular  direction  and  is  the  average  perpen¬ 
dicular  drift  velocity  of  the  cloud.  This  expression  is  the  distance  an  Alfven  wave  can  travel  in  the 
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time  the  cloud  takes  to  convect  across  the  magnetic  field  line.  If  we  use  eq.(8)  in  eq.(7)  we  obtain 
the  dielectric  constant  as 

e  ~  1  -r  2 rciou<i —  (9) 

Vj_ 

For  plasmas  with  low  neutral  density  the  collision  frequency  is  small  and  the  drift  velocity  is  of  the 
order  of  the  orbital  speed.  In  this  case,  e  ~  1  and  the  electric  field  in  the  plasma  is  the  motional 
electric  field  since  the  plasma  does  not  shield  out  the  field.  Physically,  the  parallel  current  carried 
by  Alfven  waves  is  so  large  that  no  significant  field  can  built  up  in  the  plasma  cloud  to  shield  out 
the  imposed  motional  electric  field.  For  a  plasma  where  neutral  water  density  is  high,  the  collision 
frequency  is  large  and  the  drift  velocity  is  small.  We  find  that  f  1  and  so  the  electric  field  is 
highly  shielded.  In  this  case  the  perpendicular  ion  current  is  sufficiently  large  to  give  rise  to  a 
substantial  change  in  the  motional  electric  field. 

3.6.2  Fluid  Model  Results  Suggesting  the  Need  of  a  Kinetic  Model 

The  electrostatic  potential  configuration  obtained  by  Hastings  et  al.(1989)  2-D  fluid  model  for 
a  cloud  with  an  initial  water  density  of  l010cm~3  which  decreases  radially  away  from  the  shuttle  as 
(50  m./r)J  is  shown  in  Figure  2.  In  the  model  the  inertia  terms  in  the  ion  momentum  equations  were 
neglected  and  therefore  to  zeroth  order  in  (i//f 1)  the  "fluid”  ions  drifted  along  the  equipotential 
lines  at  the  E  x  B  drift  velocity.  A  closer  inspection  of  Figure  2b,  however,  suggests  that  it  is 
possible  that  the  0+  ions  could  not  really  drifted  along  the  equipotential  lines.  This  is  because, 
from  basic  guiding  center  theory,  we  know  that  the  only  time  that  the  ion  larmor  motion  can  be 
ignored  and  the  guiding  center  drift  be  described  in  terms  of  the  local  Ex  Bj  B2  velocity  (guiding 
centers  moving  along  the  equipotential  lines)  is  when  the  electric  field  gradients  are  small  in  a 
larmor  radius  distance.  In  the  case  of  the  cloud,  however,  the  ion  larmor  radius  at  the  orbital  speed 
(Vo/n 0+)  is  about  40  meters  and  in  this  distance,  referring  to  Hastings  results  (Fig.  2b),  we  see 
sizable  gradients  in  the  electric  field.  Therefore,  it  is  questionable  whether  in  Hastings  model  it 
was  right  to  assume  that  the  ion  trajectories  would  follow  the  equipotential  lines. 

Referring  to  Fig.  2b,  we  can  see  that  if  the  ambient  0+  ions  in  front  of  the  shuttle  structure 
were  not  to  follow  the  equipotential  lines  but  instead  they  were  to  follow  straight  trajectories,  by 
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the  time  they  would  have  completed  one  gyroperiod  (0.035  sec.  in  low  earth  orbit)  they  would  have 
drifted  250  meters  at  the  orbital  velocity  and  seen  potential  increases  of  as  much  as  35  eV.  Although 
this  situation  is  impossible,  since  the  0+  ions  have  kinetic  energies  of  only  5  eV  and  therefore  do 
not  have  enough  energy  to  follow  straight  trajectories,  this  fictitious  experiment  suggests  that  if 
the  real  ion  trajectories  would  lie  somewhere  between  straight  lines  and  the  equipotential  lines  we 
would  expect  some  kind  of  bouncing  off  potential  "barrier”  or  Fermi  acceleration  to  occur. 

All  this  suggests  the  need  of  the  derivation  of  a  kinetic  model  of  the  the  cloud  where  ion  larmor 
motion  effects  and  thus  the  possibility  of  fermi  acceleration  would  be  taken  into  account.  If  fermi 
acceleration  were  to  occur  or  the  ion  trajectories  be  different  from  those  of  the  2-D  fluid  model 
this  analysis  would  help  us  to  determine  a  more  accurate  self-consistent  potential  configuration 
(self-consistent  with  the  cloud  and  ambient  species  including  the  accelerated  ions). 

3.6.3  Kinetic  Treatment  of  the  Cloud 

In  this  section  we  describe  an  analysis  of  the  cloud  which  extends  the  2-D  fluid  model  to  include 
the  kinetic  effects  of  the  ambient  ions.  As  we  discussed  in  the  last  section,  it  is  very  difficult  to 
attempt  to  derive  a  2-D  kinetic  model  of  the  cloud  that  would  include  the  shielding  of  the  motional 
electric  field.  Since  the  dominant  mechanism  for  shielding  are  the  ion-neutral  collisions  (see  sec 
5.1),  a  kinetic  model  that  would  be  useful  to  investigate  fermi  acceleration  due  to  this  shielding 
would  have  to  include  collisions  and  this  would  complicate  the  mathematics  enormously.  In  our 
discussion  of  the  previous  section,  however,  we  argued  the  possibility  of  0+  ion  fermi  acceleration 
in  a  shielded  cloud  and  that  a  kinetic  treatment  of  the  cloud  was  necessary  in  order  to  obtain 
accurate  self  consistent  models  of  the  potential  structure. 

We  have  therefore  found  it  convenient  to  use  the  2-D  fluid  model  of  Hastings  et  al.  (1989),  which 
accurately  models  ion-neutral  and  reactive  collisions  and  Alfven  currents  (crucial  for  modeling  the 
shielding),  to  model  the  cloud  ions  and  neutrals  and  ambient  neutrals  but  to  treat  the  ambient 
0+  ions  in  a  different  manner.  The  ions  are  treated  as  individual  particles  and  their  behaviour  is 
obtained  by  numerically  integrating  their  equation  of  motion.  Therefore,  ion  inertia  and  larmor 
motion  effects  are  now  taken  into  account.  These  kinetic  ambient  ions  are  then  merged  into  the 
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flu'd  cloud  model  as  follows:  First,  the  0+  equation  of  motion  is  numerically  integrated  for  various 
particles  and  the  divergence  of  the  0+  current  (VJ^-t-)  is  calculated  by  using  PIC-code  particle 
weighting  techniques.  Then,  the  fluid  value  of  VJ0+  in  Hastings  current  balance  equation  (eq 
5)  is  replaced  by  our  kinetic  value  of  V J<5+.  Finally,  the  potential  is  obtained  by  integrating  the 
resulting  current  balance  equation  which  now  includes  fermi  acceleration  effects  and  the  calculations 
are  advanced  in  time  by  using  the  algorithm  described  in  section  3.5.1. 

In  section  3.5.5  we  will  show  that  this  kinetic  analysis  should  lead  to  the  same  results  as 
those  of  the  fluid  analysis  (Hastings  et  al.  1989).  We  will  show  that  in  shuttle  cloud  models  ion 
larmor  motion  effects  can  be  ignored,  fermi  acceleration  cannot  occur  and  thus  fluid  models  give  an 
accurate  description  of  the  shielding  of  the  motional  electric  field  by  the  cloud.  In  the  next  section 
we  determine  the  conditions  necessary  for  shielding  the  motional  electric  field.  We  will  find  this 
analysis  to  be  useful  in  understanding  the  physics  of  various  plasma  clouds  which  we  will  test  in 
section  3.5.5  for  fermi  acceleration. 

3.6.4  Cloud  Effective  Radius  and  Densities  for  Shielding  the  Motional  Electric  Field 

In  order  to  determine  the  values  of  the  cloud  parameters  (ion  and  neutral  densities,  effective 
cloud  re'hus,  etc.)  that  are  needed  for  the  cloud  to  shield  the  motional  electric  field,  we  consider  a 
simple  cloud  model.  We  assume  the  cloud  to  be  spherical  in  shape  and  with  a  finite  effective  radius 
denoted  by  Rc.  Moreover,  inside  the  cloud  we  assume  the  ion  (n,c)  and  neutral  (nn7o)  densities  to 
be  uniform.  Therefore,  for  this  cloud  to  shield  the  motional  electric  field  the  following  conditions 
must  be  satisfied: 

1)  the  sheath  impedance  has  to  be  greater  than  the  Aifven  impedance 

2)  the  sheath  impedance  has  to  be  greater  than  the  resistance  across  the  cloud 

Here  we  are  referring  to  the  sheath  along  the  magnetic  field. 

If  (1)  is  violated,  the  parallel  current  carried  by  Aifven  waves  will  be  so  large  that  no  significant 
field  can  build  up  in  the  cloud  to  shield  the  motional  electric  field.  In  this  case  the  current  flowing 
through  the  cloud  will  be  limited  by  the  Aifven  impedance  rather  than  by  the  saturation  current 
of  the  ambient  plasma,  Aifven  wings  will  become  charged  up  to  the  same  potential  as  the  cloud 
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The  ions  will  therefore  flow  along  the  equipotential  lines. 

If  (2)  is  violated,  then  cloud  will  be  unable  to  short  out  the  motional  electric  field  in  the  ambient 
plasma  and  thus  the  cloud  ions  and  electrons  will  drift  downstream  of  the  structure  with  velocities 
near  the  orbital. 

Now,  the  Alfven  impedance  is  (in  cgs  units) 

=  C~1(uc,/uJfn)ambtent  (l0) 

Ci 

and  the  sheath  impedance  is  the  potential  across  the  cloud 

VqBqRc 

c 

divided  by  either  the  saturation  current  from  the  ambient  plasma  to  the  cloud,  or  the  saturation 
current  from  the  cloud  to  the  ambient  plasma,  whichever  is  less.  The  reason  being  that  the  cloud 
will  always  develop  a  potential  to  keep  the  net  current  into  the  cloud  equal  to  zero.  The  saturation 


current  into  the  cloud  is 


e(ri,  e* )  ambient  ^{^e^t)cloud  — - 


and  the  saturation  current  out  of  the  cloud  is 

e(n,  V,  )cloud  c{fteVc)amb\cnt 

The  latter  is  clearly  less  than  e(neve)cloud,  since,  as  we  shall  see,  n«  »  neo,  while  Tee  is  not  that 
much  less  than  Tea  and  thus  vea  a  vec.  Therefore,  the  sheath  impedance  is 


Vq  Bq  Rr  , 

c{2nRl)[ 


€(n,U,)c/ouj  -4“  e(  Vg  )  ambient ) 


^ -  (ii) 

cRce(nlcVtc  t"  ^ea^ea) 

Now,  since  the  largest  current  across  the  cloud  \s  the  Pederson  current,  the  resistance  across  the 


W1 ci  i  ^ci  V in  ^ 

?  n  1 

'^f,cRc  I'.n  Vc, 
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where  i is  the  ion  ( 0+ )  neutral  ( H?0 )  collision  frequency.  In  our  case, 

i>tn  —  nn7oa HiOvo+  where  oh^O  —  10~  15cm2  and  v0+  can  be  taken  as  the  orbital  velocity, 

thus  i ~  5  x  10_nn//3osec_1  (n//3o  can  range  anywhere  between  108  -  1014cm-3). 

The  first  criterion  is  then 


u>. 


Bo 


o(w p*  ) ambient  C Rc€^tlicVxc  4“  rleoVea) 

Using  c  =  3  x  1010  cm/sec,  u>c,  =  eBo/m,c  =  2  x  102/sec  for  Ot,  Bo  =  0.3  Gauss, 
(up.jamMent  =  (4*e2nea/m, ) ,  nea  ~  5  x  105  cm-3  (dayside),  vea  ~  107  cm/sec 
(for  Tta  -  0.1  eV),  V0  =  8x  105  cm/sec,  e  =  5x  10~10  statcoul, 
and  v,'e  =  5  x  104  cm/sec  for  HzO +  at  500°  K,  this  becomes 


Rc  < 


10 


13 


n,e  +  108 


(13) 


(14) 


The  second  criterion  is 


Vo  Bq 


UJ,.  ,  U!r 


cRce(n,cv,c  +  ntlJvea)  ^lxcRc  I'm  w, 


V\n  ^ 


(15) 


Using  =  47rntce2/m,  =  10s n,c,  it  can  be  shown  that  this  inequality  can  be  satisfied  only  for 
nn7o  >  I010cm~3(vxn  >  0.5sec-1),  for  n//3o  -  10ncm_3  this  becomes 


0.32 


0.08 


niC  -+-  2  x  107  nxe 


or 

nlc  >  5  x  106cm-3  ( 16) 

Finally,  by  combining  eqs.(  14)  and  ( 16)  we  can  see  that  the  two  criteria  for  shielding  will  be  satisfied 
for  n//jO  >  1010cm'3,nle  >  5  x  106cm-3,  and  Rc  <  1  Km. 

Before  comparing  these  results  with  those  from  2-D  fluid  models,  let  us  note  that  if  the  cloud  is 
elongated  by  a  factor  A  the  conductivity  across  the  cloud  would  also  increase  by  this  factor,  then 
from  eq.(12)  ntc  times  nn2o  can  he  reduced  by  A,  and  the  cloud  would  still  be  able  to  short  out  the 
electric  field,  up  to  a  factor  A  -  VAjV0  which  is  a  few  hundred  in  low  earth  orbit.  What  is  really 
relevant  is  the  scale  length,  along  the  magnetic  field,  of  times  n,e,  not  just  of  n,c  However, 

since  the  neutral  density  is  spreading  out  in  all  directions  it  is  questionable  wether  the  elongation 
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would  be  a  few  hundred.  Hastings  et  al.  (1989)  have  found  substantial  shorting  out  of  the  electric 
field  for  slightly  lower  densities  since  they  assumed  that  A  is  greater  than  VA/V0. 

It  should  be  mentioned  that  Tcc  will  be  close  to  T,c,  since  the  electron  cooling  time  by  coulomb 
collisions  with  ions,  ~  (i/eme/m,)~l  ~  6  x  10-2  sec,  is  shorter  than  the  escape  time  of  the  electrons 
from  the  cloud,  which  is  —  ~  10  1  sec;  or  at  least  they  are  comparable.  This  confirms  that 

n(acta  <  necvec,  so  the  cloud  will  develop  a  positive  potential  to  hold  electrons  in. 

3.6.5  Ion  Larmor  Motion  Effects  in  Clouds  that  Shield  the  Motional  Electric  Field 

In  this  section  we  determine  the  larmor  motion  effects  on  ambient  0+  ion  which  go  through 
various  types  of  clouds  that  shield  the  motional  electric  field.  We  do  this  by  comparing  the  fluid 
model  ion  trajectories  and  velocities  with  those  of  a  kinetic  model.  In  the  analysis  of  various  types 
of  clouds  we  consider  the  trajectory  that  is  most  likely  to  be  affected  by  larmor  motion  effects,  as 
we  shall  see,  this  trajectory  will  usually  be  that  for  ions  which  are  initially  traveling  towards  the 
center  of  the  cloud. 

Figs.  2  a  an  b  show,  respectively,  the  potential  configuration  and  equipotential  contours  of  a 
cloud  after  0.2  sec  of  being  released  with  an  initial  central  neutral  water  density  of  1010cm~3  which 
decreases  radially  as  (50  m/r)2  (Hastings  et  al.  1989).  Now,  the  dominant  ion  species  in  the  cloud  is 
the  ambient  0+  ion,  with  a  density  of  nQ+  —  2  x  105cm-3  which  is  almost  uniform  throughout  the 
cloud.  Therefore,  since  the  other  cloud  ion  species  (H20+  and  H$0+)  are  the  products  of  reactive 
collisions  between  0+  and  H^O,  they  will  also  have  densities  more  or  less  uniform  throughout 
the  cloud.  The  ion  densities  nWj0+.  ~  105cm-3  and  nH;s0+  ~  104cm~3  decrease  more  slowly 
than  (50  m./r)2,  the  water  density;  in  fact,  at  r=500  m  they  are  down  by  only  one  order  of 
magnitude  Consequently,  because  the  total  cloud  ion  density  (which  is  dominated  by  the  ambient 
ions)  decreases  radially  very  slowly,  we  can  see  that  the  resultant  potential  configuration  (Fig.  2a) 
is  smooth  and  thus  the  gradients  in  the  potential  are  small  (Fig.  2b)  even  though  the  motional 
electric  field  is  entirely  shielded  at  the  structure. 

Fig.  2c  shows  a  guiding  center  or  fluid  model  trajectory  (an  equipotential  line)  of  an  ion  which  is 
initially  in  front  of  the  shuttle.  Fig.  2d  shows  the  trajectory  of  tht.-  same  ion  obtained  by  numerically 
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integrating  the  equation  of  motion 

dV  , 

m0+  —  =  e(E  +  V  x  B)  (17) 

Therefore,  the  trajectory  in  Fig.  2d  includes  ion  larmor  motion,  the  ion  temperature  is  assumed  to 
be  0.1  eV.  Figs.  2e  and  2f  show  the  corresponding  velocities.  It  is  interesting  to  note  that  the  fluid 
(fig.  2c)  and  kinetic  (fig. 2d)  trajectories  coincide  (or  almost  coincide  since  V0  »  vt/,)  and  that  the 
ion  velocity  (or  ‘kinetic  velocity’)  (fig.  2f)  has  an  average  velocity  equal  to  the  fluid  E  x  B /  B2  drift 
(fig.  2e).  We  can  therefore  see  that  the  fluid  approximation  is  accurate  for  modeling  this  cloud 
since  the  potential  gradients  are  smooth  enough  so  that  ion  larmor  motion  can  be  disregarded. 
Obviously,  no  fermi  acceleration  occurs  in  this  case. 

In  Fig. 3  we  examine  the  case  of  a  cloud  with  larger  potential  gradients  due  to  the  shielding  of 
the  background  electric  field.  We  consider  a  very  small  cloud  with  a  central  neutral  density  same  as 
that  of  the  cloud  in  fig.  2,  1010cm-3,  but  which  decreases  much  faster,  as  (10  m./r)2.  Figs.  3c  and 
d  show,  respectively,  the  fluid  and  kinetic  trajectories  and  figs.  3e  and  f  show  the  corresponding 
fluid  and  kinetic  velocities  of  an  ambient  ion  which  is  initially  in  front  of  the  structure.  Although 
the  motional  electric  field  at  the  structure  is  almost  entirely  shielded  and  the  gradients  in  the 
potential  due  to  the  shielding  (potential  increases  as  large  as  20  eV  in  a  larmor  radius  distance, 
40  m.  at  the  orbital  velocity)  are  in  this  case  five  times  larger  than  those  in  fig.  2,  we  can  see 
that  the  ion  trajectory  (fig.  2d)  adjusts  itself  to  follow  the  equipotential  line  (fig.  2c).  Therefore, 
the  fluid  model  is  also  accurate  for  describing  this  cloud:  the  kinetic  trajectory  (Fig. 3d)  coincides 
with  the  equipotential  line  (fig. 3c)  and  the  kinetic  velocity  (fig.  3f)  has  an  average  velocity  equal 
to  the  local  E  x  B / B2  drift  velocity  (fig.  3e).  Nevertheless,  a  small  larmor  radius  effect  can  be 
noticed  in  fig.  3f,  in  the  wake  of  the  cloud  the  ion  larmor  speed  is  two  times  greater  than  upstream 
of  the  cloud  and  thus  the  ion  temperature  is  four  times  greater.  This,  however,  is  not  expected  to 
change  the  potential  configuration  significantly  since  the  fluid  quantities  (figs.  3  c  and  e)  are  still 
equal  to  the  average  of  the  kinetic  quantities  (figs.  3  d  and  f)  and  ion  neutral  collisions  will  tend  to 
keep  the  ion  temperature  down  (although  they  will  be  unable  to  keep  it  equal  to  the  background 
neutral  temperature  since  the  cooling  time  >  2Rc/Vq,  the  time  for  ions  to  travel  2 Rc).  It  is  also 
interesting  to  note  that  in  the  cloud  the  total  (drift  -^larmor)  ion  velocity  goes  up  to  15  Km  sec 
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with  corresponding  energies  of  15  eV.  This  could  therefore  explain  some  of  the  peaks  in  the  0~  ion 
energy  distribution  near  this  value  observed  by  Paterson  et  al.  (1989). 

Figs.  4  a  and  b  show,  respectively,  the  potential  configuration  and  equipotential  contours  of 
a  non  self  consistent  or  artificial  potential  of  a  cloud  with  an  outer  region  where  the  potential  is 
smooth  and  an  inner  region  which  exhibits  large  potential  gradients.  In  this  artificial  case  the 
cloud  not  only  shields  the  motional  electric  field  but,  from  the  reference  frame  of  the  earth,  the 
cloud  is  seen  to  develop  an  E-field  that  is  even  greater  in  magnitude  than  the  motional  electric 
field.  Numerical  studies  of  the  cloud  show  that  this  type  of  potential  configuration  (although  not 
as  exaggerated  as  in  fig.  4)  can  occur  in  some  particular  cases  where  Alfven  currents  are  taken  into 
account  (Gatsonis,  1989). 

Figs.  4  c  and  d  show,  respectively,  the  fluid  and  kinetic  trajectories  and  figs.  4  e  and  f  show 
the  corresponding  fluid  and  kinetic  velocities  of  the  ion  which  is  likely  to  experience  the  largest 
ion  larmor  motion  effects.  It  is  interesting  to  note  that  if  the  outer  region  of  the  cloud  (especially 
upstream  of  the  cloud)  has  a  smooth  potential  structure,  then  the  majority  of  the  upstream  ion, 
even  those  initially  in  front  of  the  structure,  will  drift  at  the  local  E  x  BjB2  velocity  along  the 
smooth  equipotential  lines  that  surround  the  inner  region.  Therefore,  only  very  few  ions  will  enter 
into  the  inner  cloud  region  and  be  fermi  accelerated  by  the  sharp  potential  gradients  there. 

The  above  analysis  of  the  three  different  clouds  (figs.  2,  3,  and  4)  shows  that  in  water  clouds 
with  central  neutral  densities  less  than  1012cm-3  the  dominant  ion  species  will  be  the  ambient 
O T  ion  with  a  fairly  uniform  density  distribution.  Furthermore,  since  the  other  ion  species  in 
the  cloud  are  the  products  of  reactive  collisions  between  the  0+  ions  and  water  neutrals,  the 
overall  cloud  ion  density  will  decrease  slowly  in  the  radial  direction  (more  slowly  than  the  neutral 
density).  Therefore,  although  near  the  shuttle  the  motional  electric  field  will  be  shielded,  the 
resulting  potential  structure  will  always  be  smooth  with  weak  potential  gradients  in  the  outer 
region  of  the  cloud.  Consequently,  the  ions  will  just  drift  according  to  the  local  E  x  BjB2  velocity 
and  no  fermi  acceleration  will  be  possible  from  the  shielding  of  the  clouds  with  nn2o  <  1012cm-3 

It  might  be  possible,  however,  that  in  a  small  dense  cloud,  with  a  central  water  density  of 
10ucm-3  which  decreases  radially  as  (10  m  ,/r)2,  the  H^O +  and  H$0+  ions  (products  of  reactive 
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collisions  between  0+  and  H2O)  could  outnumber  the  Or  ions  in  the  cloud  and  decrease  radially 
in  density  almost  as  fast  as  the  neutral  water  density.  This  in  turn  could  lead  to  large  potential 
gradients  in  the  outer  region  of  the  cloud  that  could  Fermi  accelerate  the  ambient  ions.  We  could 
not  consider  this  case  because  Hastings  et  al  (1989)  fluid  model  is  restricted  to  (nn2o)  central  <  1012 
cm"3. 

3.7  Acceleration  of  Ambient  O r  Ions  by  the  Plume  of  Primary  and  Vernier 
thrusters. 

Since  most  of  the  very  high  energy  ions  (~  100  eV)  observations  (Pickett  et  al.  1985)  were  during 
primary  and  vernier  thruster  firings,  in  this  section  we  examine  the  possibility  of  fermi  acceleration 
of  the  ambient  ions  in  the  plume  of  a  primary  thruster.  The  space  shuttle  fires  the  primary  thrusters 
in  pulses  of  80  ms.  These  thrusters  have  the  following  characteristics:  thrust  — 3925  N,  7,p  =  281  sec, 
c=2753  m/sec,  P=1  Mpa,  Nozzle  diameter^  60  cm,  mass  flow=  800  g/sec  of  LOr  —  150  g/sec 
of  LH2 ■  Now,  since  the  plume  exhaust  neutral  density  is  about  1017cm“3,  therefore  much  higher 
than  the  cloud  neutral  density,  large  number  of  exhaust  neutrals  could  be  ionized  due  to  reactive 
collisions  with  ambient  ions.  Then,  it  might  be  possible  that  the  resulting  high  plasma  density 
(~  109cm"3)  could  shield  the  motional  electric  field  inside  the  plume.  If  the  plume  axis  would  be 
directed  perpendicular  to  both  the  El-field  and  the  direction  of  orbital  motion,  sharp  gradients  in 
the  electrostatic  potential  of  the  cloud  could  exist  in  the  plume  region.  These  gradients  could  then 
fermi  accelerate  the  ambient  ions  even  for  low  water  cloud  densities.  Next,  we  determine  whether 
this  situation  could  happen. 

3.7.1  Analysis  of  the  Shielding  of  the  Motional  Electric  Field  in  the  Plume  of  the 
Primary  Thruster 

In  this  section  we  make  a  basic  calculation  of  the  shielding  of  the  motional  electric  field  E __  in 
the  primary  thruster  plume  to  determine  whether  the  density  in  the  exhaust  is  high  enough  to  short 
out  the  background  electric  field  significantly.  Assuming  that  the  axis  of  the  plume  is  perpendicular 
to  both  the  magnetic  field  and  the  direction  of  orbital  motion  and  neglecting  the  exhaust  velocity 
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(since  c  ~  ^0/4),  the  perpendicular  current  density  is 


J±  =  {^Tl 


Pg 


099  t-  o 

Tl03S^c 


-)Es 
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where  r\ott  %  10”  sec  is  the  time  an  ion  takes  to  go  out  of  the  plume,  and  the  ion  density  ir> 
the  plume  is  n,p  ss  2.5  x  109cm~3.  This  density  is  based  on  the  ionization  rate  from  collisions  of 
ambient  oxygen  atoms  with  plume  neutrals,  and  rlo,,.  The  electron  contribution  to  J __  'urns  on; 
to  be  much  less  than  the  ion  contribution,  which  is 


J  L  —  ^p,Tlc 


i>  E±.  —  ’’b pe(  Tlott) 

m, 


(19) 


Note  that  the  expression  in  parentheses  is  just  the  velocity  that  the  ion  gets  accelerated  to  in  a 
time  T[ott .  We  use  this  expression  for  Jj_i  rather  than  —  ?-£;  -  E>_  (the  pederson  current),  because 
Tio>s^ct  1-  We  then  find  that  J  ±  =  4  x  10-5  A/cm2.  By  contrast,  Jy  due  to  electrons  streaming 
out  of  the  plume  (in  the  absence  of  an  ambipolar  potential)  would  be  neeve  ss  10*  2  A/cm2,  which 
is  much  greater.  It  follows  that  V  ■  J  -  0  can  be  satisfied  only  by  having  Jy  as  0  everywhere,  i.e 
the  potential  of  the  plume  on  each  field  line  must  adjust  itself  so  that  electron  losses  on  that  field 
line  are  balanced  by  incoming  ambient  electrons  and  ion  losses.  This  is  in  contrast  to  the  case  of 
a  big  dense  plasma  cloud,  where  V  •  J=  0  is  satisfied  by  having  Jy  as  J  *_ ,  and  £N  is  shorted  out. 
If  we  assume  the  plume  electrons  to  be  Maxwellian,  the  plume  potential  will  then  have  to  be  a 
few  Tt  above  the  ambient  potential  everywhere.  This  is  to  keep  the  electrons  from  flowing  out  of 
the  plume  and  thus  dy  as  0  as  is  required.  Te  will  be  some  fraction  of  5  eV,  since  the  electrons 
are  produced  by  ionization  of  exhaust  atoms  and  incoming  ambient  oxygen  atoms  which  collide 
with  an  energy  of  about  5  eV.  Therefore,  the  plume  potential  will  be  on  the  order  of  5  V.  Whether 
or  not  incoming  ambient  0+  ions  can  reflect  from  the  plume  potential  depends  on  whether  it  is 
greater  than  or  less  than  5  V.  To  answer  this  depends  on  a  more  detailed  calculation,  but  at  least 
it  is  plausible  that  fermi  acceleration  (only  one  bounce)  occurs. 

However,  it  is  possible  that  the  plume  potential  be  much  larger  than  5eV  since  it  is  very 
sensitive  to  the  electron  temperature  in  the  plume,  and  for  that  matter  to  the  electron  distribution 
function.  Because,  it  is  not  obvious  that  the  electron  distribution  function  is  maxweilian,  even  a 
small  deviation  from  a  maxweilian,  e  g.  an  extended  tail,  can  have  a  big  effect  on  the  potential 


22 


which  will  adjust  itself  to  keep  the  electrons  from  flowing  out  of  the  plume,  if  the  potential  is  a  few 
times  greater  than  the  bulk  temperature.  This  will  depend  on  the  velocity  distribution  at  which 
electrons  are  born  (from  ionization  by  collisions  between  exhaust  molecules  and  incoming  ambient 
atoms,  and  any  other  ionization  source  that  is  important),  and  on  any  instabilities  that  occur  in  the 
initial  distribution,  and  how  these  affect  the  distribution.  Also,  ambient  electrons  coming  in  along 
the  magnetic  field,  and  accelerated  by  the  potential  may  contribute  significantly  to  the  electron 
distribution  and  instabilities. 

3.7.2  Trajectories  of  Ambient  Ions  across  a  Shielded  Primary  Thruster  Plume 

Figs.  5  a  and  b  show,  respectively,  the  potential  configuration  and  potential  contours  of  a 
plume  fired  in  the  direction  of  the  motional  electric  field  (negative  y-direction).  Figs.  6  a  and  b 
show  these  same  quantities  but  for  a  plume  which  is  fired  in  the  direction  opposite  to  the  motional 
electric  field.  These  potential  characteristics  are  (artificial)  non  self  consistent  (with  the  plume  and 
cloud  plasma)  and  the  plume  is  assumed  to  shield  the  motional  electric  field  at  a  potential  equal  to 
the  potential  at  the  structure.  Therefore,  to  the  left  of  the  shuttle  (in  the  negative  y-direction )the 
ambient  ions  will  see  a  potential  hill  at  the  plume  (fig.  5)  and  to  the  right  of  the  shuttle  they 
will  see  a  potential  well  (fig.  5).  For  simplicity,  we  have  studied  these  two  cases  separately  and 
only  considered  1-D  potential  shocks  (variations  in  the  plume  potential  along  the  y-direction  were 
ignored).  Figs  5c  and  6c  show  the  (fluid)  E  x  B  drift  0+  ion  trajectories  (which  is  of  course 
coincide  with  the  equipotential  lines)  and  figs.  5e  and  6e  show  the  corresponding  guiding  center 
drift  velocities  (equal  to  the  local  E  x  B/B2).  Figs.  5d  and  6d  show  the  particle  (or  kinetic  ) 
trajectories  of  the  ambient  Or  ions  (obtained  by  numerically  integrating  the  equation  of  motion) 
and  figs.  5f  and  6f  show  the  corresponding  particle  velocities  .  In  figs.  5  and  6  the  potential  shocks 
caused  by  the  plume  shielding  were  assumed  to  be  equal  in  magnitude  but  opposite  in  sign  and 
relatively  smooth. 

We  can  see  that  in  the  case  of  a  potential  well  (fig.  6)  the  ion  larmor  motion  does  not  affect 
the  trajectory  and  thus  a  fluid  description  of  the  ambient  ions  would  be  relatively  accurate  in  this 
case  The  small  increase  in  the  ion  temperature  downstream  of  the  plume  should  not  in  reality 
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be  that  high  due  to  ion-neutral  collisions  at  the  plume.  On  the  other  hand,  in  a  potential  shock 
equal  in  magnitude  but  opposite  in  sign  (a  potential  hill,  fig.  5)  we  can  see  that  ion  larmor  motion 
effects  become  significant.  Although  we  do  not  see  the  ion  bounce  off  the  shock  and  upstream  of 
the  shock  (up  to  the  plume  axis)  the  fluid  description  seems  to  be  accurate,  the  sharp  potential 
gradient  at  the  plume  axis  makes  the  ion  guiding  center  leave  the  equipotential  line  and  accelerate 
to  very  high  velocities. 

Finally,  in  fig.  7a  we  show  the  potential  of  a  plume  with  a  potential  hill  of  about  6  V  but 
with  large  potential  variations  in  distances  of  the  order  of  the  ion  larmor  radius  at  the  thermal 
velocity  (Ad>  as  6  V  in  distances  of  5  meters).  Figs.  7  b  and  c  show  the  ambient  ion  fluid  an  kinetic 
trajectories  and  figs.  7  d  and  e  show  the  corresponding  velocities.  In  this  plume  we  see  that  ion 
larmor  motion  and  ion  inertia  effects  are  important  since  we  see  the  ion  bouncing  off  the  plume 
(fig.  7  c)  and  thus  fermi  accelerating  to  velocities  three  times  larger  than  the  orbital  in  one  bounce. 
It  is  interesting  to  note  that  if  the  background  motional  electric  field  were  not  present  the  particle 
would  fermi  accelerate  only  to  velocities  twice  the  orbital  (fig.  1).  However,  in  our  case,  after  the 
first  bounce  the  background  motional  electric  field  does  work  on  the  ion  and  thus  the  ion  ends  up 
with  energies  of  more  than  6  times  its  initial  (not  only  4  times  as  predicted  by  fermi  acceleration). 
This  agrees  with  the  observed  peak  in  the  0+  energy  distribution  near  30  eV  (Paterson  et  ah, 
1989).  It  is  also  interesting  to  note  in  fig.  7e,  that  after  the  first  bounce,  due  to  the  work  of  the 
background  motional  electric  field  on  the  ion,  the  second  time  the  ion  encounters  the  shock  it  has 
a  velocity  of  2.3  times  V0  (corresponding  to  a  kinetic  energy  of  26  eV  ).  Nevertheless,  the  second 
time  the  ion  encounters  the  shock  it  could  still  bounce  off  it,  even  if  the  potential  increase  at  the 
shock  is  less  than  26  V,  as  long  as  the  ion  kinetic  energy  component  perpendicular  to  the  shock  is 
less  than  that  of  the  potential  increase  at  the  shock.  Therefore,  more  than  one  bounce  could  be 
possible  with  plume  shielding  potentials  not  much  larger  than  5  V  above  the  background  potential. 
If  the  ion  could  bounce  three  times  it  could  readily  be  accelerated  to  energies  of  about  100  eV’s. 
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3.8  Summary  and  Conclucions 

A  basic  analysis  of  the  shielding  of  the  motional  electric  field  by  a  spherical  plasma  cloud  shows 
that,  when  Alfven  currents  are  taken  into  account  as  well  as  the  sheath  along  the  R-field,  the 
motional  electric  field  will  be  shielded  for  nn7o  >  1010cm-3,  nlc  >  5  x  106cm-3,  and  Rc  <  1  Km 
( Rc  is  the  effective  radius  of  the  cloud  which  could  be  defined  as  the  distance  at  which  the  neutral 
water  density  is  three  orders  of  magnitude  lower  than  the  density  at  the  center).  If  the  cloud  is 
elongated  by  a  factor  A,  then  the  requirements  of  tih7o  times  n,c  decrease  by  this  factor  up  to 
VA/V o  which  is  a  few  hundred  in  low  earth  orbit. 

The  self  consistent  study  of  the  acceleration  of  ambient  ions  by  the  potential  of  a  polarized 
cloud  requires  a  kinetic  treatment.  Analytical  kinetic  models  pose  difficult  mathematical  problems 
since  collisions  must  be  taken  into  account.  These  models  are  also  complicated  by  the  fact  that 
when  considering  fermi  acceleration  problems,  regions  of  strong  field  gradients  must  be  included. 
Nevertheless,  it  is  possible  to  derive  a  constant  of  the  motion  of  the  electrons  in  these  regions  for 
studying  the  processes  that  occur  there. 

A  study  that  is  based  on  a  self  consistent  fluid  model  of  the  cloud  but  that  takes  into  account  O' 
ion  kinetic  effects  shows  that  water  clouds  with  central  neutral  water  densities  less  than  1012cm~3 
will  have  smooth  potential  distributions  for  all  values  of  /?«.,  the  effective  radius.  In  these  clouds 
the  dominant  ion  species  is  the  ambient  0+  with  a  density  distribution  that  is  more  or  less  uniform 
throughout  the  cloud.  Since  the  other  cloud  ion  species,  H20  +  and  H30* ,  are  products  of  reactive 
collisions  between  0+  and  the  cloud’s  H20  neutrals,  they  will  have  density  profiles  that  do  not  fall 
as  fast  as  the  n//3o  profiles  and  for  nn7o  <  1013cm-3  will  be  at  the  most  of  the  order  of  magnitude 
of  nQ  +  .  This  will  lead  to  smooth  potential  configurations  that  will  be  unable  to  fermi  accelerate 
the  ambient  ions. 

Large  and  medium  size  water  clouds  ( Rc  >  600  m)  with  (nH}0  jeentrai  >  1013cm“3  will  have 
potential  configurations  that  will  be  smooth  in  the  outer  region  of  the  cloud  for  the  same  reason  as 
the  less  dense  clouds.  Therefore,  no  fermi  acceleration  is  expected  to  occur  in  these  cases  either 
On  the  other  hand,  very  small  polarized  clouds  ( Rc  =;  100  m)  with  (nH7o)centrcU  >  1013cm~3  can 
have  potential  gradients  that  are  sharp  enough  in  the  outer  region  of  the  cloud  to  fermi  accelerate 
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the  ambient  ions. 

The  plume  of  a  primary  thruster  of  the  shuttle  has  a  neutral  density  of  1017cm~3  and  is  ionized 
by  reactive  collisions  with  the  ambient  atomic  oxygen  to  plume  ion  densities  of  109cm-3  Due  to 
the  fact  that  the  time  an  ion  takes  to  go  out  of  the  plume  is  small  (10~3  sec),  the  perpendicular 
current  in  the  plume  is  small.  Charge  conservation  therefore  requires  that  the  parallel  current  be 
also  small.  Consequently,  the  plume  potential  has  to  adjust  itself  to  hold  the  electrons  in  and  is  very 
sensitive  to  the  electron  distribution  function.  If  the  plume  electrons  were  Maxwellian  the  plume 
potential  would  be  about  5  V  above  the  background  potential  and  this  could  lead  to  one  bounce  or 
fermi  accelerating  the  0+  ions  to  26  eV.  However,  because  of  instabilities  and  ionization  processes 
that  take  place  in  the  plume  the  electron  distribution  function  is  not  expected  to  be  Maxwellian 
and  this  could  lead  to  much  larger  plume  potentials.  Therefore,  ambient  ions  could  bounce  two  or 
three  times  off  the  plume  and  easily  end  up  with  energies  of  about  100  eV.  (because  the  motional 
electric  field  also  does  work  on  the  ions  that  bounce  off). 

Finally,  it  is  important  to  mention  that  the  major  problem  in  the  self  consistent  study  of  fermi 
acceleration  lies  in  modeling  the  strong  potential  gradients  or  shocks  self  consistently.  If  these  shocks 
are  due  to  the  shielding  of  the  motional  electric  field,  as  in  the  case  of  polarized  clouds  and  primary 
thruster  plumes,  a  self  consistent  model  must  have  to  take  ion-neutral  collisions  into  account. 
Therefore,  since  conventional  PIC  codes  and  particle  tracking  codes  do  not  include  collisions,  they 
cannot  be  used  in  modeling  these  problems.  A  hybrid  treatment,  such  as  the  one  described  in  this 
report  (using  the  2-D  fluid  model  of  the  cloud  and  kinetic  0+  ions),  seems  to  be  a  good  way  to 
approach  the  problem  of  the  plume. 
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4  Anomalous  ionization  in  contaminant  plasma  clouds 

The  computational  simulation  of  medium  and  low  frequency  plasma  processes  can  aid  greatly  in 
developing  an  understanding  of  the  interactions  between  a  spacecraft  and  the  plasma  environment 
of  the  upper  atmosphere.  The  usual  method  of  accomplishing  time-dependent  plasma  simulations 
is  the  particle-in-cell  (PIC)  code.  These  PIC  codes  simulate  a  plasma  which,  in  reality,  may  consist 
of  between  1020  and  1025  particles  by  integrating  the  equations  of  motion  of  between,  typically, 
103  and  10s  ‘macro-particles’.  In  the  usual  explicit  electrostatic  PIC  code,  the  Lorentz  force  law  is 
used  to  ‘push’  particles  along  their  trajectories  over  some  small  time  step  (wptAt  <K  1,  where  wp  t 
is  the  electron  plasma  frequency).  The  particles  are  then  weighted  to  a  grid  to  determine  a  charge 
density  distribution  over  the  domain,  and  Poisson’s  equation  is  solved  to  determine  the  potential 
and  electric  field.  This  electric  field  is  used  to  push  the  particles  to  new  positions,  and  the  process 
repeats.  Plasma  properties  such  as  density  and  temperature  are  found  by  taking  the  appropriate 
moments  of  the  distribution  function  of  the  macro-particles. 

Numerical  stability  requirements,  however,  restrict  the  typical  explicit  PIC  code  to  time  steps  of 
wpeA  f  <  2.  Hence,  available  computer  processing  power  restricts  explicit  PIC  codes  to  simulations 
of  processes  with  frequencies  on  the  same  order  as  the  electron  plasma  frequency.  For  simulating 
low  and  medium  frequency  phenomena  (<*/  «;  wpe),  implicit  methods  must  be  used.  The  implicit 
method  uses  a  modified  form  of  Poisson’s  equation  which  results  in  an  algorithm  that  is  stable  for 
all  time  steps  and  can  therefore  be  used  to  simulate  low  frequency  processes. 

The  objective  of  the  current  research  is  to  develop  an  implicit  PIC  code  capable  of  simulat¬ 
ing  low  and  medium  frequency  plasma  phenomenon  The  code  is  to  be  developed  in  a  manner 
commiserate  with  the  ultimate  goal  of  studying  spacecraft-upper  atmosphere  interactions  Work 
towards  this  objective  has  continued  over  the  last  year  An  explicit  code  (used  as  a  baseline  model) 
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Figure  l 


B. 


Earth’s  reference  frame. 


•  Fermi-Uke  acceleration  of  ambient  O*  ions. 
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Figure  Fluid  and  kinetic  descriptions  of  an  ambient  ion  going  through  a  water  cloud 
with  a  smooth  potential  in  the  outer  region  and  sharp  potential  gradients  in  the  inner  region 


Figure  5.  Fluid  and  kinetic  descriptions  of  an  ambient  ion  going  through  a  plume  with 
smooth  potential  variations  (SI  units) 
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Figure  7.  Fluid  and  Kinet.c  descriptions  of  an  ambient  .on  going  through  a  plu 
F'vrge  potential  gradients  (SI  units). 
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and  an  implicit  code  capable  of  modelling  low  and  medium  frequency  phenomena  have  both  been 
developed  and  are  currently  functional.  Presently,  the  codes  are  doubly  periodic.  Development  of 
bounded,  explicit  and  implicit  codes  is  currently  underway.  Further  validation  of  all  the  codes  is 
also  continuing  by  comparing  the  results  of  the  codes  which  each  other  and  with  examples  from 
the  literature.  Extensions  of  the  code  to  include  more  complex  geometries  and  chemical  reactions 
are  planned  for  the  future. 

5  Low  frequency  instabilities  in  contaminant  plasma  clouds 

In  this  section,  the  abstract  of  Mr.  Mogstad’s  thesis  is  reproduced.  For  more  detail,  the  thesis 
can  be  referenced  (T.  Mogstad,  S.  M.  Thesis,  Dept,  of  Nuclear  Engineering,  MIT,  1987). 

It  has  been  observed  that  the  Space  Shuttle  in  its  ionospheric  orbit  has  a  relatively  dense  plasma 
cloud  associated  with  it.  The  Shuttle  is  releasing  water,  and  the  plasma  cloud  is  generated  when 
ambient  O'"  ionizes  the  water  through  a  charge  exchange  reaction.  Electrostatic  noise  at  frequencies 
between  0-30  kHz  has  been  recorded  in  the  Shuttle  vicinity.  We  have  performed  a  stability  analysis 
of  the  plasma  cloud  with  respect  to  the  E  x  B  instability.  We  model  the  cloud  density  distribution 
as  a  cylindrical  step  function  along  the  geomagnetic  field,  and  the  Shuttle  is  travelling  in  the 
west-east  direction.  This  modtion  with  respect  to  the  magnetic  field  generates  an  electric  field 
in  the  moving  Shuttel  frame,  which  triggers  a  destabilizing  electric  polarization  field  inside  the 
cloud.  By  applying  linearized  continuity  equations,  momentum  balance,  charge  conservation  and 
quasinuutrality  we  obtain  a  local  dispersion  realtion.  The  dispersion  relation  is  solved  analay  tically 
and  indicates  that  the  instability  growth  rate  exhibits  a  maximum  at  the  backside(0  -  180°) 
of  the  cloud,  wheras  the  oscillating  frequency  is  greatest  at  9  =  50°.  By  increasing  the  parallel 
wavenumber  kt,  the  temperature,  the  parallel  conductivity  or  the  collision  frequency,  the  instability 
is  reduced.  Increasing  the  Shuttle  velocity  makes  the  cloud  more  unstable.  In  addition  to  local 
modes,  we  consider  global  modes  and  we  find  that  exponentially  growing  eigenmodes  can  build  up 
on  the  backside  of  the  cloud. 

The  local  dispersion  relation  is  solved  numerically,  and  the  results  support  the  findings  from 
the  analytic  analysis  we  also  show  that  a  steep  density  gradient  is  inherently  unstable  in  crossed 
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E  and  B  fields  and  that  the  cloud  has  steady  state  solutions  only  when  linearized 

The  poloidal  localization  of  instability  modes  suggests  that  Shuttle  carried  sensors  operating 
inthe  very  low  frequency  regime  (i,  e.  below  100  Hz)  should  not  be  used  at  poloidal  angles  of 
9  =  ±90°.  Also,  radio  communications  at  all  frequencies  can  be  scattered  by  striations  at  the 
anisotropic  backside  of  the  cloud  where  the  instability  growth  rates  are  the  greatest,  the  growth 
rate  frequency  is  estimated  between  1-1000  Hz  and  the  oscillating  frequency  between  1-100  Hz 
The  latter  can  induce  signal  degradation  of  radio  waves  of  frequencies  less  than  100  Hz  only. 

6  Two  and  three  dimensional  fluid  theory  of  contaminant  plasma 
clouds 

Finally  in  this  report  we  include  two  published  papers  on  two  dimensional  fluid  theory  and  a 
recent  AlAA  paper  on  three  dimensional  fluid  theory. 
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l.arge  structures  in  the  ionosphere  ss ill  oulgas  or  eject  neutral  water  and  perturb  the  ambient  neutral 
environment  This  water  can  undergo  charge  exchange  with  the  ambient  oxygen  tons  and  form  a  water 
plasma  cloud  Additionally,  water  dumps  or  thruster  firings  can  create  a  water  plasma  cloud  A  simple 
model  for  the  evolution  of  a  water  plasma  cloud  about  a  large  space  structure  is  obtained  It  is  shown 
that  if  the  electron  density  around  a  large  space  structure  is  substantially  enhanced  above  the  ambient 
density  then  the  plasma  cloud  will  move  away  from  the  structure  As  the  cloud  moves  away,  it  will 
become  unstable  and  will  eventually  break  up  into  filaments  A  true  steady  state  will  exist  only  if  the  total 
electron  density  is  unperturbed  from  the  ambient  density  When  the  water  density  is  taken  to  be 
consistent  with  shuttle-based  observations,  the  cloud  is  found  to  slowly  dnft  away  on  a  time  scale  of 
many  tens  of  milliseconds  This  time  is  consistent  with  the  shuttle  observations. 


I.  IsrRoni.  vnos 

Since  the  dawn  of  shuttle  flights  to  the  ionosphere  it  has 
been  realized  that  the  ambient  ionosphere  is  strongly  per¬ 
turbed  by  the  presence  ol  a  large  body  in  space  [Green  et  til-. 
19X5)  Several  experiments  on  board  the  space  shuttle  have 
been  devoted  to  the  measurement  of  the  plasma  parameters  in 
the  vicimtv  of  the  spacecraft  [.Murphy  et  al..  1986;  Reusoner  et 
,il..  l9Sb]  It  has  been  found  that  the  shuttle  or  any  body  in 
space  affects  the  ambient  environment  through  the  interaction 
of  its  associated  contaminant  cloud  with  the  environment.  The 
contaminant  cloud  arises  from  surface  outgassing,  leaks  from 
life  support  systems,  thruster  brings,  and  water  dumps.  The 
contaminant  cloud  can  undergo  several  different  types  of  in¬ 
teractions  with  the  ambient  environment  ranging  from  physi¬ 
cal  interactions  such  as  momentum  transfer  to  chemical  inter¬ 
actions  such  as  charge  exchange  The  result  is  that  the  self- 
conststent  environment  around  a  large  body  in  space  may  be 
vignittcantly  different  from  the  ambient  environment.  Wc 
choose  tv'  put  the  emphasis  on  large  bodies  since  large  bodies 
can  have  an  associated  contaminant  cloud  which  makes  an 
fh  1 1  change  (o  the  ambient  environment.  By  an  0(1)  change 
we  mean  that  the  composition  of  the  ion  density  is  substan- 
tiallv  v hanged  For  example,  observations  around  the  shuttle 
! Lett  it  at.  1 9H < ]  indicate  that  the  observed  environment 
. tiller-  vonsiderably  Irom  the  expected  ambient  environment. 

It  i-  important  to  study  the  self-consistent  environment 
around  large  space  vehicles  such  as  the  shuttle  for  several 
reasons  fhe  maior  reason  is  to  understand  the  noise  that  will 
be  observed  on  am  sensors  carried  on  the  vehicle  f I’tekett  et 
at  I'b']  This  will  directly  a  fleet  the  utility  of  the  sensors  for 
purposes  ot  observation  Another  important  reason  is  to  build 
understanding  ol  the  space  plasma  environment  and  some  of 
the  processes  that  can  occur  in  it  Finally,  study  of  the  sclf- 
.otisisteii!  environment  will  help  in  understanding  such  haste 
Jester  -sues  as  how  to  mitigate  spacecraft  charging  in  low 
I  .o’1:  orbit 

In  this  paper  we  develop  a  simple  model  for  the  motion  of  a 
w a t et  pi. ivin. i  cloud  around  a  large  structure  in  the  tono- 
ph'-rc  In  -action  5  we  derive  the  equations  that  describe  the 
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large-scale  plasma  flow  around  the  body  In  section  3  we 
obtain  a  analytic  solution  for  the  case  where  there  is  no  elec¬ 
trostatic  modification  of  the  motional  potential  This  case  cor¬ 
responds  to  a  plasma  cloud  which  does  not  enhance  the  ambi¬ 
ent  density.  In  section  4  we  perform  a  multiple  scale  analysis 
and  obtain  the  equations  for  the  field  line  averaged  densities 
and  potential  For  the  case  where  the  cloud  density  exceeds 
the  ambient  density  we  show  that  the  equations  of  motion  do 
not  allow  steady  state  solutions  and  suggest  that  the  cloud 
may  be  unstable  to  the  growth  of  E  x  B  gradient  instabilities 
We  then  solve  the  model  equations  numerically  with  a  flux 
conserving  transport  (FCT)  numerical  method  and  show  the 
long  time  behavior  of  the  plasma  cloud.  We  obtain  the  ton 
residence  time  for  a  range  of  cloud  densities.  For  densities 
typical  of  what  is  observed  around  the  shuttle,  the  residence 
times  are  consistent  with  the  measurements.  Finally  in  section 
5  wc  conclude  with  a  discussion  of  the  significance  of  these 
results  for  measurements  from  large  space  vehicles  such  as  the 
shuttle  Although  the  equations  of  the  water  plasma  cloud  are 
generally  applicable,  we  choose  to  concentrate  on  the  dvnamic 
behavior  of  an  initial  water  plasma  'puff”  which  could  arise 
from  a  magnetoplasmadynamic  (MPD)  thruster  tiring  or  a 
plasma  source 

2  Dt  kiv  \ inis  or  t hi  Ltji  AtitiNs 
for  nu.  Plasma  Cloud 

Plasma  clouds  in  the  ionosphere  have  been  studied  both 
theoretically  and  experimentally  for  a  number  of  years  [Per¬ 
kins  et  al..  1973.  Zahuske  et  at..  1973;  .Mitchell  et  til..  1985] 
The  major  motivation  has  been  to  understand  the  dynamics  of 
such  clouds  since  they  have  been  used  for  tracing  the  magnetic 
field  lines  and  they  may  interfere  with  communications  In  this 
work  we  shall  follow  a  similar  path  to  obtain  the  equations  for 
the  plasma  and  self-consistent  electrostatic  potential  around  a 
large  body 

We  work  in  a  frame  attached  to  the  space  structure  and 
start  by  assuming  that  the  space  structure  is  emitting  neutral 
water  into  an  ambient  background  of  oxygen  ions  and  neu¬ 
trals  The  spacecraft  is  taken  to  be  moving  al  orbital  velocttv 
with  respect  to  the  ambient  environment  I  ^8  km  s  lor  low 
I  . trill  orbit)  For  an  ambient  environment  whose  neutral  den 
st t v  is  ol  the  order  of  III"  cm  '  the  mean  free  path  ol  the 
water  molecules  is  mans  kilometers  Hence  tor  length  scales  ot 
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up  to  a  few  kilometers  we  can  take  the  water  molecules  as 
freely  expanding  with  their  thermal  velocity  and  write  the 
steady  state  neutral  water  density  as 

Vo  =  Vo'V,,  IT  III 

This  assumption  about  the  neutral  density  is  reasonable  for 
low  water  densities  (water  density  much  less  than  the  ambient 
neutral  density!  but  will  have  to  be  modified  for  higher  den¬ 
sities.  However,  since  we  wish  to  concentrate  on  the  plasma 
dynamics,  we  shall  leave  more  detailed  modeling  of  the  neu¬ 
trals  for  future  work.  In  (I)  we  have  approximated  the  struc¬ 
ture  as  a  sphere  of  radius  r0  and  taken  r  as  the  radial  distance 
from  the  center  of  the  sphere.  This  approximation  will  be 
reasonable  for  any  body  shape  as  long  as  we  are  several  body 
dimensions  from  the  body  So  for  example,  the  space  shuttle, 
whose  longest  dimension  is  of  the  order  of  40  m.  can  be  rea¬ 
sonably  approximated  in  this  manner  at  distances  of  the  order 
of  100  m  and  beyond.  The  density  nHlO0  is  the  density  of  water 
at  the  surface.  The  assumption  that  the  neutral  species  is  water 
is  based  on  observations  from  the  space  shuttle,  where  the 
dominant  contaminant  species  was  observed  to  be  H  ,0  ’  [ Ca¬ 
ledonia  el  al..  1987], 

The  ionic  water  may  be  formed  by  the  charge  exchange 
reaction 

t, 

O'  +  H,0  -  O  +  HjO'  (21 

where  in  the  frame  moving  with  the  space  structure  the 
oxygen  ions  sweep  through  with  an  energy  of  up  to  5  eV.  At 
this  energy  the  reaction  rate  kCM  has  been  measured  and  is 
kel  =  6  x  W  '°  cm’s  [ Murad  and  Lai.  I986u],  and  for  O’ 
the  charge  exchange  collision  frequency  is  vti  =  nHl0ktt. 

We  assume  for  simplicity  that  the  body  is  moving  perpen¬ 
dicular  to  the  magnetic  field  with  velocity  in  a  fixed  frame  of 
V  We  define  a  wind  velocity  in  the  moving  frame  as  = 
-V  and  then  define  a  Cartesian  coordinate  system  l.x,  y,  j)  in 
the  moving  frame  by  taking  the  distance  x  to  be  along  the 
direction  of  —  V_,  the  distance  >  to  be  along  the  direction  of 
-V  x  B,  and  the  distance  r  to  be  along  the  direction  of  the 
magnetic  field  B  We  note  that  the  negative  i  direction  is 
along  the  direction  of  the  unshielded  motional  electric  field 
that  will  be  seen  in  the  far  field  from  the  moving  frame  The 
frame  is  illustrated  in  Figure  1 

We  lake  the  oxygen  and  water  ions  to  be  mainly  moving 
across  the  magnetic  field  due  to  their  large  gyroradu  while  the 
electrons  are  taken  to  dominantly  flow  along  the  magnetic 
field  due  to  their  small  gyroradu  This  assumption  on  the  ion 
motion  is  true  on  the  large  scale  where  the  ions  have  a  chance 
to  complete  their  gyroorbits  fn  (he  work  by  Caledonia  el  al 
[ 1 9x7]  the  ions  were  taken  to  be  almostly  perfectly  shielded 
from  the  motional  electric  field,  and  hence  in  the  moving 
frame  the  only  motion  available  to  (hem  was  parallel  motion 
In  this  work  the  potential  is  being  calculated  in  a  self- 
consistent  manner  so  (hai  the  ions  can  move  across  the  mag¬ 
netic  field  The  perpendicular  velocity  of  the  oxygen  ions  can 
be  obtained  from  the  steadv  state  momentum  balance  equa- 
II  on  [  k  rail  anti  I  rn  el  put  e.  1 973  ]  and  is 

vn--is,,.F  -s,,.'F  x  bill  -  N).  ‘I  ( 3 1 


and  Q0.  being  the  oxygen  cyclotron  frequency  while  \n.  H.0  is 
the  oxygen  ion-neutral  water  momentum  exchange  collision 
frequency  The  vector  F  is  given  by 

F  =  -Ic7]  eBlV  In  na.  +  Ice  B)E  -  vn.  h.o*h-o  ^o-  i5) 

In  (5).  7,  is  the  ion  temperature,  which  we  take  to  be  constant, 
and  the  first  term  is  the  force  due  to  the  ion  pressure.  E  is  the 
electric  held  seen  from  the  moving  frame  and  gives  the  electric 
force  on  the  ions  in  the  moving  frame,  and  vHO  in  the  third 
term  is  the  velocity  of  the  water,  a  term  that  expresses  the 
force  on  the  oxygen  ions  due  to  collisions  with  the  water 
neutrals  In  the  derivation  of  the  force  on  the  oxygen  ions  the 
dominant  collisions  for  these  ions  have  been  taken  to  be  col¬ 
lisions  with  the  neutral  water  molecules  irate  constant  al  2 
x  10  “  cm’  si.  Collisions  with  oxygen  neutrals  irate  constant 
a=  3  x  10  '*  cm’ s)  have  been  neglected  for  simplicity  since  we 
anticipate  little  differential  velocity  between  oxygen  neutrals 
and  ions  In  (3)  the  first  term  gives  rise  to  the  well-known 
Pedersen  conductivity  while  the  second  contains  the  E  x  B 
drift  and  the  Hall  current  We  note  that  the  Hall  parameter 
takes  very  large  values  for  parameters  typical  of  the  iono¬ 
sphere  II().  HjU  a  ’  X  10  '  s  1  and  11,,.  a  I  S  x  10;  s  '.  so 
that  a„.  ^9  x  I  O'1 1  This  suggests  that  the  dominant  motion 
for  the  oxygen  ions  will  be  the  E  x  B  drift  Hence  we  can  write 
the  continuity  equation  for  the  oxygen  ions  as 


- Vi/Jxb-Vn...  =  -  n,. .  i 
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where  the  electric  field  in  the  moving  frame  is  E  -  -  Vo  and 
the  boundary  condition  on  ( (si  is  that  for  n  •  i  the  oxygen 
density  approaches  the  ambient  density 

The  electrons  will  mainly  flow  along  the  magnetic  field  but 
will  also  have  an  E  x  B  drift  across  the  field  In  addition, 
electrons  will  be  lost  as  a  result  of  dissociative  recombination 
of  the  water  ion  This  occurs  with  a  rate  which  we  shall 

specify  later  Hence  we  can  write  the  electron  continuity  equa¬ 
tion  as 


—  -  -  V  ,;>  x  b  •  V.i  -  V  I 

,i  H  ' 


where  b  is  the  unit  vector  in  the  magnetic  field  direction  with 
x,,  being  the  Hall  parameter  civen 
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where  I  is  ihe  parallel  electron  flux  given  tn 
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and  the  potential  u,  is  defined  as 


where  the  Hall  parameter  for  the  water  ions  is 


'll,  =  0  -  -  In  <»,  <9» 

£ 

In  |9).  o  is  the  self-consistent  potential  which  is  modified  by 
the  electron  pressure  The  parallel  current  is  related  to  the 
parallel  electric  field  through  the  resistivity  It;,  ): 

i  -  -<tT  =  — l-V.,tg,)  <10* 

>1, 

where  the  resistivity  is  gisen  by 

m  > . 

<l!» 

n,a- 

In  iM  the  collision  frequency  vf  is  the  electron  momentum 
exchange  frequency  due  to  collisions  with  the  neutrals  or  the 
ions  The  physical  content  of  this  equation  is  that  the  electrons 
flow  along  the  field  lines  and  are  swept  across  it  by  the  self- 
consistent  potential  in  the  plasma  cloud  We  note  that  for  any 
steady  state  to  exist  for  this  equation  it  is  necessary  that  there 
be  a  convection  of  electrons  across  the  field  lines  This  implies 
that  the  potential  contours  must  be  open.  This  is  because 
closed  potential  contours  will  cause  the  plasma  to  rotate 
[/)ii/i  and  Whalen.  1974]  hut  will  not  lead  to  any  net  loss. 
This  implies  that  the  symmetry  breaking  part  of  the  elec¬ 
trostatic  potential  will  be  crucial  in  determining  the  motion  of 
electrons  This  part  of  the  potential  will  arise  from  the  fact 
that  the  system  is  moving  and  hence  sees  a  motional  potential. 

Once  we  have  the  equations  for  the  oxygen  ions  and  the 
electrons,  the  equation  for  the  water  ions  follows  from  the 
requirement  that  the  plasma  remain  quasi-neutral. 

V-o-  =  n,  -  no -  f>2> 

\s  an  alternative  to  using  1 6).  1 7).  and  (12)  we  could  use  1 6)  and 
i  I  2i  coupled  with  the  equation  for  the  water  ions 

gV.'-5  x  b'VnH..o-  =  ”o-v..  -  <I3) 

F  mails,  vve  need  an  equation  to  determine  the  self- 
consistent  potential  This  equation  comes  from  the  require¬ 
ment  of  quasi-neutrality  for  all  times  This  means  that  in  addi¬ 
tion  to  il2l  we  require  that  we  have  charge  conservation 
lV  j  li  where  j  is  the  current  flowing  in  the  plasma)  We 
take  the  perpendicular  current  to  be  given  by  the  perpendicu¬ 
lar  ion  llows  (equation  1 7)  and  the  analogous  equation  for  the 
pci pendicular  water  velocity)  and  the  parallel  current  to  be 
given  by  the  electron  (lows  (equation  i!0l)  and  obtain  the 
equation  [Drake  and  lluha.  l9Sb] 


In  Ous  equ.ition  lor  charge  conservation  the  perpendicular  ion 
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We  note  that  in  this  expression  all  the  E  x  B  drills  have  can¬ 
celed  out  since  thev  do  not  give  rise  to  any  net  current  The 
lirst  term  on  the  left-hand  side  of  1 14)  arises  from  the  Bede  .'sen 
current  due  to  the  electric  field,  the  second  term  dec  the 
diffusion  of  ions  down  the  density  gradients,  and  the  third  due 
to  (he  Hall  current,  while  the  fourth  is  the  electron  llow  to 
balance  the  ions.  This  equation  must  be  solved  subject  to  die 
boundary  condition  that  in  the  far  field  the  electric  field  in  the 
moving  frame  is  the  motional  electric  field.  This  requires  tiiat 
for  lr|  -  /.  we  have  V  o  =  -  V  x  B  e.  We  can  interpret  the 
first  two  terms  in  the  potential  equation  (equation  1 1 4 1 )  as 
giving  the  shielding  from  the  plasma  cloud  while  the  last  term 
is  the  charge  neutralization  by  the  parallel  flow  of  electrons 
from  far  away  in  the  ionosphere. 

3.  Analytic  Solutions  ok  thk:  Cloud 
Equations 

The  set  of  equations  for  the  potential  and  densities  cannot 
in  general  be  solved  analytically.  However,  it  is  possible  to 
obtain  a  special  solution  and  draw  some  general  conclusions 
from  the  structure  of  the  equations. 

If  the  density  of  the  cloud  is  very  low.  then  we  expect  that 
the  motional  potential  will  be  only  weakly  shielded  by  the 
cloud.  In  this  case  we  approximate  the  electric  field  every¬ 
where  by  the  boundary  condition  on  (14)  and  write  V  o  = 
-V  x  B/c  for  all  space  If  we  require  that  d>  =  if>|r  ki  e  .  we  do 
not  allow  any  parallel  electric  fields),  then  in  the  absence  of 
recombination  we  have  from  (7)  in  steady  state  that 

—  V  •  V  n,  =  0  (17) 

which  has  only  two  possible  solutions,  either  that  nt  -  njn>h„„, 
everywhere  or  that  nt  varies  in  only  one  direction  in  space, 
which  is  physically  unreasonable.  The  oxygen  equation  is 

-V  •  V  n().  -  -n0.viX  118) 

which  has  the  solution 


This  shows  the  product  of  the  water  column  density  traversed 
bv  an  O  *  ion  and  the  charge  exchange  cross  section  Hence 
the  electron  density  is  constant  everywhere  while  the  oxvgcn 
ion  density  is  depleted  and  partially  replaced  by  the  water  ion 
density  The  contours  of  constant  water  ion  density  relative  to 
the  ambient  electron  density  are  shown  in  Figure  2  for  a 
typical  case  of  an  outgassing  structure  which  is  emitting  water 
from  a  central  region  of  100  m  in  radius.  We  choose  to  look 
on  these  large  length  scales  since  there  is  evidence  from 
plasma  diagnostics  package  iPDP)  measurements  around  the 
vhuttlc  of  water  ions  out  to  a  kilometer  from  the  shuttle  (G 
Murphy,  private  communication.  I9H7).  We  see  that  most  of 
the  depletion  for  the  oxygen  density  occurs  on  the  wake  side 
of  the  structure,  and  for  this  case,  somewhat  less  than  half  of 
the  oxygen  ions  arc  converted  to  water  ions  Hence  for  this 
case  the  water  plasma  cloud  is  actually  a  depletion  relative  to 
the  ambient  electron  densitv 
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I-  i(i  2  Steady  state  water  ion  density  contours  for  a  uniform  electron  density  Water  ion  density  ts  normalized  to  the 

electron  density 


4  Dirivation  and  Ncmurical  Solution 
or  thl  Avtraoi.d  Cloud  Eolations 
Since  the  plasma  floss  along  the  magnetic  held  is  extremely 
rapid  compared  to  the  perpendicular  flow,  we  can  simplify  the 
equations  describing  the  plasma  by  doing  a  formal  multiple 
scale  analysis  Physically,  we  expect  the  cloud  to  take  a  cigar 
shape  along  the  magnetic  held  The  equations  for  the  held  line 
aseraged  density  and  held  line  aseraged  potential  can  be  ob¬ 
tained  by  averaging  the  density  and  potential  equations  where 
the  averaging  operator  is  given  by 


1  t  V 
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Hence  Irom  this  point.  V,.  VH  An..  and  't>  refer  to  the  held 
line  averaged  electron  density,  water  ion  density,  etc  The  held 
line  average  operator  annihilates  the  parallel  terms  in  the  elec¬ 
tron  densitv  and  potential  equations  and  gives  for  the  density 
equations  ifrom  ihl  and  ("’ll 

-  V  'll  X  b  •  V\„.  -  v„.i  1 20) 
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In  obtaining  these  equations  we  have  used  the  multiple  scale 
assumption  to  say  that  the  densities  and  potential  have  a 
strong  dependence  on  the  perpendicular  direction  and  a  weak 
dependence  on  the  parallel  direction  This  is  a  reasonable  as¬ 
sumption  for  the  plasma  quantities  but  cannot  be  applied  to 
the  neutrals  Therefore  in  (20l.  (21 1.  and  (221  we  have  evaluated 
the  neutral  water  density  in  the  plane  r  -  0.  This  assumption 
enables  us  to  obtain  a  closed  set  of  equations  lor  the  lield  line 
averaged  quantities  and  still  contains  the  appropriate  physics 
We  note  that  in  1 22)  we  have  neglected  the  parallel  current 
carried  bv  Alfven  waves  This  point  will  be  addressed  in  future 
work 

This  set  of  equations  describes  the  dynamic  evolution  ol  a 
water  plasma  cloud  as  seen  from  a  moving  Irame  and  with  the 
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potential  calculated  self-consistenlly  We  shall  use  these  equa¬ 
tions  to  studs  the  evolution  of  a  water  plasma  cloud  where  the 
initial  densitv  may  he  larger  than  the  ambient  electron  density 
As  pointed  out  by  Caledonia  el  al.  [19x7],  such  an  enhanced 
density  cannot  arise  self-consistently  solely  by  chemical  pro¬ 
cesses  such  as  charge  exchange  Therefore  in  those  cases  where 
the  water  plasma  density  exceeds  the  ambient  densitv  we 
assume  that  the  water  plasma  has  been  injected  at  the  initial 
time  due  to  some  mechanical  mechanism  such  as  the  firing  of 
a  plasma  source. 

This  set  of  field  line  averaged  equations  for  a  plasma  cloud 
around  a  space  structure  still  cannot  be  solved  analytically, 
but  it  can  be  shown  that  the  only  steady  state  solution  al¬ 
lowed  is  the  one  previously  discussed  where  the  electron  den¬ 
sitv  is  uniform  everywhere  The  proof  of  this  statement  is 
given  in  the  appendix  and  follows  a  similar  proof  in  work  by 
Duihiei  [  1958]  Since  in  the  far  field  the  electron  density  must 
match  the  ambient  density,  we  conclude  that  plasma  ciouds 
whose  density  exceeds  the  ambient  density  cannot  be  steady 
state  structures  Therefore  the  enhanced  densities  seen  around 
the  shuttle  [Picker r  er  al..  1985]  must  be  transient  phenomena. 
The  fact  that  a  plasma  cloud  whose  density  exceeds  the  ambi¬ 
ent  will  continuously  evolve  is  also  suggested  by  the  following 
phvsical  arguments:  Such  a  plasma  cloud  will  have  a  density 
distribution  which  will  probably  peak  toward  the  center  of  the 
cloud  The  moving  cloud  will  have  Pedersen  currents  in  it 
which  act  to  modify  the  imposed  motional  field  These  Ped¬ 
ersen  currents  arc  density  dependent  since  they  arise  from  the 
bulk  motion  of  ions  due  to  the  potential  gradients.  Therefore 
we  expect,  and  (22)  confirms  this,  that  the  highcst-density  re¬ 
gions  of  the  cloud  will  modify  the  motional  potential  most 
effectively  Hence  the  low -density  cloud  edges  will  see  a  sclf- 
consistcni  F  x  B  drift  different  than  the  higher-density  re¬ 
gions  This  leads  to  a  velocity  shear  through  the  cloud  which 
causes  ii  to  continuously  distort  Therefore  no  steady  stale  will 
exist  The  time  scale  on  which  this  distortion  will  take  place 
will  be  a  lunction  of  density  It  can  be  estimated  on  the  basis 
ol  the  arguments  given  by  Perkins  et  al  [1975]  if  we  assume 
that  the  qominant  elVcci  of  the  cloud  is  to  shield  the  potential 
It  the  cloud  has  density  n  ■>  »ilmBlct„.  then  the  electric  held  in 
the  cloud  is  E ,  =  n  where  t,,  -  V  x  B  <  =  L0e. 

from  V  x  K  =  0  we  obtain 

-  -  ""mimhi  ”  ‘/v,T„ 

The  v loud  v eloeitv  is  V  -  ,  F  x  B  B'  which  gives 
f,  -  *  1  >'  1 

i  —  -Tti  n  "  min  v  i 

\X  ,  -ee  th  it  ui  the  moving  tramc  the  cloud  will  drift  backward 
with  a  veiovitv  which  depends  on  the  densitv  and  will  drift 
sidewavs  with  a  velocitv  which  depends  both  on  the  densitv 
an,;  on  h  wv  sharp  the  dcrisilv  gradient  is  in  the  direction  of 
lb,'  ■ , i ; ; . 1 1  a  mil  I  rom  these  estimates  we  sail  see  111, it  distof 
ti. hi  whi  vein  i  -  it  a  time  'sale  ol  /  _  (  I  I n  n  1 1  w  here  /  „ 

i,  the  dciisiiv  leneth  scale  I  or  the  ease  ot  the  space  shuttle  if 
we  take  -he  densitv  length  scale  to  be  the  shuttle  length  t  ~  'll 
mi  and  :  ;sc  a  water  ion  cloud  ol  10  cm  '  in  an  ambient 
deti-itv  '  |n*  ,.m  .  then  lor  a  velocity  of  >  km  s  the  dislor- 

1 ;  o 1  'ill’.'  i  -  1 1  fi  c 


We  have  solved  the  averaged  cloud  equation.  ,r  .  aa',. 

1 20).  (211.  and  (22l|  for  a  5  km  by  '  km  square  region  using  a 
uniform  101  by  101  Cartesian  mesh  The  center  o!  the  -au.t.e 
was  fixed  on  the  structure  With  VM,0.iv.  i)  and  V(l  n  j 
given  at  some  time  we  calculate  the  potential  ftoin  i22l  and 
use  that  result  to  obtain  the  (low  held  velocitic-  Then  that 
information  is  used  in  1 201  and  ( 2 !  I  to  advance  the  ;er:  it,.  ,  by 
one  time  step 

Successive  point  over  relaxation  tSPOR)  was  nsgu  '•  . r 

the  elliptic  current  balance  equal  m  This  scheme  introduces  a 
relaxation  parameter  into  the  equation,  which  then  all.  w  ,  an 
iterative  solution  of  the  problem  ( Roaeh.  I97f. ’•  The  dc,.  va¬ 
lues  were  approximated  by  second-orde:  (irate  diflcreiices. 
and  iterations  were  carried  out  until  '.he  maximum  irroi  .a  . 
less  than  10  '  between  twv>  successive  iterations  Neumann 
boundary  conditions  were  applied  in  both  the  x  and  i  direc¬ 
tions.  We  required  the  potential  to  match  the  following 
boundary  conditions  in  the  far  field: 

V  x  B 

V  ip  - - - - 

c 

in  the  V  x  B  direction  and 

V  0  =  0 

in  the  direction  of  motion.  The  potential  of  the  central  struc¬ 
ture  was  taken  to  be  4  V,  which  is  consistent  with  measure¬ 
ments  taken  from  the  shuttle 

For  the  convective  equations  we  used  the  two-dimensional 
flux  correction  method  (FCTl  of  Zalesak  [1979]  The  high- 
order  scheme  was  a  leapfrog  trapezoidal  [Roaeh.  1976]  with 
fluxes  calculated  with  the  flux  formulae  developed  by  Zalesak 
(  1984],  while  a  donor  cell  scheme  [  Roach.  1976]  was  used  lor 
the  lower  scheme  to  complete  the  FCT  algorithm  Symmcliic 
boundary  conditions  were  applied  on  the  density  distributions, 
and  the  flux  limiter  was  applied  on  every  iteration  No  as¬ 
sumptions  were  made  about  the  symmetry  of  the  solutions  so 
that  the  computations  were  carried  out  through  the  whole 
mesh  The  oxygen  ion  density  approached  the  ambient  ion 
density  (taken  to  be  2  «  10'  cm  ')  at  the  edges  of  the  square 
The  water  ion  density  was  assumed  to  have  an  initial  Gaus¬ 
sian  profile  falling  oil  on  a  length  scale  of  500  m 

i  lr  -  U 

VH,o-  ~  Vo  I'  -  Ol  exp  (  — —  J 

where  :r,l  50  m  and  K,  500  m  This  simulates  a  large 
scale  ionic  water  injection  The  initial  water  ion  prohle  for  a 
typical  simulation  is  shown  in  Figure  where  the  central 
density  for  this  case  was  \„.,,.lx  0)  -  Hf  cm  1  la  this 
ligure  ihe  x  direction  corresponds  to  the  negative  wind  veloci¬ 
tv  direction  Thai  is.  the  neutral  wind  is  blowing  from  positive 
v  to  negative  \  in  llus  liguie  The  negative  i  direction  is  ihe 
duection  ol  the  motional  elcMric  held 

In  I  igutes  4  h  we  show  i lie  initial  distribution  ol  potential 
b.r  a  densitv  profile  as  eiven  in  Future  '  and  lor  three  initial 
i  ri  water  densities  In  I  igure  4  the  initial  water  density  was 
\„  ,,.(v  -  0|  =  1  O'*  cm  '  In  this  figure  we  see  the  shielding 
'electric  held  decreases!  that  arises  when  ihe  total  ion  densitv 
evceeds  ihe  ambient  oxvgcn  ion  densitv  and  the  cloud  behaves 
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Initial  water  ion  density  contours  for  numerical  simulations  This  corresponds  to  an  initial  puff  of  water  plasma 
Water  ion  density  is  normalized  to  the  far  field  electron  density 
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Fig  7  Water  ion  density  contours  at  0  1  s  for  an  initial  central  water  ion  density  of  106  cm  J.  Water  ion  density  is 

normalized  to  the  far  held  electron  density 


as  a  classical  dielectric.  In  Figure  5  the  central  water  density 
was  taken  to  be  VH]0.(x  =  01  =  I0fc  cm  " J  Here  we  see  exact¬ 
ly  the  opposite  tendency  from  Figure  4  The  electric  field 
inside  the  cloud  increases  rather  than  decreases  This  occurs 
because  the  density  gradient  terms  in  the  current  balance 
equation  are  large  enough  to  overcome  the  natural  polariza¬ 
tion  which  occurs  due  to  the  imposed  motional  electric  field. 
The  electric  field  increases  since  a  large  mobility  driven  elec¬ 
tric  llut  is  necessary  to  balance  the  density  gradient  driven  ion 
flux  This  tendency  is  confirmed  in  Figure  6.  where  we  give  the 
potential  distribution  for  ,V„ =  ()|  =  10"  cm  ’  The  po¬ 
tential  contours  here  can  be  well  modeled  as  due  to  a  uniform 
motional  field  and  a  positive  charge  on  the  positive  y  side  of 
the  cloud  and  a  negative  charge  on  the  negative  t  side  of  the 
cloud  This  suggests  that  the  picture  of  the  plasma  cloud 
around  large  objects  tn  the  ionosphere  as  being  a  dielectric 
shield  [Kai:  el  al..  19X4]  is  too  simplistic  a  model  The  impor¬ 
tant  effect  of  the  plasma  density  gradients  must  also  be  taken 
i n ; v *  Plenum 

\  typical  simulation  is  shown  in  Figures  7  13  for  the  initial 
water  ton  pulf  given  in  Figure  3  and  with  an  initial  density  of 
N,,  ,,  lx  -=  11)  10'’  cm  '  Since  an  enhanced  density  of  this 

magnitude  would  have  to  be  caused  by  some  means  such  as  a 
plasma  source  which  would  probably  emit  warm  electrons,  we 
took  the  recombination  rate  for  the  water  as  =  3 

■  10  cm  1  s  [  Murad  and  Lae  19866],  which  corresponds  to 
an  electron  energy  of  1  eV  The  neutral  water  is  assumed  to 
have  the  form  given  in  ill  In  Figure  7  we  show  Ihe  water  ton 
denstts  contours  at  0  1  s  after  the  start  of  the  simulation  This 
corresponds  to  16  gvropertods  which  is  enough  time  so  that  a 
substantial  F.  x  B  drift  can  occur  The  water  ton  cloud  is  seen 
lo  he  drifting  backward  and  to  he  undergoing  some  distortion 


on  the  backside  of  the  cloud.  The  central  region  of  the  cloud  is 
seen  to  be  breaking  into  two  regions  with  very  steep  gradients 
on  the  backside  of  the  central  region  The  fact  that  most  of  the 
distortion  is  concentrated  on  the  backside  of  the  cloud  can  be 
understood  by  noting  that  the  high-densitv  core  of  the  cloud 
will  drift  more  slowly  than  the  rear  fringes  of  the  cloud  so  that 
the  back  edge  of  the  cloud  will  collide  with  the  central  core, 
thus  leading  to  the  observed  distortion  Also,  around  the 
source,  tonic  water  is  still  being  formed  due  to  charge  ex¬ 
change  and  then  immediately  drifting  backward  so  that  a  verv 
steep  density  gradient  in  the  ionic  water  is  observed  near  the 
source  The  distortion  on  the  backside  of  the  cloud  may  be  the 
de  elopment  of  the  K  x  B  instability  In  barium  cloud  studies 
where  the  equations  [Perkins  er  al .  1973]  for  the  hanum 
cloud  are  very  similar  to  the  ones  used  here  the  cloud  was 
observed  to  break  up  into  filaments  both  experimentally  and 
numerically  /.abmkv  et  at.  197.1]  For  barium  clouds  it  is 
believed  that  the  filaments  may  be  associated  with  the  forma¬ 
tion  of  E  x  B  instabilities  [Drake  and  Huba.  19X6]  This  leads 
us  to  suggest  that  the  filaments  may  be  due  to  a  similar 
reason  This  important  topic  of  possible  instabilities  is  (he 
subject  of  ongoing  research  and  will  be  reported  in  a  future 
publication  This  also  suggests  that  long-wavelength,  low- 
frequency  electrostatic  noise  may  be  inextricably  linked  with 
the  large  water  densities  observed  around  large  objects  In 
Figure  X  we  show  the  associated  oxygen  density  contours  at 
the  same  time  This  region  delimits  the  area  where  new  ionic 
water  is  being  created  from  the  area  where  ionic  water  exists 
because  it  h.is  come  from  the  initial  conditions  Comparison  of 
this  figure  with  the  last  figure  suggests  that  part  of  the  break 
up  of  Ihe  central  core  may  be  due  to  the  creation  >1  ionic 
water  there  m  an  asvmmctric  manner  In  this  tigurc  we  also 
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a co  the  very  sharp  density  gradients  in  the  oxygen  density 
along  the  edges  of  the  oxygen  hole.  In  0  1  s  the  region  of 
oxygen  depletion  extends  back  behind  the  structure  a  distance 
of  approximately  0.7  km  This  suggests  that  the  center  of  the 
water  cloud  is  moving  with  a  drift  velocity  of  approximately  7 
km  s  This  is  less  than  the  velocity  at  which  an  individual  ion 
would  drift  backward  if  it  saw  only  the  motional  e.cctric  field 
and  occurs  because  the  ions  in  the  cloud  are  being  affected  by 
the  self-consistent  electric  held  and  not  the  motional  electric 
held  In  Figure  9  we  show  the  potential  contours  at  0.1  >  The 
fact  that  the  cloud  has  moved  relative  to  the  structure  is  dear¬ 
ly  Aeen  in  this  figure  sinee  the  region  where  the  potential  con¬ 
tours  are  distorted  has  moved  backward.  For  Figures  10  13 
we  show  the  time  development  of  the  water  ion  cloud  at  0.15. 
ii  3.  O  3 N  and  0  3  s  The  cloud  continues  to  develop  a  tadpolel- 
ike  shape  in  the  direction  of  cloud  motion,  and  the  central 
breakup  and  gradients  become  more  pronounced.  The  insta- 
bilitv  on  the  backside  of  the  cloud  has  entered  a  nonlinear 
state  where  distinct  lingers  have  formed  on  the  backside  of  the 
doud  This  has  also  been  seen  in  barium  cloud  releases  in  the 
ionosphere.  The  doud  distribution  over  this  period  of  time 
suggests  that  measurements  of  the  doud  density  from  the 
structure  will  see  a  density  that  is  both  time  dependent  and 
spatially  anisotropic.  This  suggests  that  density  measurements 
from  the  space  shuttle  must  be  interpreted  with  care 

The  time  that  ions  spend  in  the  vicinity  of  the  structure  is 
critical  to  determining  the  range  of  chemical  reactions  that 
they  can  undergo  This  question  is  of  substantial  interest  since 
there  is  evidence  that  chemical  reactions  can  occur  over  sur¬ 
faces  in  space  [ Green  el  al .  1985].  In  Table  I  we  give  an 
estimate  of  the  residence  time  of  a  water  ion  in  the  vicinity  of 
ihe  space  shuttle  as  a  function  of  the  initial  central  water  ion 
densitv  The  residence  time  :  wa>  defined  as  t  =  I.  i  p  where  L 
is  ihe  length  of  the  shuttle,  taken  as  50  m.  and  i  p  is  the  doud 
drift  velocity  determined  from  the  simulation  We  see  that  the 
residence  time  of  the  ions  is  in  the  millisecond  to  tens  of 
milliseconds  range  and  furthermore  the  higher-density  cases 
lose  urns  substantially  faster  than  the  low-density  cases  This 
suggests  that  measurement  of  the  decay  rate  of  a  pull  of 
plasma  will  not  give  a  unique  answer  for  the  ion  residence 
lime  ( Sasaki  cl  al.  1985]  Rather  the  ion  residence  time  is  a 
function  of  the  density,  as  we  would  expect  for  a  loss  process 
which  is  not  linear  For  ihe  low-density  case  the  residence 
time  is  vonsistent  with  measurements  from  the  shuttle  Caledon¬ 
ia  ,  t  al .  1987]  In  the  work  by  Caledonia  ei  al  [1987]  it  was 
shown  that  a  signal  of  n„.().  rtn.  of  0  I  was  consistent  with  u 
loss  trne  for  the  water  ions  of  40  ms  We  also  note  that  these 
times  .r-  consistent  with  results  from  Spacelab  3.  where  en¬ 
hanced  ionic  densities  were  no!  seen  and  substantial  shielding 
of  ihe  motional  potenf  ll  was  not  observed  (J  Rant,  private 
communication.  19861 
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5  Cost  u  skins 

We  have  formulated  a  simple  model  for  the  plasma  cloud 
around  a  large  structure  in  the  ionosphere  One  such  structure 
is  the  space  shuttle,  and  another  may  be  a  system  like  the 
space  nation.  The  model  has  been  applied  to  an  ionic  water 
puff,  and  the  dynamic  behavior  of  the  cloud  followed 

\  number  of  important  conclusions  can  be  drawn,  first  is 
that  the  cloud  possesses  no  steady  state  if  the  electron  densitv 
is  enhanced  above  the  ambient  density.  The  second  is  that  the 
effect  of  density  gradients  in  the  cloud  may  actually  enhance 
the  electric  held  that  the  cloud  sees,  so  that  the  picture  of  the 
cloud  always  acting  as  a  shield  against  the  motional  potential 
is  too  simplistic.  The  third  is  that  the  cloud  may  be  sublet  to 
electrostatic  instabilities  which  will  grow  and  saturate  with 
distortion  of  the  cloud  structur’  This  will  mean  that  elec 
trostatic  noise  and  turbulence  may  always  accompany  these 
clouds  Finally,  we  conclude  that  the  cloud  structure  will  be 
highly  anisotropic,  which  suggests  that  measurements  of  the 
cloud  density  from  the  structure  must  be  treated  with  care 
Future  work  will  concentrate  on  elucidation  of  the  elec¬ 
trostatic  instabilities  associated  with  the  cloud  and  on  resolu¬ 
tion  of  the  ion  residence  times  with  '.he  observed  ehcmistrv 
around  the  shuttle. 

These  conclusions  indicate  that  sensors  or  board  -y stems 
like  the  space  shuttle  or  a  space  station  will  have  to  cope  with 
a  background  which  is  spatially  anisotropic,  temporalis  \ ars¬ 
ing.  and  also  electrostalically  noisy.  These  issues  will  also  be 
studied  in  future  work 

Appendix  Analytic  Proof  That  Clous 

Eolations  Havi  No  Stfady  Staff 
We  consider  equations  CO),  (31).  and  (32l  with  recombi¬ 
nation  in  the  electron  equation  ignored  for  simplicitv  We 
define  two  potentials  by 
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Wtr  ass  ume  that  the  density  equations  for  the  oxygen  and 
water  possess  a  steady  >taie  and  that  the  plasma  cloud  is 
monng  vs. ? t h  nome  constant  velocity  V  With  this  assumption 
v\e  can  a  nte 

i‘nn . 
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We  substitute  i25i  and  O61  into  (20i  and  the  equivalent  equa¬ 
tion  for  the  water  ions  and  use  the  definitions  of  the  potentials 
in  1 i  and  f 24i  to  obtain 
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■a here  H  ana  G  arc  con>unt  vector*  related  to  V  and  V'PH.1). 
lit:  V'H,,.  are  arbitrary  functions  of  the  water  ion  densitv  and 
ovvecn  ion  densitv,  respeetivelv  II  we  substitute  these  ex¬ 
pressions  lor  the  potentials  in  i22i.  we  obtain 

*  T„.0.  i  --  M  •  •  N-Vnll;1).  1 3 1 ) 

where  the  sectors  '1  and  N  are  independent  of  the  densities 
We  integrate  i’ll  over  a  volume  I'  bounded  b>  a  surface  of 
lonvtant  electron  densitv  to  obtain 

|  ST-|*f<>.  -  t  ./l  -  i)  (321 

Now  it  we  muitipiv  1 3 1  >  bv  4',,.  *  4/w  integrate  over  the 
same  volume,  and  use  t  32 1  and  Green's  theorem,  we  obtain 

I  iV4'„.  -  V4'H. ,,.i:  I\  =  0  (33) 

Jl 

F his  in  v'niv  possible  tf 

lW{)  V4*H  n  .  i  =  0  1 34) 

fh>  indicates  that  the  electron  densitv  equation  must  be  of 
the  form 

V  •  Vrt  -  0  1 35) 

f  he  oniv  nonconstant  solution  to  this  equation  is  for  the  elec¬ 
tron  lieiiMtv  to  vary  :n  only  one  direction  in  space,  which  is 
phvsic.illv  unreasonable  Therefore  we  can  conclude  that  if  the 
electron  densitv  has  closed  contours  then  it  has  no  steady 
-late  and  if  there  is  a  stead>  state  then  the  electron  densitv 
nv.ot  he  the  ambient  densitv  everywhere 
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The  Motion  of  Contaminant  Water  Plasma  Clouds 
About  Large  Active  Space  Structures 

D.  E.  Hastings  and  N.  a.  Gatsonis 

Department  of  Aeronautic*  and  Astronautics,  Matiach usett*  Institute  of  Technology,  Cambridge 


Large  structures  in  low  earth  orbit  will  release  neutral  water  or  ion»  through  outgaanng,  water 
dump*  or  thruster  firing*  and  thus  perturb  the  ambient  tonoaphere.  Neutral*  and  ion*  within 
the  perturbed  environment  will  undergo  chemical  reaction*  to  form  a  contaminant  cloud  It  u 
a**umed  that  the  contaminant  cloud  con*i*t*  of  ion*  luch  a*  0  +  ,  HjO'r,  HaO+  and  neutral* 
*uch  a*  O,  H.  OH  and  HjO  A  two-dimenaional  model  for  the  motion  of  the  contaminant  cloud 
perpendicular  to  the  magnetic  field  line*  ia  developed.  Numerical  aolution  of  the  derived  model 
equation*  examine*  the  effect*  of  Alfven  wave  coupling,  neutral  water  density,  ion  temperature 
and  initial  condition*  on  the  cloud  motion.  In  low  density  neutral  water  cloud*  (<  10*  cm'1)  the 
ahielding  of  the  electric  field  ia  *mall.  In  neutral  water  cloud*  with  den»itie»  of  interest  for  shuttle 
condition*  (—  10*°  cm-3)  the  shielding  is  predicted  to  be  of  the  order  of  the  motional  electric 
field  and  the  ratio*  of  the  line  averaged  densities  of  the  ions  is  consistent  with  experimental 
data  The  effect  of  ion  temperature  is  negligible  for  the  range  of  densities  and  temperatures 
con»tdered.  It  is  predicted  that  for  symmetric  itial  condition*,  the  drifting  cloud*  will  develop 
fingerlike  instabilities  Plasma  depletion*  ar  predicted  in  the  wake  of  the  moving  structure  with 
enhancements  in  the  ram  direction.  Asymmetric  initial  conditions,  such  that  could  result  from 
a  thruster  firing,  rotate  the  plasma  cloud  which  undergoes  a  highly  asymmetric  distortion. 


l  Introduction 

With  the  advent  of  large  active  structure*  in  space  such 
as  the  Bpace  shuttle  and  the  planned  space  station  it  has 
been  realized  that  these  structures  will  significantly  perturb 
the  ambient  ionosphere  by  the  generation  of  a  contaminant 
plasma  cloud.  The  plasma  cloud  may  arise  either  by  ioniza¬ 
tion  of  gases  released  from  the  structure  due  to  sucii  mech¬ 
anisms  as  outgasaing  or  thruster  firings  or  may  be  due  to 
the  deliberate  release  of  a  plasma  cloud  by  the  firing  of  a 
piasma  thruster.  It  is  of  interest  to  understand  the  dynamic 
evolution  of  these  clouds  for  three  rea-ons.  One  reason  is 
to  give  of  information  on  the  ambient  ionosphere  and  its 
chemistry.  Another  reason  is  to  understand  the  environ¬ 
mental  impact  that  large  space  structures  will  have  on  the 
ionosphere  and  thirdly  to  determine  the  impact  that  con¬ 
taminant  clouds  may  have  on  the  operation  of  the  space 
structure  itself  especially  if  it  contains  sensitive  sensors. 

The  contaminant  plasma  cloud  about  the  space  shuttle 
nas  been  studied  experimentally  by  several  authors  Pickett 
’l  it.,  1985,  Murphy  et  at.,  1986;  Rtasoncr  et  at.,  1986j. 
Caledonia  et  at.  [1987;  showed  that  their  measurements  were 
consistent  with  an  ion  residence  time  in  the  vicinity  of  the 
shuttle  of  40  ms.  Theoretical  work  on  the  plasma  cloud 
about  the  space  shuttle  started  with  tb»  suggestion  that 
the  plasma  cloud  could  be  highly  polarized  leading  to  a  long 
residence  time  for  the  ions  in  the  shuttle  frame  Katz  et  at. 
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1984] .  A  model  for  the  unsteady  motion  of  the  plasma  cloud 
was  developed  by  Haetinge  et  at.  [1988],  where  it  was  shown 
that  for  the  densities  around  the  shuttle  the  polarization 
of  the  cloud  was  highly  dependent  upon  the  ion  density. 
Another  study  of  the  plasma  cloud  which  takes  into  account 
the  effect  of  Hall  and  Pederson  currents  has  been  undertaken 
by  Ecclcs  et  at.  [1988],  The  Hall  currents  are  shown  to  rotate 
the  polarization  field  so  that  it  is  not  antiparallel  with  the 
motional  electric  field. 

In  this  work  we  extend  the  model  of  Hastings  et  at.  '1988 
to  include  the  coupling  of  the  cloud  to  Alfven  waves  which 
can  carry  current  parallel  to  the  magnetic  field.  The  cloud 
equations  are  studied  systematically  as  a  function  of  the 
parallel  interaction  length,  neutral  density,  ion  temperature, 
and  the  initial  conditions. 


2  Development  of  the  Cloud  Equations 

We  assume  that  the  space  structure  moves  with  the  or¬ 
bital  velocity  of  V,  =  8  km  > s  perpend icu iar  to  the  magnetic 
field  lines  and  creates  a  neutral  water  cloud.  In  a  reference 
frame  attached  to  the  structure  (structure  frame  1 ,  the  ambi¬ 
ent  O'1-  ions  sweep  through  the  water  neutrals  with  energies 
up  to  3  eV.  Under  such  conditions  ionic  water  is  formed  ty 
the  charge  exchange  reaction 

CG  -  HjO  O  -  HjO*  ;  l ' 

For  simplicity  we  shall  consider  a  plasma  cloud  consisting 
of  three  components,  namely,  0  +  ,HjO*  and  HjO*  .nter- 
mixed  with  a  neutral  cloud  composed  of  O,  HjO,  OH  ar.d  H 
The  ions  and  neutraj  are  allowed  to  interact  ineiastica..\ 
through  the  following  reactions, 
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y?i) 


HjO+  +  HjO  H30+*0H  (2) 

j^recom 

h2o+  +  «  h,-S+  'oh  +  h  (3) 

j^rtcon 

HjO+'  +  e  Hs-+  'H20  +  H  (4) 

The  charge  exchange  reaction  rate  is  taken  to  be  krx  =  2.4  x 
10~9  cm3/<  [ Bolden  and  Twiddy,  1972]  and  the  bimolecular 
reaction  rate  is  taken  a*  *^'0’+  ,HjO  =  17  X  10'“  CmV» 
'Caledonia  et  al.,  1987].  The  electron  recombination  reaction 
rates  are  taken  as  *h,'o+.«  =  3  x  10~7  cm3/s  [Murad  and 
Lai,  1986]  and  fc”'o"  e  =  2.4  x  10-7  cm 3/e  [Mitchell  and 
McGowan,  1983]. 

We  work  in  a  frame  of  reference  where  the  space  structure 
is  stationary,  the  magnetic  field  B  is  in  the  1  direction  and 
the  ambient  wind  in  the  negative  x  direction  with  a  speed 
of  8  km/s  (Figure  1).  In  the  structure  frame  and  in  the 
absence  of  any  other  electric  field  within  the  plasma  pertur¬ 
bation  the  observer  will  see  the  unshielded  motional  electric 
field  E„,  =  V.  x  B  in  the  negative  y  direction  (Figure  1); 
in  that  frame  the  guiding  centers  of  the  ions  will  drift  with 
the  -V,  -elocity.  However,  within  a  plasma  density  pertur¬ 
bation  the  electric  field  and  consequently  the  drift  velocity 
will  be  detei  mined  by  the  balance  between  Pedersen  and 
Hall  currents,  ambipolar  diffusion  currents,  “pickup”  cur¬ 
rents  as  well  as  parallel  currents  carried  by  the  Alfven  wings. 

Lloyd  and  Hacrendcl,  1973;  Goertz,  1980;  Neubauer,  1980). 
Charge  buildup  due  to  ion-neutral  collisions  or  “pickup” 
processes  will  create  a  polarization  field  which  will  shield 
the  imposed  motional  electric  field.  Consequently,  the  elec¬ 
tric  field  within  the  plasma  cloud  will  be  |E|  <  |Em|  and 
the  ions  will  drift  backward  (in  the  structure  frame)  with 
a  speed  which  is  less  than  the  orbital.  In  the  case  that  the 
motional  electric  field  is  entirely  shielded,  the  plasma  cloud 
will  be  stationary  or  it  will  be  moving  with  the  structure 
velocity  in  a  fixed  frame  [ Katz  et  al.,  1984]. 

The  continuity  equation  for  the  «-th  ion  with  density  N, 
.?  given  bv 

d-~L  +  v  (U0./V0.)  =  -rt0+k€XxHto  (5) 


z 


Fig  1  Coordinate  nystem  and  computational  domain 
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where  the  right-hand  side  terms  are  the  ionic  chemical 
source  and  loss  terms  obtained  from  the  chemical  re ac t i  - n 
equations  (Equations  (l),  (2),  (3),  and  (4)). 

The  drift  velocity  of  each  ion  species  can  be  obtained 
from  the  ion  momentum  equation  ( Hastings  et  al.,  1988,  Ec- 
clei  et  al.,  1988].  For  the  t-th  ion  and  if  the  ion  inertia  Is 
neglected  (which  implies  that  the  drift  velocity  of  the  ions 
is  smaller  than  the  thermal)  the  steady  state  momentum 
equation  is 


0  =  -T.VN,  +  eN,  (e  +  B)  -  R* 

6 

—m,Nj  !/,.„( U;  -  U„)  -  m,N,  t/,., (TJ,  -  UJ  (8) 

»  i 

where  T,  is  the  temperature,  N,  is  the  density,  U,  is  the  bulk 
velocity,  U„  is  the  velocity  of  the  neutral  n  and  u,.,  „  is 

the  collision  frequency  of  momentum  transfer  between  the 
species  t  and  ]  or  n.  We  included  the  term  Rf  to  account 
for  momentum  transfer  due  to  reactive  collision  of  the  ion 
1  in  reaction  b.  In  general  for  the  ion  »  created  or  lost  in  a 
reaction  4,  the  momentum  transfer  operator  has  the  form 

R‘  =  S‘m,( U,  -  U‘)  (9) 

where  5*  is  the  rate  of  production  or  depletion  of  the  i-th 
ion  in  the  reaction  b  and  the  U‘  is  the  velocity  with  which 
the  ion  is  created  or  lost  in  the  reaction.  Clearly,  with  this 
definition  of  the  momentum  operator,  there  is  no  average 
momentum  lost  for  those  ions  which  are  depleted  in  the 
reactions  considered,  since  these  ions  are  lost  with  thei.  av¬ 
erage  velocity.  However,  for  the  ions  which  are  created  it  is 
important  to  establish  the  velocity  with  which  they  leave  the 
reactions.  In  the  charge  exchange  reaction  we  assume  that 
the  water  ions  are  created  with  the  water  neutral  velocity, 
so  =  Uh,o-  In  the  bimolecular  reaction  we  assume 

that  the  created  hydronium  ions  gain  the  bulk  velocity  of 
the  water  neutrals,  so  that  we  can  write  =  Uh,<>. 

A  complete  model  for  the  momentum  transfer  operator  in 
a  chemically  reacting  plasma  is  beyond  the  scope  of  this 
paper  [Burge n,  1969;  Li,  1966].  However,  detailed  numeri¬ 
cal  simulations  indicated  that  the  dominant  mechanism  for 
momentum  transfer  are  the  elastic  ion-neutral  collisions  and 
not  the  reactive  collisions.  We  can  write  then  the  momen¬ 
tum  transfer  operators  due  to  reactive  collisions  as  follow 


R 


f  X 

H,0* 


r;,\  -  0 

-  uHl.O 


(10) 

(11) 
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written  as 


jbimol 

U 

(12) 

jrtcnm  A 

1HjO+  ~  U 
pbimol 

(13) 

►  —  Uh,o) 

(14) 

1H,0+  -  U 

(IS) 

ion  i  then 

can  be 

We  expand  the  ion  drift  velocity  in  the  expansion  pa¬ 
rameter  v' /Q,, which  is  the  collisionality,  to  obtain  to  lowest 
order 

u.‘x  *B  (23) 

and  to  first  order 


U<‘> 


/  j/  \  c 

FI,  x  B 

f;  <  b 

\\ r  y +  b3 

e;V. 

e  A\ 

24) 


u.-A  = 


A’  x  B  A-  /  i/ 
B2  +  "  '  ~ 


b  \n,J  i  +  (v'/n,) 


A‘  X  B 

B2 


1  +  («/•/«.) 


-  l 


(16) 


where  the  vector  A'  is  given  by 

A’  =  _eV  (*  +  t)  +  ikr*  +  7kr' 


(17) 


In  (17)  the  potential  is  <j>,  the  ion  temperature  is  T,,  the  ion 
density  is  JV,.  The  neutral  force  densities  which  arise  from 
elastic  and  inelastic  ion-neutral  collisions  are 


+mH,o+  ^o+|/  Uh,o  (19) 


pHjO  + 

*  N 


(20) 


where  is"  =  1„IVh,o  and  t'b,mo1  =  H)0Nn,n.  The 

comparable  force  density  for  the  t-th  which  arises  from  ion- 
ion  collisions  is 


FJ  =  rn.N,  ^  ^U,.  (21) 


The  total  collision  frequency  for  momentum  transfer  w'  is 
defined  as 

=  y  W  „  +  *  y:  Vi.,  (22) 

r» 

where  is),.CI,00  is  the  rate  of  momentum  transfer  due  to 
chemical  reactions  where  an  ion  is  created  or  lost.  For 
O*  we  have  i/^t.no  =  °.  for  HJ°*  we  have  = 

fyV,,»/fV„,()-) is"  and  for  H,G+  we  have  is”’.”,,,,,  = 

(J*^H)<i*  1  n> 


where  the  pressure  modified  potential  0,  is 

t/>.  =  4>  +  —  In  N,  (25) 

In  using  the  collisionality  as  an  expansion  parameter  we  re¬ 
strict  ourselves  to  neutral  densities  below  10M  cm-3  and 
neglect  the  Hall  terms  which  give  rise  to  the  rotation  seen 
in  the  work  of  Ecciet  et  ad.  (1988).  However,  small  collision- 
ality  allows  us  to  keep  finite  ion  temperature  in  a  consistent 
manner,  since  the  effect  of  the  ion-ion  collisions  appears  to 
first  order  in  the  velocity  which  we  neglect  for  the  density- 
equations.  In  the  work  of  Eccles  et  al.  (1988)  the  collision¬ 
ality  is  not  taken  to  be  small  and  the  ion  temperature  is 
taken  to  be  zero.  We  are  interested  in  keeping  a  finite  ion 
temperature  since  we  wish  to  model  initial  conditions  such 
as  thruster  firings. 

The  perpendicular  ion  current  due  to  collisions  is 


j..  =  T.’N-V‘" ‘Z-br-  * 


where  the  collision  frequency  V  is 

p‘  =  •'*'  -  H ut  ■> 


(26) 

(27) 


;*• 


We  note  that  all  the  ion-ion  collision  terms  have  dropped 
out  of  the  net  current,  although  they  have  not  dropped  out 
of  the  the  individual  ion  currents  .  This  is  due  to  ion  mo¬ 
mentum  conservation. 

The  perpendicular  ion  current  due  to  the  polarization 
drift  is 


c3  1  dVd> 
v\  4rr  dt 


(28) 


where  the  Alfvfin  velocity  is 


g3 

4*  N,m, 


(29) 


An  expression  for  the  parallel  current  can  be  derived  from 
the  theory  of  Alfven  wings  [ Goertz ,  1980],  The  parallel  cur¬ 
rent  is 


Jn  = 


: — - - V  E  ^ 

4x«a  B 


(30) 


We  use  (26),  (28),  and  (30)  in  the  charge  conservation 
equation,  keep  the  largest  time  derivative  terms,  neglect  the 
parallel  electric  field  compared  to  the  perpendicular  electric 
field  and  obtain 
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ion  current  is  sufficiently  large  to  give  rise  to  a  eub-t  1  r, t •  1 
change  in.  the  motional  electric  field. 

The  density  equation  for  the  t-th  ion  is  given  l  y 


V/V, 


5,  - 


where  d/ds  is  the  derivative  along  the  Alfven  characteristic 


U±  ±  uaB/S 
s  =  - 


Va 


We  can  integrate  (31)  along  the  Alfven  characteristics  to 
obtain  the  potential  equation 


where  the  righthand  side  terms  are  the  ionic  chein,.  -  ouri  > 
and  loss  terms  (see  also  (5),  (6),  (7)).  In  this  equation  we 
neglect  the  parallel  flow  terms  relative  to  the  perpendicular 
flow  terms.  This  can  be  justified  by  noting  that  the-  parallel 
scale  length  is  expected  to  be  much  longer  than  the  perpr-r 
dicular  scale  length  and  if  the  electric  field  is  not  perfectly 
shielded  [Hastings  et  a J.  1 988] ,  then  the  perpendicular  ion 
velocity  is  much  larger  than  the  parallel  ion  velocity  which 
is  in  magnitude  of  the  order  of  the  ion  acoustic  velocity 
[Caledonia  et  al.,  1987], 

The  neutral  densities  (for  water,  oxygen,  hydroxyl  and 
atomic  hydrogen)  are  obtained  as  the  solution  of  equations 
of  the  form 

dN 

-—■  +  V  {N.V „)  =  Sn~Ln  (38) 


where  v*  is  the  Alfven  velocity  of  the  t-th  ion  species.  The 
current  carried  by  Alfven  waves  is  represented  by  the  first 
term  in  f32).  As  can  be  seen  from  the  second  term  the 
Alfven  term  is  dominant  for  low  collision  frequencies  while 
for  higher  collision  frequencies  the  Pederson  terms  are  dom¬ 
inant.  The  charge  conservation  equation  can  be  seen  to  be 
of  the  form 

V  - (cVd))  =  (33) 

for  some  dielectric  constant  £  and  some  external  charge  den¬ 
sity  p,«,.  The  dielectric  constant  has  the  form 


The  parallel  interaction  length  to  perform  the  parallel  inte¬ 
gration  over  can  be  estimated  as 

~  2rrtollH^  (35) 

where  *•, i,,,, ,i  is  the  cloud  radius  in  the  perpendicular  direc¬ 
tion  and  v±  is  the  average  perpendicular  drift  velocity  of 
the  cloud.  This  expression  is  the  distance  an  Alfven  wave 
an  travel  in  the  time  the  cloud  takes  to  convect  across  a 
m.ignrtit  field  line.  If  we  use  (35)  in  (34)  we  obtain  the 
lieiectric  constant  as 

£  ~  1  f  2rri,,u,1  —  (36) 

vi 

For  plasmas  with  low  neutral  water  density  the  collision  fre¬ 
quency  is  small  and  the  drift  velocity  is  of  the  order  of  the 
orbital  speed.  In  this  case,  £  ~  1  and  the  electric  field  in  the 
plasma  is  the  motional  electric  field  since  the  plasma  does 
not  shield  out  the  field.  Physically,  the  parallel  electron  cur¬ 
rent  carried  by  Alfven  waves  is  so  large  that  no  significant 
tteid  car.  build  up  in  the  plasma  cloud  to  shield  out  the  iro- 
posed  motional  electric  field.  For  a  plasma  where  neutrai 
water  density  is  high,  the  collision  frequency  is  large  and  the 
drift  velocity  is  small.  We  find  that  t  ut>  1  and  so  the  elec¬ 
tro  fipi  i  is  highly  shielded.  In  this  rase  the  perpendicular 


where  S„,Ln  are  the  appropriate  neutral  source  and  loss 
terms.  Consistent  with  our  low  coLlisionality  limit  we  take 
the  neutrals  to  be  ballistic  and  therefore  take  the  neutral 
velocity  as  the  initial  velocity  with  which  they  enter  the  sys¬ 
tem.  For  the  oxygen  and  hydrogen  this  is  the  orbital  speed 
in  the  negative  x  direction,  i.e  U0  h  =  -V,.  This  stream 
of  neutrals  constitutes  the  neutral  wind  in  the  structure 
frame  of  reference.  For  the  water  and  hydroxyl  we  assume 
that  their  velocity  has  a  radial  direction  from  the  moving 
structure  with  magnitude  equal  to  the  thermal  speed,  i.e., 
Uh)0»i  oh  =  \/2Tn/m„e,. 

3  Numerical  Analysis  op  the  Cloud  Equations 

The  density  and  potential  equations  have  been  solved  as 
in  Gatsonis  [1987]  and  Hastings  et  al.  [1988|.  We  briefly 
review  the  procedure.  The  equations  are  solved  for  a  5  km 
by  5  km  square  region  using  a  uniform  101  by  101  cartesian 
mesh.  The  center  of  the  computational  domain  was  fixed  on 
the  structure  which  is  assumed  to  be  circular  with  a  radius 
of  50  m.  The  magnetic  field  is  in  the  »  direction  and  the 
motional  electric  field  is  in  the  negative  y  direction  (Figure 
1).  With  /VHj0+(i,y),/VHjOa(i,y)  and  Nn+( z,y)  given 
at  some  time  we  calculate  the  potential  from  (32)  and  use 
that  result  to  obtain  the  flow  field  velocities.  Then  that 
information  is  used  in  the  ion  density  equations  ((37)),  to 
advance  the  densities  by  one  time  step.  The  neutral  density 
equations  ((38))  were  advanced  in  a  similar  manner. 

A  successive  point  over  relaxation  (SPOR)  scheme  was 
used  to  solve  the  elliptic  current  balance  equation  ((32)). 
Neumann  boundary  conditions  were  applied  in  both  the  x 
and  y  directions.  We  required  the  potential  to  match  the 
following  boundary  conditions  in  the  far  field 
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in  the  V,  x  B  direction  and 

V±<f>  =  0 

in  the  direction  of  motion.  The  potential  of  the  central 
structure  was  taken  to  be  4  Volts  which  is  consistent  with 
measurements  taken  from  the  Shuttle.  A  potential  of  this 
order  would  arise  from  spacecraft-ionosphere  sheath  effects. 

For  the  convective  equations  ((37),  (38))  we  used  the 
two-dimensional  flux  correction  method  (FCT)  of  Zalesak 
[1979|.  The  high  order  scheme  was  a  leapfrog-trapezoidal 
with  fluxes  calculated  with  the  flux  formulae  developed  by 
Zalesak  [1984]  while  a  donor-cell  scheme  was  used  for  the 
lower  scheme  to  complete  the  FCT  algorithm.  Symmet¬ 
ric  boundary  conditions  were  applied  on  the  density  distri¬ 
butions  at  the  outer  boundaries  while  solid  wall  boundary 
conditions  on  the  interior.  The  flux  limiter  was  applied  on 
every  iteration.  All  simulations  were  run  to  a  time  of  0.2  s 
(~  6  ion  gyroperiods)at  which  time  the  cloud  has  moved  a 
significant  way  towards  the  edge  of  the  computational  grid. 

Symmetric  Initial  Conditions:  An  Outgassed  Cloud 

We  assumed  that  the  water  neutrals  are  outgassed  ra¬ 
dially  from  the  moving  body  with  their  thermal  velocity. 
Then  their  initial  density  can  be  written  like 

M.,o(r)  =  Ar„,„(r  =  0)(^)3  (39) 


where  r„  =  50  m  is  the  radius  of  the  moving  structure  and 
r  is  the  distance  from  its  center.  The  central  neutral  water 
density  was  taken  as  a  variable  in  the  simulations. 

For  the  water  ions  we  assumed  that  they  follow  an  initial 
distribution  given  by 


^H,o+(r't  —  °)  — 


10  4  +  AfH,o+(*>  =  0)exp 


(40) 


The  above  formula  represents  a  Gaussian  distribution 
falling  off  at  a  distance  rL  =  500  m  on  the  top  of  a  flat 
distribution.  The  central  water  ion  density  was  taken  to  be 
,VH  (r)  =  103  cm-3.  For  the  hydronium  ions  we  assume 
an  initial  flat  distribution  with  fVHj()+(r,f  =  0)  =  10”4 
cm"3.  The  neutral  ambient  species  were  assumed  to  follow 
uniform  distributions  with  lV0(r,f  =  0)  =  8  x  10*  cm-3, 
.VH(r,  t  =  0)  =  105  cm-3,  and  Noh (t,  t  =  0)  =  10“4  cm'3. 

For  the  ion  temperatures  needed  in  the  potential  equa¬ 
tion  we  assumed  that  the  outgacsed  contaminant  cloud  is 
at  a  temperature  of  300°  K,  so  that  TH)0+  =  0.025  eV  and 
7^,0+  =  0.025  eV.  For  the  oxygen  ions  however  the  tem¬ 
perature  was  that  of  the  ambient  LEO  environment  taken 
to  be  Tn*  =  0.1  eV. 

The  ion-neutral  collision  rates  were  calculated  assuming 
a  hard-sphere  collision  model.  For  collisions  between  HjO* 
or  HjO+  and  O,  H,  OH  or  HjO  the  ion  temperature  was 
taken  to  be  0.025  eV.  For  collisions  between  O +  and  O  or 


H  we  assumed  that  T0+  =  0.1  eV,  since  most  of  the  colli¬ 
sions  will  be  between  the  ambient  0+  which  is  the  domi¬ 
nant  ion  and  the  background  neutrals  O  and  H.  However, 
for  collisions  between  ionic  oxygen  and  neutral  water  or  the 
hydroxyl  radicals  it  was  assumed  that  the  relative  energy 
w as  5  eV.  For  the  time  scales  that  we  are  interested  in,  the 
outgassing  occurs  at  a  constant  rate  giving  rise  to  a  spher¬ 
ically  expanding  neutral  water  cloud.  Hence  the  outgassed 
neutrals  will  always  see  a  streaming  0+  with  energies  al¬ 
most  equal  to  the  ram  energy,  i.e.,  5  eV.  We  anticipate  also 
that  the  OH  species  upon  their  creation  follow  the  water 
neutrals  and  interact  with  the  0+  in  a  similar  manner  like 
HjO. 

We  performed  simulations  assuming  central  neutral  wa¬ 
ter  density  of  10*  cm'1,  10®  cm'3  and  101D  cm”3  for  both 
the  limits  of  the  parallel  interaction  length  of  50  and  500 
km.  The  lower  limit  of  the  parallel  length  corresponds  to 
a  500  meter  in  radius  cloud,  moving  with  drift  velocity  of 
approximately  8  km/s  in  an  environment  where  the  Alfven 
velocity  is  about  400  km/s.  The  upper  limit  of  the  paral¬ 
lel  length,  taken  to  be  500  km,  is  the  distance  along  the 
lines  of  force  from  the  cloud  in  low  earth  equatorial  or¬ 
bit  down  to  the  iJ-layer  of  the  ionosphere;  this  is  a  sit¬ 
uation  where  the  cloud  should  experience  strong  shielding 
and,  thus  should  drift  with  much  slower  than  the  orbital 
velocities.  Figures  2,  3  and  4  show  results  for  the  three 
neutral  densities  and  for  the  maximum  parallel  interaction 
length  of  500  km. 

The  three  cases  shown  in  Figures  2-4  mark  the  transi¬ 
tion  from  a  practically  unshielded  cloud  (Figure  2)  to  a 
cloud  with  moderate  shielding  (Figure  3)  to  a  cloud  with 
substantial  shielding  (Figure  4).  This  behavior  is  depicted 
very  clearly  in  the  pictures  of  the  potential.  The  potential 
contours  start  from  almost  a  flat  distribution  for  the  low 
density  neutral  cloud  to  end  in  a  highly  shielded  distribu¬ 
tion  for  the  high  density  neutral  water  cloud  case. 

In  Figure  2  where  the  central  neutral  water  is  10*  cm"3 
the  contaminant  cloud  drifts  backward  quickly.  The  front 
side  of  the  water  ion  cloud  remains  smooth  while  its  back 
develops  fingerlike  disturbances.  New  ionic  water  is  formed 
with  maximum  of  1.39  x  103  cm”3  at  the  wake  of  the 
movi  structure.  The  hydronium  ions,  initially  flat,  de¬ 
velop  a  cloud-structure  as  well.  Hydronium  ion  formation 
is  predicted  to  be  confined  at  the  wake  region  of  the  mov¬ 
ing  structure.  The  oxygen  ion  distribution  completes  the 
above  picture,  showing  a  depletion  region  which  extends  al¬ 
most  2000  m  behind  the  moving  structure.  This  depletion  is 
confined  at  the  region  where  charge  exchange  and  recombi¬ 
nation  reactions  are  important,  i.e.,  the  region  of  new  ionic 
formation. 

In  Figure  3  the  central  neutral  water  density  is  10s  cm'3 
and  shielding  effects  are  important.  In  that  case  the  cen¬ 
tral,  dense  region  of  the  cloud  is  moving  slower  than  the 
edges,  producing  a  tadpole  shaped  formation.  In  the  back 
side  of  the  cloud  a  steep  gradient  is  formed,  since  the  back 
side  moves  more  quicly  and  catches  up  with  the  slowly  mov¬ 
ing  center.  The  front  side  however  elongates  and  the  whole 
cloud  deforms  into  a  “dumbbell’  shape.  The  front  side  of 
the  moving  cloud  remains  smooth  while  its  back  side  devel¬ 
ops  finger  like  disturbances.  Overall,  there  is  a  substantial 
formation  of  ionic  water,  with  a  maximum  of  1.5x  10*  cm'3. 
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Kig  2  Ion  density  contours  and  potential  distribution  in  a  plasma  cloud  with  symmetric  initial  conditions  and  central  neutral  water 
density  of  10*  cm"3  at  a  time  of  0  2  s  after  release  The  interaction  parallel  length  is  500  km  and  the  contaminant  H3Of  and 
HjO"  temperature  is  0  025  eV  (top-lefl)Water  ion  contours  Normalisation  is  to  the  initial  central  value  of  10 ' 3  cm  3  (top-right) 
Hydronium  ion  contours  Normalisation  is  to  the  ambient  value  of  10"  4  cm-3  (bottom-left)  Oxygen  ion  contours  Normalisation  is 
to  the  ambient  value  of2xl01cm"3  ( bottom -right )  Potential  distribution  in  volts  as  measured  in  the  frame  of  the  moving  structure 
which  is  1  "cated  at  (0.0)  km 


an  -  rder  of  magnitude  increase  to  the  initial  central  density. 
The  hydronium  ion  species  develop  a  very  similar  'dumb¬ 
bell"  structure  with  a  maximum  at  the  wake,  predicted  to  be 
at  1 ,4  <  103  cm  "3.  The  oxygen  ions  deforms  in  a  cloud  struc¬ 
ture  fuming  a  depletion  region  in  the  wake  with  a  'hole' 
which  levels  off  at  0.7  x  101  cm'3  In  the  ram  direction, 
however,  an  "enhancment"  is  predicted  due  to  the  density 


pile-up.  The  morphological  picture  of  the  cloud  deforma¬ 
tion  presented  here  is  very  similar  to  what  was  predicted 
.n  previous  work  in  the  discussion  of  the  behavior  of  dielec¬ 
tric  clouds  [ Perktru  et  of,  I973|.  The  potential  figures  show 
the  shielding  behavior  of  the  cloud  to  be  confined  at  places 
where  neutral  density  is  high  A  very  important  feature  is 
shown  in  the  wake  region  of  the  moving  structure  where  the 
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Fig  3  Ion  density  contours  and  potential  distribution  in  a  plasma  cloud  with  symmetric  initial  conditions  and  central  neutral  water 
density  of  10°  cm"3  at  a  time  of  0  2  seconds  after  release  The  interaction  parallel  length  is  500  km  and  the  contaminant  H^O *"  and 
HsO^  temperature  is  0.025  eV  (top-left)  Water  ion  contours  Normalisation  is  to  the  initial  central  value  of  IQ3  cm-3  (top-nght) 
Hydronium  ion  contours  Normalisation  is  to  the  ambient  value  of  10-4  cm-3  (bottom-left)  Oxygen  ion  contours  Normalisation  is 
to  the  ambient  value  of  2  x  105  cm-3  (bottom-right)  Potential  distribution  in  volts  as  measured  in  the  frame  of  the  moving  structure 
which  is  located  at  (0,0)  km 


potential  contours  are  shown  to  steepen  in  contrast  with 
the  ram  region  where  they  smoothly  fall  off  to  match  the 
boundary  conditions.  This  behavior  of  the  potential  is  also 
depicted  in  Figure  2  and  Figure  4.  The  explanation  lies  in 
the  equation  of  the  potential  and  the  associated  currents.  In 
the  wake  the  water  neutral  density  is  low,  compared  with 
the  ambient,  while  there  is  an  excessive  build-up  of  ionic 


species.  At  high  neutral  density  regions,  collision  frequen¬ 
cies  are  large,  the  drift  velocities  are  small  and  the  electric 
field  is  shielded.  However,  in  regions  where  collision  fre¬ 
quencies  are  low  currents  due  to  ion  density  gradients  act 
with  the  motional  electric  field  and  short-circuit  the  shield¬ 
ing  currents.  These  regions  in  turn  move  very  fast  with 
large  parallel  currents  carried  by  Alfven  waves.  A  similar 
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temperature  is  0  025  eV  (top-left)  Water  ion  contours  Normalitation  is  to  initial  central  value  of  103  cm-3  (top-right)  llydromum 
1..11  contours  Normalisation  is  to  the  ambient  value  of  10~4  cm  3  (bottom-left)  Oxygen  ion  contours  Normalisation  is  to  ambient 
value  f  2  <  103  cm-3  (bottom-right)  Potential  distribution  in  volts  as  measured  in  the  frame  of  the  moving  structure  which  is 
.  t  ate  I  at  ( 0,0)  km 


behavirr  has  been  observed  before  in  cases  where  the  ion 
lensity  was  larger  than  the  ambient  Hastings  et  a/.,  1587:. 

Figure  4  represents  a  cloud  with  central  neutral  density 
f  101"  :m~3  The  shielding  effect  is  dominant  at  the  re¬ 
gion  of  high  neutral  densities.  The  water  ion  cloud  drifts 
backward  mere  slowly  than  the  previous  cases.  The  overall 
picture  :  f  the  cioud  resembles  that  of  the  Figure  3  but  on  an 


earlier  state  of  development.  At  0.2  s  the  cloud  just  started 
to  develop  the  ‘dumbbell”  shape.  There  is  a  substantial 
formation  of  ionic  water  due  to  much  longer  residence  time 
if  oxygen  ions  with  a  maximum  at  142  x  103  cm  3  The 
Sngerlike  disturbances  are  just  about  to  start  forming  in  the 
central  regions,  due  to  the  differential  streaming  of  the  ions. 
The  hydronium  ion  also  shows  significant  formation  and  a 
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Normalized  line  average  H20+  densities  Normalized  line  average  H20+  densities 


Theta(deg)  Theta(deg) 

Fig  5  Normalized  line  average  HgO  +  densities  as  function  of  the  azimuthal  angle  for  different  neutral  water  densities  Normalization 
is  to  the  initial  [O  +  ](0,t  =  0)  =  5  x  10lu  cm-2  !,  Normalised  (HjO  +  |  density  at  the  time  of  release  A,  Normalized  i H j O ^ J  density 
at  a  time  of  0.2  s  after  the  release  for  an  initial  central  neutral  water  density  of  10*  cm-1  B,  Normalited  [HjO^]  density  at  a  time  of 
0  2s  after  the  release  for  an  initial  central  neutral  water  density  of  JO®  cm-3  C,  Normalited  (HjO't|  density  at  a  time  of  0  2  a  after 
the  release  for  an  initial  central  neutral  water  density  of  1010  cm-3  (left)  The  interaction  parallel  length  is  Li  =  50  km  (right)  The 
interaction  parallel  length  is  L\\  =  500  km 


cloud  Btucture.  The  maximum  at  the  wake  is  3.4x  104  cm-3. 
There  is  also  a  large  depletion  of  the  ionic  oxygen  with  a 
“hole"  at  0.3  x  103  cm-3. 

The  simulations  in  the  case  of  an  outgassed  cloud  were 
performed  to  study  the  effects  of  the  initial  central  water 


density,  the  parallel  interaction  length  and  the  ion  temper¬ 
ature.  For  comparison  of  this  information  we  calculated  the 
line  average  densities  as  measured  from  the  center  of  the 
moving  structure  as  function  of  the  azimuthal  angle.  For  a 
species  i*  the  line  average  density  is  defined  as 


0  60  120  ISO  240  300 

Thata(deg) 


Normalized  line  average  H30+  densities 


Fig  6  Normalised  line  avera/e  fr  jO  *  densities  as  function  of  the  aiimuthal  angle  for  different  neutral  water  densities  N'ormaliiati  n 
is  to  the  initial  [0*1(0. t  0)  -  5  x  10,n  cm  1  I,  Nor  in  allied  j  H  i  O  *  [  density  at  the  time  of  release  A.  Normalized  IliO*  '  density 
at  a  time  ofO  2  s  after  the  release  for  an  initial  central  neutral  water  density  of  10®  cm  1  B,  Normalized  :  H  -|  O  *  density  at  a  tur.e  -1 
i  2  s  after  the  release  for  an  initial  central  neutral  water  density  of  ]fj®  cm'3  C.  Normalized  'HaO“'  density  at  a  time  f  0  2s  after 
the  release  for  an  initial  central  neutral  water  density  of  liy1"  cm  3  (left)  The  interact. on  parallel  length  is  L  SO  km  trightlTle 
interaction  parallel  length  is  L.  =  500  km 
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Fig  7.  Normalized  line  average  0+  densities  as  function  of  the  acimuthal  angle  for  different  neutral  water  densities  Normalisation 
is  to  the  initial  [O  +  ](0,t  =  0)  =  5  x  1010  cm"1.  1,  Normalised  [0  +  ]  density  at  the  time  of  release.  A,  Normalised  [0  +  ]  density  at  a 
time  of  0.2  s  after  the  release  for  an  initial  central  neutral  water  density  of  10*  cm'3.  B,  Normalised  |0  +  ]  density  at  a  time  of  0.2 
s  after  the  release  for  an  initial  central  neutral  water  density  of  10s  cm-3.  C,  Normalised  |0  +  ]  density  at  a  time  of  0.2  s  after  the 
release  for  an  initial  central  neutral  water  density  of  1010  cm"3  (left)  The  interaction  parallel  length  is  L\\  =  50  km.  (right)  The 
interaction  parallel  length  is  L||  =  500  km 


(l,8,t)dl 


(41) 


where  R  is  the  radial  distance  taken  to  be  2.5  km  and  8  is  the 
azimuthal  angle  measured  from  the  direction  of  the  velocity 
vector  counterclockwise.  For  all  ions  the  line  average  density 
at  time  0.2  s  after  the  release  was  normalized  with  that  of 
the  O’"  at  the  time  of  the  release.  Using  the  definition  of 
the  average  and  the  initial  conditions  one  can  show  that 
rC )r’(8,t  =  0)  =  5  x  1010  cm-3. 

The  [HjO  +  |  densities  are  shown  initially  and  at  time  0.2 
s  for  L i,  -  50  km  in  Figure  5A  and  for  Ltl  =  500  km  in 
Figure  5B.  For  all  the  central  water  densities  considered 
there  is  ionic  water  formation  in  the  wake  of  the  moving 
body.  For  the  Nh,o(t  =  0)  =  10*  cm-3  neutral  cloud  there 
is  no  difference  in  the  levels  of  ionic  water  between  the  two 
parallel  lengths.  The  shielding  is  not  important  for  the  low- 
density  cloud.  However,  for  high-density  neutral  clouds  the 
difference  becomes  more  evident,  since  shielding  slows  down 
the  cloud,  thus  making  the  charge  exchange  reaction  more 
efficient  in  producing  ionic  water. 

The  ( H •> O |  densities  are  shown  initially  and  at  time  0.2 
s  for  L,\  =  50  km  in  Figure  6A  and  for  Ln  =  500  km  in 
Figure  6B.  In  all  cases  there  is  production  of  hydronium 
ions.  The  levels  of  hydronium  ions  increase  with  the  cen¬ 
tral  neutral  density  for  both  the  interaction  lengths.  Again, 
as  for  the  water  ions,  the  low-density  case  shows  no  differ¬ 
ence  between  the  two  interaction  lengths.  For  high  neutral 
densities  however  the  large  interaction  length  allows  more 
hydronium  ion  creation  due  to  larger  water  ion  residence 
times. 


The  [0  +  ]  densities  initially  and  at  time  0.2  s  are  shown 
in  Figures  7A  and  7B,  for  interaction  lengths  of  50  and  500 
km,  respectively.  The  0+  depicts  a  depletion  region  ex¬ 
actly  where  new  ions  are  formed,  that  is,  at  the  wake  of  the 
moving  structure.  The  depletion  increases  with  increasing 
neutral  water  density  for  both  the  parallel  lengths.  For  the 
dense  cloud  with  AfH,o(r  =  0)  =  1010  cm3  the  depletion 
extends  throughout  the  structure  while  for  the  less  dense 
clouds  it  is  confined  between  8  =  120°  and  8  =  240°.  For 
f^Hto(r  =  0)  =  1010  and  L||  =  500  km  the  ratio  of  ionic 
water  to  oxygen  and  the  ratio  of  the  hydronium  ion  to  water 
ion  is  ~  0.1  at  the  wake  ,  which  is  consistent  with  measure¬ 
ments  indicated  by  CaUdoma  cl  al.  |1987|. 

From  the  above  results  we  can  conclude  that  at  low  neu 
tral  water  densities  (<  10®  cm"3)  there  is  no  difference 
in  the  ionic  water  formation  between  the  two  interaction 
lengths.  The  potential  field  shows  no  effectively  shielding 
and  the  ions  that  are  formed  are  “trapped’  in  the  mag¬ 
netic  field.  For  high  neutral  water  densities  (>  10®  cm'3), 
shielding  occurs,  especially  for  the  large  parallel  interac¬ 
tion  lengths.  These  cases  represent  slow  moving  clouds 
where  large  parallel  current  build  up  at  the  edges  of  the 
cloud,  thus  diminishing  the  motional  electric  field.  The  re¬ 
actions,  then,  basically  charge  exchange  and  recombination, 
build  up  more  ionic  species  *  ‘eplete  the  ambient  oxygen 
ions.  However,  as  our  results  ..dicate,  the  dominant  factor 
which  determines  the  shielding  is  not  the  parallel  length  but 
the  initial  central  neutral  water  density.  In  the  high  den¬ 
sity  neutral  clouds  the  currents  due  to  collisions  are  large 
enough  to  shield  the  motional  electric  field.  In  Figure  8 
it  is  shown  the  magnitude  of  the  electric  field  for  a  clouJ 
with  ArH,o(r  0)  --  10lu  cm”3  and  i(  --  500  km,  initially 
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NORMALIZED  MAGNITUDE  OF  ELECTRIC  FIELD 


NORMALIZED  MAGNITUDE  OF  ELECTRIC  FIELD 


Fig  8  Magnitude  of  the  electric  field  as  measured  in  the  moving 
frame  initial  (top)  and  at  a  time  of  0.2  •  after  the  release  (bot¬ 
tom)  The  initial  condition*  are  symmetric  with  central  neutral 
water  density  of  1010  cm-3,  the  interaction  length  is  500  km 
and  the  contaminant  HjO+  and  HaO*  temperature  is  0  025  eV 
Normalisation  is  to  the  motional  electric  field  of  0  279  V/m. 


and  at  a  time  of  0.2  a  after  the  release.  The  magnitude 
of  the  electric  field  ha*  been  normalized  with  that  of  the 
motional  taken  to  be  |£m|  ~  0.28  V/m.  Shawhan  it  al. 
[  1 984 j  from  measurement*  on  STS-3  anticipated  reductions 
on  the  potential  of  the  order  of  the  motional.  Furthermore, 
Katz  it  al.  j  1984)  showed  that  the  ions  near  the  vicinity  of 
the  space  shuttle  could  be  almost  stationary  with  respect  to 
the  moving  frame.  As  Figure  8  shows,  the  electric  field  is 
strongly  diminished  in  the  near  space  of  the  body.  Initially, 
the  shielding  is  symmetric  in  both  directions,  a  direct  result 
of  the  symmetry  of  the  initial  ion-neutral  distributions.  At 


time  of  0.2  s,  however,  the  symmetry  of  the  electric  field 
in  the  motional  direction  no  longer  holds,  a  direct  result  of 
the  evolution  of  the  cloud;  the  increase  in  the  magnitude  of 
the  field  occurs  at  the  wake  of  the  body,  and.  as  we  already 
speculated,  this  may  be  attributed  to  the  high  ion  density 
region  which  is  located  at  the  wake.  In  general  the  electric 
field  disturbances  are  shown  to  be  both  highly  localized,  as 
anticipated  from  Shawhan  it  al.  [1984),  and  time  dependent. 
Finally  ,  all  cases  show  fingerlike  deformations  located  in  the 
backside  of  the  drifting  clouds  in  agreement  with  previous 
theoretical  and  experimental  data.  These  disturbances  may 
be  attributed  to  E  x  B  instabilities. 

To  study  the  effect  of  ion  temperature,  we  performed 
simulations  using  the  same  set  of  initial  conditions  as  the 
above  but  with  with  a  higher  ion  temperature  [T,  =  0.1 

tn 

ev  j. 

In  Figure  9  the  line  average  densities  for  water  ion  and 
hydronium  ion  are  shown  for  two  simulations  with  different 
ion  temperatures  but  with  NHao(r  =  0)  =  10°  cm"3  and 
L\\  =  500  km.  The  [HaO+]  are  shown  to  be  identical,  while 
[H30  +  ]  shows  a  slight  increase  with  increasing  ion  tempera¬ 
ture.  The  line  average  densities  for  oxygen  ion  were  identical 
for  the  two  cases.  The  magnitude  of  the  higher  ion  tem¬ 
perature  was  chosen  to  be  consistent  with  the  temperature 
expected  in  a  thruster  firing.  The  weak  dependence  on  the 


Normalized  line  average  ion  densities 


Theta  (deg) 


Fig.  9.  Normalised  line  average  HjO*  and  HiO+  densities  aa 
function  of  the  asimuthal  angle  for  different  ion  temperatures 
Normalisation  is  to  the  initial  [O  +  ](0,t  =  0)  =  5  *  I0l"  cm-2 
The  initial  conditions  are  symmetric  with  central  neutral  water 
density  of  10®  cm“3  and  the  interaction  length  is  50  km  IA, 
Normalised  (HjO^]  densities  at  the  time  *  f  release  IB,  Nor¬ 
malised  [H3O  +  ]  densities  at  the  time  of  release  A,  Normalised 
densities  at  a  time  of  0.2  a  after  the  release  The  contam¬ 
inant  HaO*  and  H3O*  temperature  is  0.025  eV  B,  Normalised 
(HsO  +  ]  densities  at  a  time  of  0  2  a  after  the  release  The  con¬ 
taminant  and  temperature  is  0  025  eV  AA,  Nor¬ 

malised  [HiO^]  densities  at  a  time  of  0  2  ■  after  the  release  The 
contaminant  H3O  +  and  H3O*  temperature  is  0  1  eV  BB  Nor¬ 
malised  (H3O'*]  densities  at  a  time  of  0  2  a  after  the  release  The 
contaminant  HjO+  and  HsO’’'  temperature  is  0  1  eV 
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ion  temperature  occurs  because  in  the  potential  equation 
the  diffusion  terms  (which  contain  the  temperature)  are  al¬ 
ways  subdominant  for  the  ion  density  ranges  explored.  The 
electric  field  is  primarily  determined  as  the  V,  x  B  field 
with  a  modification  due  to  the  dielectric  shielding.  For  ion 
densities  higher  than  the  ambient  the  behavior  is  different 
and  is  illustrated  by  Hastings  et  at.  (1988).  The  other  major 
effect  of  the  ion  temperature  is  on  the  collision  frequencies. 
However  for  the  density  ranges  considered  the  dominant  col¬ 
lision  frequency  is  that  due  to  oxygen  atoms  colliding  with 
ions.  The  relative  velocity  in  this  collision  frequency  is  de¬ 
termined  by  the  orbital  velocity  of  the  oxygen  which  is  only 
very  weakly  dependent  on  the  ion  temperature.  However, 
we  must  state  that  the  conclusions  about  the  temperature 
are  limited  from  the  range  of  temperatures  considered.  The 
newly  created  ions  within  the  plasma  cloud  could  have  tem¬ 
peratures  much  higher  than  0.1  eV.  The  possible  tempera¬ 
ture  for  the  HjO+  produced  in  the  charge  exchange  reaction 
is  possibly  leV  or  more.  There  are  large  source  terms  in  the 
ion  energy  equation  as  a  result  of  the  large  differential  mo¬ 
tion  of  the  reactants  but  this  consideration  is  beyond  the 
limits  of  the  isothermal  model  assumed  in  this  study. 


Asymmetric  Initial  Conditions:  A  Thruster  Cloud 
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Fig.  10  Asymmetric  (ellipsoidal)  initial  conditions  for  a  cen¬ 
tral  water  density  of  10°  cm-3  at  initial  time  and  for  a  parallel 
interaction  length  of  50  km  Ion  temperature  is  0.1  eV 
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In  order  to  simulate  a  directed  contaminant  release,  such 
that  could  result  from  a  thruster  firing  we  chose  asymmetric 
initial  conditions  for  the  water  ions  and  water  neutrals.  The 
neutral  water  was  modeled  as 


-■VH,i>*  (r  =  °)exP 


Mt)<>(r,  f  =  0)  =  lCT" 


where  x"  and  y"  are  the  coordinates  in  the  rotated  frame, 
il  is  the  decay  distance  in  the  x”  direction  and  y[  in 
•he  y"  direction.  The  decay  distances  were  chosen  so  that 
X[‘  ■■--  ISO  in  in  the  positive  x",  while  zl  =  450  m  in  the  neg¬ 
ative  x'  direction.  The  y^.  =  1 50  tn  in  the  y’  direction.  The 
initial  profile  described  by  the  above  is  an  asymmetric  ellip¬ 
soidal  distribution  on  the  top  of  a  flat  profile  and  is  shown 
in  Figure  10.  The  central  water  was  taken  to  be  109  cm-3. 
The  water  ions  were  following  a  similar  distribution  with 
a  central  density  of  103  cm*  .  The  ion  temperature  was 
taken  at  0.1  eV,  which  is  consistent  with  what  is  expected 
in  a  thruster  firing  and  the  parallel  interaction  length  was 
50  km. 

[n  Figure  11  we  show  the  normalized  density  distribu¬ 
tions  of  ionic  water,  hydronium  and  oxygen  ions  at  time  of 
0.2  s  after  the  release.  The  ion  density  peak  has  rotated  due 
to  the  E  ■  B  motion  of  the  ions.  In  addition,  both  densities 
•how  c  nsiderable  asymmetry,  while  hydronium  ion  has  bro¬ 
ken  up  into  two  discernible  cloud  structures.  The  ion  clouds 
lrift  backward  and  remain  smooth  in  their  front  sides  with 
instabilities  developing  in  the  back  side.  In  the  same  figure 
it  a  time  f  0  2  s  the  oxygen  ion  contours  depict  a  depletion 

t  Yl 


in  the  wake  of  the  moving  body,  while  the  potential  con¬ 
tours  shows  the  shielding  and  the  assymetry.  The  rotation 
of  the  cloud  is  confirmed  from  Figure  1°  where  we  show 
the  line  averages  for  the  ions  at  0.2  seconds  as  compared 
to  the  initial  line  averages.  The  oxygen  ion  density  is  seen 
to  undergo  a  slight  reduction  behind  the  body.  The  water 
ion  density  increases  by  about  an  order  of  magnitude  from 
its  initial  value  and  shows  a  broad  peak  mainly  behind  the 
body.  The  hydronium  ion  density  also  shows  asymmetry 
behind  the  body  with  the  line  averaged  density  having  two 
asymmetric  peaks  at  ~  180°  and  ~  210°  .  This  asymmetry 
is  due  to  the  nonlinear  coupling  of  the  density  and  potential 
fields  and  suggests  that  the  initial  conditions  for  a  water  re¬ 
lease  are  very  important  to  understanding  the  subsequent 
evolution  of  the  water  plasma  cloud.  This  also  suggests 
that  mass  spectrometer  measurements  taken  from  a  space¬ 
craft  such  a*  the  shuttle  will  see  a  substantial  dependence 
on  the  look  angle  relative  to  the  ram  direction.  Of  course 
one  reason  for  such  a  dependence  is  that  the  collection  of 
ions  is  much  more  efficient  in  the  ram  direction  than  in  other 
directions  due  to  the  fact  that  ions  from  the  ram  direction 
can  directly  enter  the  instrument.  Here  we  point  out  that 
another  reason  is  the  existence  of  nonlinear  distortion  of  the 
plasma  density  contours  which  is  a  strong  function  of  angle. 
From  Figure  12  we  see  that  typically  there  is  a  factor  of  2 
to  5  between  the  line  averaged  densities  for  angles  less  than 
x  and  for  angles  greater  than  x. 


•i  Summary  and  Conclusions 

A  model  has  been  developed  for  the  motion  of  contami 
nant  water  plasma  clouds  about  large  space  structures.  The 
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Fig  11  Ion  density  contours  and  potential  in  a  plasma  cloud  for  assymetric  ellipsoidal  initial  conditions  and  central  neutral  water 
density  of  10°  cm*3  at  a  time  of  0  2  s  after  release.  The  interaction  parallel  length  is  50  km  and  the  ion  temperature  is  0  1  eV 
(top-left)  Water  ion  density  contours  Normalisation  is  to  the  initial  central  value  of  103  cm-3  (top-nght)  Hydronium  density- 
contours  Normalisation  is  to  the  ambient  value  of  10-4  cm-3  (bottom-left)  Oxygen  ion  contours  Normalisation  is  to  ambient 
value  of  2  x  103  cm-3  (bottom-nght)  Potential  distribution  in  volts  as  meaau  ed  in  the  frame  of  the  moving  structure  which  is 
located  at  (0,0)  km 


model  includes  emission  of  Aifven  waves  from  the  cloud  to 
carry  current  along  the  magnetic  field.  The  ions  exhibit 
E  x  B  drifts  which  is  determined  by  the  current  balance 
within  the  plasma  cloud.  A  parallel  interaction  length  scale 
is  derived  which  represents  the  coupling  of  the  cloud  with 
the  ionosphere  via  parallel  currents  carried  by  the  Alfv^n 
waves  emitted  by  the  moving  cloud.  The  cloud  equations 
have  been  studied  systematically  as  a  function  of  the  paral¬ 


lel  interaction  length,  neutral  density,  ion  temperature  and 
initial  conditions.  We  conclude  that  for  low  density  neutral 
water  clouds  (<  10°  cm-3)  the  amount  of  shielding  is  small 
and  the  electric  field  is  the  motional.  In  this  case  the  cloud 
drifts  with  almost  the  orbital  speed  in  the  structure  refer¬ 
ence  frame.  However,  for  neutral  water  densities  of  interest 
for  shuttle  conditions  10 |r>  cm*1)  shielding  is  very  signif 
icant  in  the  vicinity  of  the  moving  structure  and  is  predicted 
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Fig  12  Ion  line  average  densities  as  a  function  of  azimuthal 
angle,  for  initial  and  at  time  of  0  2  s  conditions  The  interaction 
length  is  50  km,  the  ion  temperature  is  0  1  eV  and  the  initial 
neutral  water  density  is  10®  cm'3  Ionic  line  averages  are  nor- 
mahred  to  the  initial  line  average  for  ionic  oxygen  Asymmetric 
(ellipsoidal)  initial  conditions. 


to  be  of  the  order  of  the  motional  electric  field,  in  accordance 
with  experimental  data.  This  represents  a  case  in  which  the 
plasma  cloud  drifts  almost  with  the  moving  structure.  For 
the  highest  density  neutral  cloud  the  ratios  of  the  line  av¬ 
eraged  signals  is  consistent  with  experimental  data.  The 
effect  cf  ion  temperature  is  negligible  for  the  range  of  densi¬ 
ties  considered.  The  effect  of  different  (asymmetric)  initial 
conditions  is  mainly  to  allow  the  plasma  cloud  that  is  formed 
to  rotate  under  th  crossed  electric  and  magnetic  fields  and 
to  undergo  a  highly  asymmetric  distortion.  In  all  the  cases 
considered  it  is  predicted  that  the  backward  moving  clouds 
will  develop  fingerlike  instabilities  in  the  back-side  of  the 
clouds;  the  appearance  of  these  instabilities  is  in  agreement 
with  theoretical  and  experimental  data. 

Future  work  will  be  to  develop  a  three-dimensional  model 
■f  the  plasma  cloud  and  will  attempt  to  elucidate  the  low- 
frequeruy  stability  ■  >(  the  cloud 
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A  THREE-DIMENSIONAL  MODEL  FOR  AN  ARTIFICIAL  PlASMA 
CLOUD  ABOUT  A  SPACECRAFT  IN  LOW  EARTH  ORBIT 
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We  developed  a  fully  three-dimensional  model  for  an  artificial  plasma  cloud  in  the  ionosphere.  Such 
a  cloud  could  be  a  result  of  a  deliberate  release  or  contamination  about  a  spacecraft.  We  took  into 
account  finite  perpendicular  and  parallel  lengths,  all  elastic  collisions  between  the  considered  species, 
finite  temperature  effects,  variable  ionospheric  densities,  variable  neutral  winds,  ambient  electric 
fields  and  gravity.  The  velocities  of  the  charged  species  were  written  with  the  use  of  transport 
coefficient  tensors.  We  derived  an  equation  for  the  self-consisted  perturbation  potential  due  to  the 
plasma  perturbation  which  is  a  three-dimensional,  non-self  ajoint  elliptic  type  of  equation  with  very 
dissimilar  coefficients.  The  current  balance  equation  was  solved  numerically  in  order  to  obtain  the 
perturbation  potential  at  the  time  of  the  release.  We  examined  the  effects  of  the  central  density  of 
the  ions  and  the  dimensions  of  the  cloud.  The  perturbation  potential  increases  with  increasing  ion 
density.  With  densities  of  1014  m'3  the  cloud  shields  out  the  motional  electric  field.  Ion  clouds 
with  smaller  dimensions  result  in  weaker  electric  field  perturbations.  In  all  cases,  the  magnetic 
field  lines  are  not  equipotential  lines  within  a  region  bounded  by  the  density  perturbation.  The 
perturbation  along  the  magnetic  field  lines  is  stronger  for  the  case  of  denser  and  larger  clouds. 
The  effect  of  the  presence  of  neutral  was  also  investigated  through  the  use  of  a  gaussian  neutral 
cloud.  The  perturbation  potential  is  stronger  when  neutrals  are  present.  Due  to  the  screening  of 
the  motional  electric  field  the  cloud  is  expected  to  drift  with  differential  speeds  which  will  res  t It  in 
finger-like  shaping. 


Introduction 

The  presence  of  a  moving  body  in  the  ionosphere  repre¬ 
sents  a  disturbance  in  the  ambient  environment.  The  Space 
Shuttle  Orbiter  (SSO)  flights  provided  clear  evidence  of  the 
severe  .nteractions  between  a  space  vehicle  and  the  am 
blent  ionosphere.  Cangnan  and  Miller  j  1 9831 ,  Shawhan  et 
al  119841,  Picket  et  af. !  1085]  Murphy  at  u/.il986'i.  Caledonia 
et  al.  10871. 

The  study  of  the  induced  environment  about  a  space  vehi 
fie  ,s  extremely  important  as  a  design  parameter.  Especially 
f  r  :>  v-  hl'ie  which  is  expected  to  remain  in  space  for  a  1  ng 
peri'il.  .ike  the  proposed  space  station,  it  is  important  to 
unuerst and  the  characteristics  of  the  induced  environment 
if  .*  •  ru.  .at  to  be  able  to  model  this  environment  for  frag 
calculations,  prediction  of  surface  degradation  and  space- 
raft  barging,  .ontrol  of  electrostatic  noise  and  airgi  ,w. 
estimate  (  the  impact  on  the  ionosphere  and  the  scale  f 
iisi.irbance  that  it  represents. 

*  ’■ l  r  J  tin  t  e  R  i»if  wf  h  -tu-lrnt  Member  A(AA 

Ain'.  i  itf  Profe**or.  Member  AlAA 


In  this  Btudy  we  will  focus  our  attention  to  the  plasma 
cloud  which  surrounds  a  space  vehicle  in  Low  Earth  Orbit 
(LEO).  Hastings  et  at.  ;1988),  Cafjonij.;]987  .Eccles  et  al. 
jl988|,//os(inyj  and  Gatsonis  [1989'.  The  schematic  of  the 
process  that  we  study  is  as  follows:  we  consider  a  vehicle  in 
LEO  which  creates  a  neutral  cloud  around  it.  The  neutrai 
cloud  expands  while  undergoing  chemical  reactions  which 
produce  ions.  The  motion  aDd  composition  of  this  plasma 
cloud  is  then  a  result  of  complicated  chemical  and  eiectrody- 
namical  processes.  mhe  schematic  that  we  consider  is  quite 
general  and  thus  can  be  applied  to  processes  like  outgassing 
or  a  thruster  firing.  It  can  also  be  applied  to  another  class 
of  events,  those  ...f  ion  cloud  experiments.  The  purpose  of 
this  study  is  to  obtain  a  theory  and  a  model  Dr  the  motion 
of  an  artificial  high  speed  plasma  cloud  in  the  .  :.  sphere. 

A  S-Dimensional  Model  for  an  Ion-Neutral  Cloud 
in  the  Ionosphere 

In  the  present  model  we  assume  that  there  are  two  neutral 
and  two  ion  species;  s  for  the  ejected  neutral*  and  a  for 
the  ambient  ionospheric  neutrals.  Similarly,  j*  iesigr.ates 
the  contaminant  ions  and  a*  the  ambient  ,  ns  Emphasis 
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will  be  given  on  water  neutrals  and  water  ions  since  those 
two  species  have  been  identified  as  the  major  contaminants 
about  large  space  structures.  There  are  subsequently  two 
neutral  winds  which  are  taken  into  account.  The  ambient 
r.eutral  wind,  [?.  and  the  contaminant  neutral  wind,  U.. 

We  adopted  a  simple  chemistry  model.  The  ejected  neu¬ 
trals  interact  with  the  ambient  stream  in  ions  via  charge 
exchange  reaction  to  form  ions  and  produce  neutrals.  The 
charge  exchange  reaction  is 

u  ’  +  a  —  a  +  s+  (1) 

The  expansion  of  the  neutral  gas  in  a  rarefied  environment 

liv  itself  a  very  difficult  analytical  and  computational 
problem-  One  expects  the  neutral  cloud  to  pass  through 
various  flow  regimes  and  be  highly  dependent  upon  the  con¬ 
ditions  of  its  release.  Initially  we  will  use  simple  spherically 
expanding  neutral  clouds.  This  will  enable  us  to  develop 
and  validate  the  ion  cloud  model. 

Fluid  Model  for  the  Plasma  Motion 

Coordinate  System 

Given  the  latitude  of  the  release  of  the  neutral  cloud 

we  cor.stract  an  orthogonal  coordinate  system,  the  fixed  to 

■ 


field  aligned  system  •  x.y.l) 

earth  f  e  jj.  J  figure  1)  In  that  system  the  E  points  to  the 
east.  V  points  the  north  and  /?  refers  to  the  altitude. 
Tr  ;j  system  will  1  ring  the  spatial  variability  of  the  runn- 
ent  ...  n  vspheric  parameters.  The  coordinate  system  when 
.j  out.i!  ie  f  r  analytical  purposes  is  the  field  aligned  system 
Given  the  latitude  and  altitude  of  the  release  this  system 
has  its  origin  at  the  point  •  f  the  release  and  is  orthogonal. 
The  i  axis  j  parahel  to  the  local  magnetic  field  line,  th«*  y 
axis  j  parade!  to  the  E  direction  and  the  f  direction  is  the 
pel  pend icu.ar  t  -■  h  th  the  i  and  y  directions. 


Continuity  and  Momentum  Equations 

Assume  that  the  plasma  cloud  consists  of  contaminant 
ions  with  density  n,+  and  velocity  Vm,  T! «.  '  ;.s.i>  ' 

ambient  ions  is  rt.+  and  their  velocity  is  Vn .  .  Electrons  are 
shown  with  a  density  n«  and  velocity  Vr.  The  continuity 
equation  for  any  species  t  =  a+ ,  s  ,  t  is 

1h  +v  I"'17')  =  p' v  C  m 

The  terms  on  the  r.h.s  of  the  continuity  equation  come  from 
the  production  or  loss  of  the  species  due  to  chemical  reac¬ 
tions  are 

P.  +  =  0  L.v  =  -*i„r  kc,n. 

P.c  =  L.v  -  0 

P.  =  0  L.  -  0  (2) 

Ln  the  momentum  balance  we  take  into  account  electric 
fields,  pressure  forces  ,  gravity  and  elastic  collisions.  As¬ 
suming  that  the  time  scales  of  interest  are  much  larger  than 
collision  times  the  unsteady  term  can  be  neglected.  This  as¬ 
sumption  restricts  the  analysis  to  the  low  frequencies.  The 
inertia  term  VvV  can  be  neglected  as  long  as  the  thermal 
velocity  is  much  greater  than  the  drift  velocity.  The  mo¬ 
mentum  equations  can  be  written  in  the  following  form 

—  +  x  fl)  -  (p„  +»/,.  +  V  u,AV,  =  0  (4) 

TTit  J71  f 

}** 

where  the  force  for  the  species  t  is 

F,  =  q,H  +  ^  m 

VP, 

- +  m,g  (5) 

4| 

Here  m,  is  the  mass  of  the  species,  E  is  the  electric  field, 
Vt,  is  the  collision  frequency  for  momentum  transfer,  P,  is 
the  pressure  and  g  is  the  gravitational  acceleration.  The 
ambient  neutral  wind  is  denoted  by  V0  while  the  neutral 
wind  from  the  ejected  neutrals  is  V,. 

The  momentum  equations  written  in  this  form  allows  as 
to  define  the  total  collision  frequency  for  momentum  trans¬ 
fer  due  to  elastic  collisions  for  any  of  the  plasma  species 

as 

+  •'.x  .  r  .s  +  t'.r, 

v.r  -  ".t.  -  ix...  r  -r  1 /,♦, 

V*  -  IX,.  *-  is,.  »  lxr„  ,  a  (f.) 

Implicit  Solution  of  the  Momentum  Equation 

We  have  shown  that  the  momentum  equation  for  any 
species  can  be  written  in  the  form 

ixt7  -  -lit'  -  li:  ,  £  ( 7 ; 

m  m 
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Here,  V7  is  the  velocity,  m  is  the  mass,  q  is  the  charge,  v  is 
the  total  collision  frequency  and  F*  is  the  force.  Take  now 
a  coordinate  system  with  its  z  axis  aligned  to  the  magnetic 
held  direction  B.  For  any  vector  A  we  can  write  A  -  V.  -V, 
With  the  use  of  vector  identities  we  can  obtain  the  velocity 


components 


V 


— - —  Fl  +  - z—:Fl  x  B 

qB  1  +-  x3  qB  1  +-  x3 


(») 


where  0  is  the  gyrofrequency  defined  as  fl  =  qB/m  and  k  = 
^  the  ratio  of  the  gyrofrequency  to  the  total  collision 
frequency. 


1 . 

07 
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07  j 
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3‘h  0'r  j 

2. 
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Motion  on  the  Perpendicular  Plane 


The  diffusion  coefficients  of  ions  a*  associated  with  density 
gradients  of  a  '  are 


Electron  Perpendicular  Flow 

For  the  perpendicular  component  of  the  electron  velocity 
we  cafi  write 

V,j_  =  kr.F.x  +  kw.F.x  *  6  (9) 

This  expression  can  be  further  simplified  by  examining  the 
magnitude  of  the  collisionaljty  ratio  for  electrons  at  alti¬ 
tudes  between  150  -  500  km.  For  both  the  ambient  ,md 
the  H  ud  conditions  the  collisionaljty  ratio  is  very  large  i.e 
k.  >  1.  The  Pedersen  term  in  the  perpendicular  electron 
vei  iifity  can  be  neglected  since  it  will  be  much  smaller  than 
the  Hall  term.  Furthermore  since  we  expect  that  for  both 
ambient  and  cloud  conditions  ly*a  '  <£  l  we  can  neglect  all 
the  contribution  from  the  ion  and  neutral  collisions  as  well 
ls  fr  ra  gravity  so  that  the  velocity  can  be  written  finally 


p  -i 

p 

‘ 

r 

P.. 

0 

-M*h 

o. 

.  V" 

Mh 

0 

.  . 

p 

" 

o^e  -D’H, 

t  7, 
o,  J* 

D‘„.  D'r . 

1  ^ng 
»  a,  d, 
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where  ~  n  18  l^e  <,|ectron  tial1  mobility  (in 

r  \  3  ,fl  j  l).H  s  -  bdw  .  JU  T  jg  the  Hail  electron  lif 

!  tv  :  ue  t  :•  electron  density  gradients  fin  m3  'si. 


^a*r  -  [-eiltm,,i'l+,ti,,,4  *7*r(7*r  *  h) 


(12) 

[kl~  r  -  H) 

~  ka+  H  |  f  i  jCf>*r  p  (13) 

The  denominator  C7  is  given  by 
Cl  =  ’n;,m',v’,lti/;,j,(i3a,r  +  k\r  H  )<k].r  * 

*:m,,mJ.i/,,1,i/1,1,  -1  (14) 


The  diffusion  coefficients  of  ions  a*  associated  with  iens.ty 
gradients  of  s  *  are 

D°*r  =  (k].r  -  k],  H) 

(  7  T  T  “  7  »  H  )  -  rrtaet'a*.*  (  7  *  H  7  '  M  — 
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The  diffusion  coefficients  of  ions  a*  associated  with  tensity 
gradients  of  electrons  are 

07  -  (  ;m3,  m, ,  a-, .  „  .  .k,,  H  (k3*  ,.  -  k],  H) 
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7b,.  s.  ;ution  .f  the  system  of  the  ion  momentum  •■qua- 
t.  us  vie., is  the  on  velocities.  After  some  manipulator. 

,  iar  velocity  of  the  ions  V,  x  where  t  i’,i’  an 
he  written  the  form 


^ P 

-Mm 

1 

-  rn.  .  m,-  i. , .  „  .  u, {k,.  ,,  (,  f  r  ) 

-m3.  -n>;,  v..  (7-  r  -*.“♦«)( *7  r  *  7*  ♦  H  ) 


1  *-  All  ' 

l!  1  -  g,  1 


18) 


67 


The  diffusion  coefficient*  of  ions  a*  associated  with  den¬ 
sity  gradient*  of  a  r  are 

D\~  r  —  -m..m,.t'.4,4*',4»4*.4p(i,,r 

+*.rrj[ir,,!/Ci  (19) 

D*.  H  =  (*.»p  *  *•»  h  ) 

!«*.»: /cx  (20) 

The  diffusion  coefficient*  of  iona  s'1-  associated  with  denaity 
gradients  of  a+  are 

T?*t  p  =  1 — m. 4  m.  +  u.4 . 4  z*.  4 . 4  (k#*p 
(•t’-p  '  k«4H)  -  m.'.‘',-r.+  “  *.4p*.rp)l 

(*r.,)/c.  (2i) 

C*4H  =  ;-m, (fc.T  p^.4  h  i.4pi.4n)l 

(*r.»)/Cx  (22) 

The  diffusion  coefficient*  of  ion*  ar  associated  with  density 
gradients  of  electron*  are 

D*p  =  {  ~  fc.4„) 
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Finally  the  following  relationship*  hold  between  the  tiffu 
gion  coefficients  and  the  coefficients  f  r  gravity 


Ph  —  D%  H 


Parallel  Plasma  Motion 

The  parallel  velocity  for  any  species  f  is  given  in  the  im¬ 
plicit  solution  of  the  momentum  equation  in  the  form 

?.»  =  (33) 

m,i/i 

From  the  parallel  momentun  equations  we  can  write  the 
solution  for  the  parallel  velocity  of  a  species  t  —  a  r ,  a  ^ ,  c  in 
the  form 


.  rs  V,,na4  ri,  V,,n.4  V,n, 
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The  parallel  diffusion  coefficients  of  ambient  ions  due  to 
density  gradients  of  a*,  s+  and  e  are 
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(24> 

A  carefull  examination  of  the  obtained  transport  coefficient* 
reveals  very  important  and  simple  relations  among  them. 
Fiist,  a  generalized  Ftrutem  relation  hold*  between  the  mo¬ 
bility  and  the  diffuaion  coefficient* 


(35) 

C/ 1|  trig 

The  denominator  Ch  1*  given  by 

C,|  =  I'.fpl'p  -  ^.^4, 4^.4  0  4  ~  t',.  4^.4,  l/,  4 

~  l/..+  */.+  «*'«  +  'Wffirp  (36) 

The  parallel  diffuaion  coefficient*  of  the  contaminant  ions 
due  to  density  gradient*  in  ambient,  contaminant  and  elec¬ 
tron*  are 
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The  ambient  neutral  wind  transport  coefficients  relate  with 
the  iiff  .sicn  coefficient*  as  follow* 
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0.4  r  tr<t  (2'! 

c:  i.milar  exp  res*  i  .r.»  noid  between  the  contaminant  neu- 
t.-.i.  wina  transp  rt  coefficient*  and  the  diffusion  coefficient* 
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The  parallel  diffusion  coefficients  of  the  electrons  due  to 
density  gradients  in  ambient  and  contaminant  i  ■•ns  an  1, 
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A»  in  the  case  of  the  perpendicular  transport  coefficients, 
simple  relations  hold  between  the  paralel  transport  coeffi¬ 
cients  as  well.  A  generalized  Eimtcin  relation  holds  between 
the  parallel  mobility  and  the  parallel  diffusion  coefficients 

=  +  +  d'«W.  (39) 

Also,  the  following  relation  holds  between  the  diffusion  co¬ 
efficient  and  the  ambient  and  contaminant  neutral  wind  co¬ 
efficients 


W'  =  D'« 
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Finally,  the  diffusion  coefficients  and  the  gravitaional  drift 
coefficient  relate  as  follows 


3\  =  D'rn 


Here,  we  have  defined  the  total  direct  conductivities  in  the 
Pedersen,  Hall  and  parallel  directions,  taking  into  account 
all  the  plasma  species.  They  are  give  by 

«■> 

C 7  p  —  &  p  C  p 

&h  =  o  h  +■  e 

f  ■+  f  f 

crij  =  cr,,  <r,,  a.  ,'<51 

Note  also  that  there  is  no  electron  current  tensity  ,n  the 
Pedersen  direction,  since  we  assumed  that  electrons  mainly 
exhibit  drifts  in  the  Hail  direction. 

The  total  diffusion  current  density  due  to  gradients  of 
species  t  —  is  given  by 


D,r  —D,h 
D,h  Dir 
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Electric  Currents 


The  current  balance  within  the  plasma  cloud  is  the  key  to 
the  eiectrodynamical  interactions  taking  place.  The  current 
closure  between  the  plasma  cloud  and  the  ambient  iono¬ 
sphere  will  determine  the  self-consistent  potential  of  the 
cloud  and  consequently  will  determine  the  dynamical  be¬ 
havior  and  evolution  of  the  cloud.  As  a  first  approach  to 
the  problem  we  will  not  include  polarization  and  Alfven  cur¬ 
rents  in  the  model.  The  current  balance  in  this  model  will 
determined  by  direct  currents  ,  diamagnetic  (or  diffusion) 
currents,  currents  due  to  both  ambient  and  neutral  winds 
and  finally  currents  due  gravitational  drifts. 

The  current  is  given  by  direct  substitution  of  the  velocities 
of  the  piasroa  species  Eqa.  (11,  10,7?)  into  the  definition  of 
the  current  density. 

J  -  .  V,  ,  ,  +  ntq,V.±  (42) 

We  car.  write  for  the  total  current  density 

/  -  -  j*‘J  -  Jt''  -  J?"  -r  -t-  (43) 

wi.i  re  we  have  taken  into  account  the  direct  current,  the 
fi'f-js,.  r.  ,  urrer.t  due  to  density  gradients  of  the  source  ions, 
the  imr  ent  ions,  the  electrons  and  the  current  due  to  the 
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where  the  total  diffusion  coefficient  tensor  due  to  ions  *f  Is 
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The  total  diffusion  coefficient  tensor  due  to  ions  a*  ,s 
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The  total  diffusion  coefficient  tensor  due  to  electrons 
D,p  =  —  Dir  - 
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The  current  density  due  to  the  contaminant  (source)  neutral 
wmd  is 
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Fi.'.ally,  the  total  current  density  due  to  gravitational  drifts 

is 

P.P  -P'H  0 

=  e  P.H  P'P  0 

0  0  P.H 

Equation  for  the  Current  Balance 

We  need  now  to  derive  an  equation  for  the  potential 
within  the  plasma  cloud  .  We  make  here  the  assumption 
that  the  plasma  is  quasineutral,  that  is, 

+  -I -n.+  q,-,  =  n,q,  (55) 

This  assumption  implies  that  the  Debye  length  is  small  com¬ 
pared  with  the  length  scales  of  interest.  Also,  we  assume 
that  the  electric  fields  are  electrostatic,  that  is  E  =  -V<£, 
where  0  is  a  scalar  potential.  This  in  turn  implies  that 
dB/dt  =  0  which  means  that  the  self-consistent  plasma 
currents  produce  negligible  variations  in  the  magnetic  field 
B  Charge  conservation  then  implies 

V  /  =  0  (56) 

In  the  analysis  before  we  denoted  by  E  the  self-consistent 
electric  field.  We  take  now  the  electric  field  to  be 


potential  4>  It  is  important  to  recast  the  above  equation  :u  i 
form  of  a  steady-state  advection-diffusion  equation  for  the 
potential  <f>  as  follows 


where,  £  =  E„  +  From  inspection  we  see  that  all  the 
terms  which  are  not  <t>  dependent  can  be  regarded  as  source 
terms.  Note  also  that  it  is  only  the  direct  (  conductive)  cur¬ 
rent  which  depends  on  the  potential.  Thus,  we  can  consider 
the  divergence  of  all  the  other  current  densities  as  source 
terms,  in  an  non-self  ajoint  elliptic  equation  for  the  poten¬ 
tial.  This  car.  be  written  as 


Poh  d<t>  Poh  _ 

t- - —  — —  —  — - - —  +  O  —  (J 

dy  dz  dz  dy 

where  <rp,<rw,<rz  and  5  are  functions  of  z,y,c. 


(60) 


E(z,  y,  z)  =  E„  -r  f-  E„  (57) 

Here  E: ,  ,s  the  ambient  self-consistent  electric  field,  (  is  the 
self-consistent  potential  due  to  the  presence  of  the  plasma 
perturbation  and  E„  s  the  motional  electric  field  used  to 
account  for  any  change  of  reference  frame  that  we  might 
consider  The  equation  for  charge  conservation  is  now 

v  r‘"  -  v  sy  ~  v  sy  -r  v  yy  +  v  yy 

*v  y?**  +  v  ./?•*  + v  /"*  =  o  (ss) 

If  we  .assume  that  the  electric  fields  are  electrostatic  then 
the  perturbation  potential  Is  <p  such  that  t  =  -  Vo.  We  will 
concentrate  m  the  divergence  of  the  direct  current,  since  this 
s  the  term  which  will  give  the  differential  dependence  of  the 


Current  closure  -  Boundary  Conditions 

To  solve  the  equation  for  the  potential  one  must  supply 
appropriate  boundary  conditions.  In  turns  -ne  needs  to 
discuss  a  current  closure  model  .  In  the  beginning  we  wi'd 
describe  a  closure  model  which  is  applicable  to  the  high  l.iti 
tude  lonospher  outside  of  the  equatorial  region.  The  density 
distribution  of  the  ambient  ionosphere  is  shown  in  Figure  2. 
The  conductivities  calculated  from  our  model  are  shown  in 
Figure  3.  For  comparison  we  calculated  conductivities 
based  on  the  formulas  and  the  ambient  conditions  given  in 
Hatrrndel  et  o/.[l967j  We  can  see  that  the  oond uctivit^s 
derived  from  our  model  compare  well  with  those  presented 
in  Table  1,  the  differences  are  due  to  the  adopted  model 
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Figure  2:  The  ambient  ionoepheric  densities  (m  J) 


Figure  3:  The  ambient  ionospheric  conductivities  (in  S/m) 

ionosphere  and  the  models  for  collision  frequencies.  From 
Figure  3  it  is  shown  that  that  the  maximum  of  the  tranverse 
conductivities  occurs  at  ~  180  km.  Above  and  bellow  that 
altitude  both  Pedersen  and  Hall  conductivities  decrease  in 
magnitude  significantly.  The  parallel  conductivity  increases 
with  altitude  and  is  between  5  and  12  orders  higher  than 
the  tranverse  conductivities.  In  high  latitudes  the  magnetic 
field  lines  will  go  through  the  magnetosphere  before  reaching 
the  conjugate  ionosphere.  Thus,  one  can  consider  that  the 
upper  F-region  to  be  a  noncunducting  boundary  since  in  ef¬ 
fect  is  ccllisionless.  The  lower  E-region  can  be  also  regarded 
as  an  insulator.  Given  that,  one  can  bound  the  ionospheric 
plasma  between  two  insulators  iocated  at  heights  hi.  and 


hu  respective!]/. 

For  the  disturbed  ionosphere  we  assume  that  at  the  upper 
and  lower  boundaries  the  flux  will  be  mainly  in  the  longitu¬ 
dinal  direction,  and  that  the  electron  fiux  will  be  equal  to 
the  ambient  ion  fiux. 

F||.  =  rtl.v 

r„.+  =0  z  =  At,  hi;  (61) 

If  we  substitude  for  the  flux  from  equations  (11,  10,34)  we 
can  write 

(<7,11  -  <T.+  n)^J  -  ~eiD<*\\  ~  (621 

But  from  the  definition  of  the  electric  field  it  will  be 

${x,y,z)  =  4>„  -l-  *>(x,y,z)  (63) 

In  the  undisturbed  ionosphere  the  equality  of  longitudinal 
fluxes  holds  as  well  so  that  after  manipulating  the  above 
equations  one  gets 

dd> 

—  =  0  z  —  hr.,  hu  (64) 

In  other  words,  we  impose  the  condition  that  no  current  will 
flow  peprendicular  to  the  bounding  insulators.  However,  for 
the  equatorial  regions  the  closure  occurs  between  the  two 
conjugate  E-regions  of  the  ionosphere  and  the  magnetic 
field  lines  which  cross  the  perturbation  do  not  neccecerely 
pass  through  the  magnetosphere;  a  closure  model  applica¬ 
ble  to  this  case  is  under  development. 

For  the  tranverse  direction  we  require  that 

0(x,y,z)  =  O  (x,  y)  — »  oo  (65) 

Numerical  Solution  of  the  potential 

The  equation  for  the  potential  (60)  is  a  non-self  adjoint 
elliptic  equation,  three  dimensional,  with  highly  dissimilar 
coefficients.  It  is  this  inherent  difficulty  that  prohibited  the 
solution  of  the  fully  three-  dimensional  problem  in  the  past, 
and  resulted  in  the  well  known  two-  dimensional  theories. 
The  difficulties  can  be  revealed  further  if  one  recasts  the 
potential  equation  in  the  form  of  a  steady-state  advection 
diffusion  equation.  From  a  physical  and  consequently  nu¬ 
merical  point  of  view  one  has  to  deal  with  time  scales  wich 
are  very  disimilar,  in  the  the  tranverse  and  parallel  direc¬ 
tions.  One  more  source  of  difficulty  lies  in  the  N'eumroan 
condition  which  J  applied  in  the  parallel  direction,  which 
is  th p  direction  of  the  highest  conductivity,  thus  making 
convergence  very  difficult  to  achieve. 

Discrete  Potential  Equation 

We  descritized  the  potential  equation  on  a  rectangular 
gTid  with  Ax  ,Ay  and  Az  being  the  grid  sues  in  the  x,  y  and 
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z  ujrections  respectively.  In  the  discrete  space  the  indices 
refer  to  the  computational  coordinates,  aligned  with 
the  z.y  and  z  directions  and  such  that  z,  =  tAi,  y,  =  }Ay 
ana  Zk  =  kAx.  For  a  grid  point  as  shown  in  Figure 

'41  the  following  finite  difference  operator  was  used 


Figure  4:  Grid  geometry 
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(66) 


Both  -A  and  F  are  functions  of  z,y,z.  The  quantities  at  the 
interfaces  of  a  cell  is  defined  as 


Ai+i  ,.k  +  A.  ,  t 


Similar  operators  applied  in  the  j  and  k  directions.  Using 
the  above  and  expanding  in  all  three  directions  the  finite 
difference  equation  for  the  potential  <#.,>,*  at  the  grid  point 
t,;,k  can  be  written  as 


where  5,.,,*  is  the  descrete  source  term  obtained  through 
the  finite  difference  of  all  the  source  terms  in  the  potential 
equation.  Lf  we  work  on  a  grid  with  nz,ny  and  nz  grid 
points  in  the  x,y  and  z  direction  respectively,  applying  the 
above  discretization  in  all  the  interior  grid  points  we  er.d  up 
with  a  matrix  equation  of  the  form 

l-At-l  {#*}  =  (Sr)  \b'») 

The  subscript  p  refers  to  a  natural  ordering  which  orders 
the  three  dimensional  array  into  a  one  dimensional  vector 
given  by  the  transformation 

P  =  (k  -  2)(nz  -  2)(nz  -  2)  +  (;  -  2)(nz  -  2)  +  (i  -  1)  (70) 

Consequently  /  and  m  of  the  two-dimensional  array  denotes 
the  row  and  column  destination  of  the  elements  envolved 
in  the  finite  difference  equation  at  each  grid  point.  The 
matrix  [Ai„|  is  an  .V  x  N  sparse,  where  N  =  (nz  -  2)  x 
(ny  -  2)  x  (nz  -  2),  with  only  seven  nonzero  elements  in 
each  row.  In  order  to  make  the  problem  computationaiy 
tractable  we  developed  a  special  storage  system.  We  denote 
as  the  East, West, North  and  South  grid  points  as  shown  in 
Figure  4.  Then  for  each  grid  point  t,j,k  we  store  the  seven 
nonzero  components  along  with  the  pointer  denoting  the 
column  destination  in  the  uncompressed  array  A(m.  In  do¬ 
ing  so  the  matrix  is  compressed  into  an  N  x  7  matrix  result¬ 
ing  thus  in  O(N)  saving  in  storage  requirements.  Solution 
then  of  the  matrix  equation  was  obtained  by  a  conjugate 
gradient  algorithm,  after  applying  suitable  preconditioning 
on  the  original  matrix. 


The  unperturbed  electric  fields 


To  obtain  the  equilibrium  electric  fields  we  solve  equation 
(68).  The  problem  is  to  find  the  potential  such  that  the 
ambient  density  distribution  is  in  steady  state.  An  approx¬ 
imation  to  the  unperturbed  electric  field  can  be  obtained 
by  assuming  steady  state  ion  and  electron  continuity  equa¬ 
tions.  With  the  absence  of  source  and  sink  terms  for  the 
ions  and  the  electrons  it  is  required  that 

V  f.  =  V  I\  =0  (71) 


If  we  neglect  the  ambient  neutral  wind  and  gravity  and  un¬ 
der  the  assumption  that  there  is  no  shear  in  the  perpendic¬ 
ular  electric  field,  the  above  relation  becomes 


Wit  some  manipulation  and  applying  the  condition  of  equal 
fiuxes  cf  ions  and  electrons  at  the  boundaries  we  get 


D,  -  D,  Jn 
m'  -  (if)  Jz 


(73) 


72 


XV, b  numerical  solution  is  shown  in  Figure  l  i.t  c ' tr.p. r: 

w,ch  the  analytic  Boltiman  potential  derived  from  Equa¬ 
tion  ■ 


Figure  5:  Electric  field  obtained  from  the  numerical  solution 
(JO;id  line;  and  the  analytical  Boltiman  potential  (symbols] 


Neutral  Cloud  Model-IaUlal  Condition* 


In  most  cases  of  artificial  ion  clouds  a  neutral  cloud  in 
present.  It  is  beyond  the  scope  of  this  study  to  model  com¬ 
pletely  the  neutral  cloud.  Instead  we  assume  that  the  basic 
processes  is  self-diffusion  of  the  neutrals  with  losses  due  to 
charge  exchange  reaction.  We  also  assume  that  the  neutral 
can  be  d  .cribed  initially  by  a  gaussion  profile  given  by 


n,(r,  C  -  0)  =  n.„  exp  ^ 


x3  + 


exp 


(7<) 


win  re  z,y,z  are  the  distance  from  the  point  of  release  and 
n,  .  ;j  the  density  at  the  center  of  the  cloud.  With  this 
ni ...lei  the  mass  of  the  released  material  is  given  by  M,  - 
<r“  : ,r‘r  m.  For  the  ion  cloud  we  assume  initialy  a  profile 
g,-,.  •;  ‘  y  i  -.milar  gaussian  distribution. 


Numerical  Results 


We  stu.lied  numerically  the  release  of  an  ion  and  a  neu¬ 
tral  cloud  at  an  altitude  of  230  km  and  high  latitude.  We 
consulted  water  ions  and  neutrals,  since  water  was  found  to 
be  the  major  contaminant  about  space  structures.  However, 
the  analysis  could  be  easily  apllied  to  any  other  rele.ase  The 
computational  grid  was  orthogonal  with  20  x  20  x  64  grid 
points  in  the  i,  y  and  z  directions  respectively  The  upper 


ar.a  iowr:  ooundaries  .  eaten  .appropriately  •  -  :;.a:  there 
would  be  no  effect  on  the  obtained  numericai  solution.  A 
uniform  mesh  was  generated  with  Ax  =  Ay  =  3  km  and 
A z  =  4  km.  The  choice  of  almost  equal  grid  sues  prooved 
to  be  necessary  since  unequal  gTid  sizes  would  aiter  and  re¬ 
duce  the  physical  difference  between  the  tranverse  and  par¬ 
allel  directions.  The  temperature  of  the  contaminant  ions 
was  taken  to  be  =  0.025  eV.  The  ambient  ionosphere 

was  modeled  by  isothermal  oxygen  ions  with  Ta,  eV  and 
isothermal  electron*  with  T.  =  0.2  eV  .  In  all  the  simula¬ 
tion*  we  neglected  neutral  winds  and  gravity  The  results 
are  presented  in  a  reference  frame  moving  with  the  orbital 
speed  of  8km/sec,  so  that  in  this  Orbiter  frame  there  is  a 
motional  electric  field  of  |E„|  =  0.35  V/m.  The  computa¬ 
tional  reference  frame  is  shown  in  Figure  (6).  In  that  frame 


Figure  6:  Computational  grid 


positive  z  direction  points  towards  lower  altitudes,  the  posi¬ 
tive  y  direction  points  to  the  east  and  is  the  direction  of  the 
motional  electric  field,  and  the  negative  x  direction  points 
to  the  south  and  is  the  direction  of  the  orbital  motion.  In 
this  series  of  numerical  experiments  we  investigated  the  ef¬ 
fects  of  the  central  density,  the  cloud  dimensions  and  the 
presence  of  a  neutral  cloud  on  the  electrodynamic  behavior 
of  the  cloud. 


Ion  cloud 

In  Figures  ]7,8.9)  we  present  results  from  a  low  density- 
ion  cloud.  The  central  density  is  two  orders  lower  than 
the  ambient  and  thus,  corresponds  to  a  very  weak  density 
perturbation.  The  falling  distances  of  the  ion  cloud  are 
rj_  —  4.5  km  and  r,,  --  r.kn  and  the  mass  of  the  release  wa¬ 
ter  ions  is  A-fWj0*  =  0.002  kgT.  The  perturbation  potential 
at  the  plane  of  the  release  is  shown  in  Figure  7.  The  poten¬ 
tial  shows  a  dipole  structure  and  the  perturbation  electric 
field  is  directed  oposite  to  the  motional  and  tries  to  screen  it 
out.  In  the  center  of  the  cloud  the  perturbation  electric  field 
is  (  —  4.1  x  IQ"4  V/m.  Thus,  the  total  electric  field  in  the 
Orbital  reference  frame  will  be  almost  the  motional.  Note 
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that  the  maximum  and  minimum  of  potential  are  symmet¬ 
rical  since  diffusion  currents  are  very  small  and  effects  from 
Hall  currents  are  not  important.  The  y  =  0  plane  shows 
the  very  important  effects  of  the  three-dimensionality;  the 
magnetic  field  lines  are  no  longer  equipotential  lines,  an  as¬ 
sumption  present  in  all  the  previous  two-dimensional  cloud 
models.  The  non-equipotential  region  seems  to  be  confined 
within  the  cloud  boundaries.  Beyond  a  radius  of  ~  9  km 
in  the  perpendicular  plane  and  ~  12  km  in  the  parallel  di- 
-ection  the  magnetic  fild  lines  become  equipotential  again. 

:  .  weak  perturbation  can  be  attributed  to  the  very  slight 
aerification  of  the  conductivities  as  it  is  shown  in  Figures 
8,9). 

As  the  central  density  increases  the  potential  perturbs 
ti  n  .ncreases.  In  Figure  10  the  central  density  is  of  the 
rder  of  the  ambient  and  the  perturbed  electric  field  at  the 
.enter  is  now  c  =  4.1  x  10_1  V/m.  The  equipotentiality 
of  the  B  lines  holds  outside  of  the  cloud  boundaries.  The 
conductivities  are  altered  in  the  cloud  region  as  it  is  shown 
in  Figures  (11,12).  Both  the  Pedersen  and  the  parallel  con¬ 
ductivities  increase,  while  the  effect  of  the  ion  cloud  on  the 
Hail  conductivity  is  not  so  evident. 

In  Figures  (13,14,15)  results  are  shown  for  a  very  dense 
.loud  with  central  ion  density  of  10IS  m-J  and  Mmo+  = 
15  kgr.  The  perturbation  electric  field  is  now  0.349  V/m, 
which  shows  that  the  motional  electric  field  is  practicaly 
shielded  at  the  cloud  center.  There  is  also  a  rotation  in  the 
perturbation  electric  field  due  to  the  Hall  currents.  The  Hall 
currents  however  are  flowing  to  the  Em  x  B  direction  which 
implies  that  the  ion  Hall  conductivity  dominates  over  the 
electron.  As  Figure  14  shows  the  Hall  conductivity  becomes 
negative  (with  our  definition  of  conductivities  it  is  positive 
when  is  dominated  by  the  electron  conductivity  and  would 
imply  Hall  currents  in  the  B  x  E„  direction).  In  Figure  15 
the  Pedersen  conductivity  shows  an  increase  of  almost  three 
order  of  magnitudes  and  thus  ,  once  expects  a  significant 
perturbation  to  propagate  along  the  B  field.  The  y  =  0 
plane  plot  in  Figure  13  shows  that  the  non-  equipotential 
region  extends  now  to  almost  18  km  along  the  magnetic 
field. 

We  investigated  also  the  effects  of  the  cloud  dimensions 
on  the  electrodynamic  behavior  of  the  cloud.  In  Figure  16 
the  perpendicular  and  parallel  lengths  are  ri.  =  3  km  and 
n,  =  4  km  and  the  perturbation  field  becomes  e  =  0.341 
V/m  at  the  cloud  center.  The  rotation  of  the  electric  fieid 
denotes  Hail  currents  in  the  Em  x  B  direction.  In  Figure  17 
the  Hail  conductivity  is  shown  to  be  negative  in  the  outer 
regi on  of  the  cloud  but  remains  always  smaller  than  the 
Pedersen.  As  the  cloud  dimensions  decrease  (Figure  19)  to 
r  p4,p  1  5  km  and  r =  2  km,  the  perturbation  electric 
fieid  at  the  center  of  the  cloud  drops  to  0.327  V/rn.  The 
iipoie  structure  of  the  potential  prevails  again  but  in  this 
smail  radius  there  is  an  asymmetry  in  the  center  caused 
!.y  the  diffusion  currents,  which  are  comparable  there  with 
•.fie  direct  currents.  The  non-equipotential  region  has  also 


become  smaller  and  now  extends  to  ~  9  km  along  the  B 
lines,  and  to  a  radius  of  ~  3  km  in  the  perpenicular  plane. 

Ion  and  Neutral  Clouds 

We  studied  the  effects  of  the  presence  of  contaminant  neu¬ 
trals  along  with  the  ions.  We  considered  a  neutral  cloud 
with  central  density  of  1C17  m~3  and  falling  off  distances 
of  rp,rf  —  3  km  and  r jj  =  4  km.  This  corresponds  tc  a 
release  of  15  kgr  of  H?0+  and  1500  kgr  of  H?0.  Results 
are  shown  in  Figures  (22,23,24).  The  perturbation  electric 
fieid  is  c  =  0.3496  at  the  center  which  is  larger  than  the 
case  where  only  the  ion  cloud  is  present.  Collisions  with 
the  contaminant  neutrals  are  very  frequent  arid  the  o . 
ductivites  are  much  larger  than  the  case  where  only  ions 
are  present  as  we  see  by  comparing  Figures  17  and  23.  The 
effects  of  Hall  currents  are  also  less  visible  and  the  Hail  con¬ 
ductivity  is  dominated  by  the  electrons.  The  maximum  and 
minimum  of  the  perturbation  potential  are  also  larger.  The 
region  however  of  nonequipotential  B  lines  seems  to  extend 
at  ~  18  km  along  the  B  lines,  although  the  disturbance  itself 
appears  to  be  stronger  (compare  Figures  22  and  16). 

Conclusion  -  Future  work 

We  developed  a  fully  tree-dimensionnal  model  for  an  ar¬ 
tificial  pjlasma  cloud  in  the  ionosphere.  Such  a  cloud  could 
be  a  result  of  a  deliberate  release  of  neutrals  or  ions  or  a  re¬ 
sult  of  contamination  about  a  spacecraft.  The  neutral  cloud 
consists  of  the  ejected  neutrals  s  and  the  ambient  neutrals  a. 
The  ejected  neutrals  undergo  charge  exchange  reaction  with 
the  ambient  ions  to  create  contaminant  ions.  It  is  assumed 
that  the  cloud  has  finite  perpendicular  and  parallel  lengths. 
We  take  into  account  all  elastic  collisions  between  the  con¬ 
sidered  species  and  include  finite  temperature  effects.  We 
model  the  ionosphere  and  take  into  account  variable  densi¬ 
ties  for  ambient  neutrals  and  ions,  altitude  dependent  ambi¬ 
ent  neutral  winds,  ambient  electric  fields  and  gravity.  In  the 
momentum  balance  of  the  charged  species  we  include  elec¬ 
tric  fields,  forces  due  to  elastic  collisions,  pressure  forces  and 
gravity.  The  plasma  species  are  assumed  to  be  isothermal. 
An  analytic  solution  to  the  system  of  momentum  equations 
yields  the  velocities  of  the  charged  species  in  both  the  per¬ 
pendicular  and  parallel  plane.  The  velocities  were  written 
with  the  use  of  transport  coefficient  tensors.  These  trans¬ 
port  coefficients  account  for  the  mobility,  the  diffusion  due 
to  density  gradient  of  all  the  charged  species,  gravity  and 
ambient  and  contaminant  neutral  winds.  The  transport  co¬ 
efficients  depend  drastically  on  the  direction  and  are  greatly 
modified  by  the  presence  of  the  plasma  cloud.  The  values 
for  the  undisturbed  ionosphere  calculated  from  our  model 
com  pared  with  simplified  models  and  found  to  be  in  good 
agreement.  It  is  found  that  generalired  Einstein  relations 
hold  between  the  mobility  and  the  diffusion  transport  co¬ 
efficients.  Simple  relations  hold  between  the  diffusion  anil 
the  rest  of  the  transport  coefficients. 
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The  current  balance  between  the  cloud  and  the  ambient 
plasma  determine*  the  dynamical  behavior  of  the  plasma 
cloud.  It  is  assumed  that  the  plasma  is  quasineutral  and 
th3t  the  electric  fields  are  electrostatic.  On  applying  charge 
conservation  an  equation  for  the  self-consisted  potential  due 
to  the  plasma  perturbation  is  derived.  The  current  balance 
includes  conductive  currents,  diffusion  currents  due  to  den¬ 
sity  gradients  in  the  ambient  and  contaminant  ions  as  well 
as  electrons,  currents  due  to  the  ambient  and  contaminant 
neutral  winds  and  gravitational  drift  currents.  The  coupling 
o(  the  plasma  cloud  to  the  ambient  plasma  is  examined  with 
the  use  of  a  closure  model  in  which  insulators  are  placed  in 
the  upper  and  lower  boundaries  and  coupling  is  obtained 
via  classical  currents  only.  The  equilibrium  electric  field 
were  obtained  and  compared  with  simple  analytical  solu¬ 
tions.  The  current  balance  equation  was  solved  numerically 
in  order  to  obtain  the  perturbation  potential  at  the  time  of 
the  release.  We  examined  the  effects  of  the  central  density  o! 
the  ions  and  the  dimensiona  of  the  cloud.  The  perturbation 
potential  increases  with  increasing  ion  density.  With  densi¬ 
ties  of  101&  m-s  the  cloud  shields  out  the  motional  electric 
field.  Ion  cloud*  with  smaller  dimensions  result  in  weaker 
electric  field  perturbations.  In  all  cases,  the  magnetic  field 
lines  were  not  equipotential  lines  within  a  region  bounded 
by  the  density  perturbation.  The  perturbation  along  the  B 
lines  was  stronger  for  the  case  of  denser  and  larger  clouds. 
The  effect  of  the  presence  of  neutral  was  also  investigated 
A  gaussian  neutral  cloud  was  considered  with  the  same 
dimensions  as  the  ion  cloud.  The  perturbation  potential 
was  stronger  than  the  case  where  only  ions  were  present. 
The  electric  fields  will  determine  the  velocities  of  the  ion 
species.  It  is  expected  that  for  there  will  be  a  differential 
speed  within  the  cloud,  which  will  result  in  the  finger-like 
shaping  of  the  cloud.  Currenlty  a  numerical  model  for  the 
solution  of  the  three-dimensional  equations  for  the  densities 
is  under  development.  This  will  allow  to  follow  in  time  the 
motion  of  the  release. 
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